EAH: Sen, rie 
ATRVES i! JESKE 


THE BELL SYSTEM 
TECHNICAL JOURNAL 


DEVOTED TO THE SCIENTIFIC AND ENGINEERING ASPECTS 
OF ELECTRICAL COMMUNICATION 


Microwave Repeater Research H. T. Friis 183 


The Measurement of Delay Distortion in Microwave 
Repeaters D. H. Ring 247 


Frequency Shift Telegraphy—Radio and Wire Applications 
J. R. Davey and A. L. Matte 265 


Reflections from Circular Bends in Rectangular Wave 
Guides—Matrix Theory................. S. O. Rice 305 


The Approximate Solution of Linear Differential Equations 
Marion C. Gray and S. A. Schelkunoff 350 


Potential Coefficients for Ground Return Circuits 
W. Howard Wise 365 


Abstracts of Technical Articles by Bell System Authors... 372 


Contributors to this Issue 


AMERICAN TELEPHONE AND TELEGRAPH COMPANY 
, NEW YORK 





50¢ per copy $1.50 per Year 





THE BELLSYSTEM TECHNICAL JOURNAL 
Published quarterly by the 


American Telephone and Telegraph Company ~ 
195 Broadway, New York, N.Y. 


See Bre Gee See Bee Hoc GeeGesSe Geog 


EDITORS 
R. W. King J. O. Perrine 


EDITORIAL BOARD 


W. H. Harrison O. E. Buckley 
O. B. Blackwell M. J. Kelly 
H. S. Osborne A. B. Clark 
J. J. Pilliod F. J. Feely 


BoB 01 HeeG org eGo Gee hoacQagerk 


SUBSCRIPTIONS 


Subscriptions are accepted at $1.50 per year. Single copies are 50 cents each, 
The foreign postage is 35 cents per year or 9 cents per copy. 


Copyright, 1948 
American Telephone and Telegraph Company 


PRINTED IN U.S.A, 


PROPERTY OF. 
STEVEN JESKE 


The Bell System Technical Journal 


Vol. XXVIT April, 1948 No.2 





Microwave Repeater Research 
By H. T. FRIUS 
INTRODUCTION 


T WAS some 80 years ago that Maxwell and Hertz demonstrated that _ 
free space is a good transmission medium for electromagnetic waves. 
Since this fundamental contribution, the radio art has advanced tremen- 
dously and a decade ago it had progressed to the point where it was possible to 
construct equipment suitable for quantitative propagation studies of micro- 
waves. Such studies were made and they indicated that normal propagation 
over “line-of-sight” paths of signals of 10 to 20 centimeters wavelength 
was characterized by free space attenuation and freedom from atmospheric 
interference. These results, together with the facts that in this wavelength 
range wide bands of frequencies are available and it is possible to design 
small antennas having high directivity, encouraged us to start more com- 
prehensive research work on microwave repeater circuits. This paper gives 
the present status of the work which was interrupted by our war efforts and 
resumed at the end of the war with the construction of an experimental 
New York-Boston system as an initial objective. 

The first section will describe our propagation studies. It will be followed 
by sections on repeater circuit planning, antennas, radio frequency channel 
filters, the construction and testing of the repeater amplifier, and a conclud- 
ing section on the whole repeater. 


I. PRopAGATION STUDIES* 


That portion of the radio frequency spectrum represented by. wave- 
lengths shorter than about five meters has long been considered as the proper 
domain for point-to-point communication links, local broadcasting, and 
mobile radio communication. Since these ultra-short waves are not re- 
flected by the ionosphere, their effective range is not much greater than the 
horizon distance and it therefore becomes possible for a number of stations, 
properly separated, to operate in the same frequency band; for the same 
reason, atmospheric interference is not an important factor. in this wave- 

* This section was prepared. by A. B. Crawford who, with W. M. Sharpless, is at present 
engaged in microwave propagation studies. 
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length range. Also, as the wavelength decreases, it becomes possible to 
construct antennas large in comparison with the wavelength so that high 
antenna gains are obtained and the corresponding directivity further re- 
duces the interference areas. 

Since about 1930, with the exception of the war years, we have conducted 
fundamental studies in radio propagation, taking advantage of advances in 
the art to extend the wavelength range from about four meters (ultra-short 
wave region) in the beginning to 1.25 centimeters (microwave region) at 
the present time. A considerable portion of the effort of those engaged in 
propagation studies has, of necessity, been devoted to the development of 
measurement techniques and reliable measuring apparatus. The present 
discussion, however, will be concerned with the results of experiments rather 
than with a description of the apparatus and methods.- Most of these 
results have been described in the literature; the following is a review in- 
tended to show the development of the background leading to the present 
field trial of a microwave repeater circuit. 

The object in making propagation studies has been to evaluate and to 
understand the effects of the terrain and of the lower atmosphere upon the 
transmission of ultra-short-wave and microwave signals. The evaluation 
is usually obtained by amassing sufficient data on a particular transmission 
experiment so that a statistical analysis can be made. Efforts to under- 
stand the transmission phenomena usually take the form of experiments 
involving specially designed apparatus. These experiments are varied 
from time to time as information is obtained or as it becomes desirable to 
check the validity of such theories as may be devised. The hope is always 
present that an understanding of the phenomena may suggest a means for 
reducing the transmission difficulties. 

The absence of ionospheric reflections at these frequencies suggested at the 
start that propagation studies would probably be concerned mainly with 
phenomena familiar in optics, namely: reflection, refraction and diffraction. 
Two of the early papers! ? treated ultra-short-wave propagation from this 
viewpoint. It was soon observed that diffracted signals tended to be un- 
stable in the shadow region; furthermore, as the wavelength is decreased the 
shadows cast by obstacles such as hills or the bulge of the earth itself become 
more sharply defined. For these reasons, a considerable part of our ex- 
perimental work has been done on paths for which a line-of-sight exists be- 
tween transmitter and receiver. The chief interest, therefore, has been in 
ground reflections and the effect of the atmosphere. 

1J. C. Schelleng, C. R. Burrows and E. B. Ferrell, ‘‘Ultra-Short-Wave Propagation”’, 
Proc. I. R. E., vol. 21, pp 427-463; March 1933. 

2C. R. Englund, A. B. Crawford and W. W. Mumford, ‘Some Results of a Study of 


Ultra-Short Wave Transmission Phenomena’’, Proc. I. R. E., vol. 21, pp 464-492; March 
- 1933. x 
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GROUND REFLECTIONS 


Some of our first experiments with ultra-short-waves showed that regular 
reflections were obtained locally from open, relatively flat fields. The 
reflection coefficients were in good agreement with theory. Later, measure- 
ments of propagation between a transmitter located on a hill top and a 
_ receiver carried in an airplane? showed that for near-grazing angles of in- 
cidence, the irregular and wooded terrain, typical of the New Jersey country- 
side, could give rise to regular reflections at wavelengths as short as four 
meters. The depth of the minima in received signal strength, caused by 
wave interference between the direct and ground reflected components, cor- 
responded to a reflection coefficient of about 0.9. In 1939, unpublished 
results obtained over the 39-mile Beer’s Hill-Lebanon optical path (See map 
of Fig. I-1) indicated that for a wavelength of 30 centimeters the reflection 
coefficient was still large, about 0.8. 

More recently, microwave propagation studies have been made over the 
same type of terrain at wavelengths of 3.25 centimeters and 1.25 centimeters 
and the situation in regard to ground reflections seems to have changed 
somewhat. Experiments were conducted over the 12.6 mile Beer’s Hill-Deal 
path in which the height of the transmitting terminal was varied and which . 
also made use of narrow-beam scanning antennas to separate the direct wave 
from a possible ground reflected component. The results showed the 
apparent reflection coefficient to be of the order of 0.2 at 3.25 centimeters 
and to be even less at 1.25 centimeters. Figure I-2 shows typical curves 
of signal level versus transmitter heights for wavelengths of 3.25 and 1.25 
centimeters. Actually, the shapes of the curves can be accounted for better 
by diffraction, for which the hill about two miles from Deal is considered to 
be a straight edge, than by reflection from an assumed average ground 
plane. The true picture is probably a combination of reflection and 
diffraction effects. 

In an effort to minimize ground reflection, over-water paths were avoided 
in the layout of the New York-Boston microwave repeater circuit and as a 
final check a number of variable antenna-height tests* were made in the 
preliminary survey of all sites. A few curves obtained at a wavelength of 
7 centimeters are reproduced in Fig. I-3. Similar results were observed 
during a survey of sites between Chicago and Milwaukee. 

It is concluded, therefore, that although in the wavelength range down to 
30 centimeters, at least, the effects of ground reflection must be taken into 
account in the choice of sites for an optical path radio circuit, in the lower 
microwave range, below say 10 centimeters, scattering and absorption of the 
reflected wave by rough terrain and vegetation usually results in sub- 
stantially free-space propagation under normal conditions when the line of 


* F, F. Merriam was in charge of this work. 
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sight is well clear of intervening obstructions. In order to have a rule-of- 
thumb as to the amount of path clearance desirable, we have suggested that 
the first Fresnel region should be clear of all obstacles. The first Fresnel 
region for a given transmitter and receiver is bounded by points for which the 
length of the path, transmitter to point to receiver, is greater by one-half 
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Fig. I-2.—Variable antenna-height tests on Beer’s Hill-Deal Path. 


wavelength than the direct path from transmitter to receiver; its cross- 
section by any plane perpendicular to the direct path is the first Fresnel 
zone in the sense used in optics. A wave can be transmitted with practically 
no loss through an opening whose area is of the order of the first Fresnel 
zone. Also, in the case of a smooth reflecting surface between transmitter 
and receiver, the first Fresnel zone clearance provides a maximum in re- 
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ceived field strength since the half wavelength path difference plus the 180- 
degree phase change at reflection causes the direct wave and the reflected 
wave to arrive in phase at the receiver. In Fig. I-4, the first Fresnel region 
is sketched on the profile map of a typical microwave link for wavelengths 
of 3 meters and 3 centimeters. 

It should be emphasized that the above remarks on ground reflections 
apply only for rough terrain and for the case of reflection at a distance from 
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the terminals. Variable height experiments involving short distances over 
open flat fields reveal the presence of almost perfect ground reflections at 
wavelengths as short as 1.25 centimeters. For transmission paths over 
water, strongly reflected components are often observed. Reports? of 
experiments in the Arizona desert indicate a strong ground reflection at 
wavelength of 3 centimeters. In such locations, and most likely in the 
plains regions, the presence of substantial ground reflected components may 
prove to be troublesome. nk 


3 Report No.6. Electrical Engineering Research Laboratory, The University of Texas, 
February 1, 1947. 
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ATMOSPHERIC REFRACTION 


As in the case of ground reflection, the refractive effect of the atmosphere 
has been found to play a somewhat varying role, depending upon the wave- 
length employed. In some of the early work on ultra-short-wave propaga- 
tion,! ? the concept of average atmospheric refraction was found to bring 
about better agreement between observed and calculated results. Due to 
the variation of temperature and water vapor content of the atmosphere 
_ with height above ground, the dielectric constant of the atmosphere nor- 
mally decreases with height. The effect of this negative dielectric constant 
gradient -is to cause the path of a radio wave to be bent slightly downward 
toward the earth, thus effectively increasing the horizon distance. It 
has been suggested that a good approximation for average refraction was to 
assume the radius of curvature of the ray to be four times that of the earth. 
This condition is used at the present time to describe a “standard atmos- 
phere.” 

It was soon found, however, that atmospheric refraction could vary be- 
tween rather wide limits depending chiefly on the gradient of water vapor 
with height.4 Refraction effects were found to be greater in summer than 
in winter since the air contains a higher percentage of water vapor in the 
summertime. A diurnal variation in refraction was also observed on over- 
land transmission paths. During the day, rising convection currents and 
surface winds, caused by surface heating of the earth, usually produce a 
well mixed atmosphere near the earth so that “standard” atmospheric 
conditions prevail. On clear nights, however, particularly if the wind veloc- 
ity is low, radiation cooling of the earth may cause a temperature inversion 
in the lower atmosphere; if, also, the water vapor decreases with height, the 
combined temperature and water vapor effects may add to produce a steep 
negative gradient in the dielectric constant. Stormy weather and over- 
cast skies usually result in standard atmospheric conditions. 

Most of the signal variations observed during a two-year study of propaga- 
tion of two and four-meter waves over the 39-mile over-land optical path 
between Beer’s Hill, N. J. and Lebanon, N. J.° could be explained satis- 
factorily on the basis of wave interference between direct and ground- 
reflected radiations; the relative path lengths, and hence the phases, of these 
two components of the received field varied with the refractivity of the 
atmosphere. The fading on the two wavelengths was usually similar in 
major detail as might be expected from the geometry of the path. On the 

1,2 Loc. cit. . | 

+ Englund, Crawford and Mumford, “Further Studies of Ultra-Short-Wave Transmis- 
sion Phenomena”, B.S. T. J., vol. 14, pp 369-387; July 1935. 


5 Englund, Crawford and Mumford, “Ultra-Short-Wave Transmission over a 39 mile 
‘Optical’ Path”, Proc, I. R. E., vol. 28, pp. 360-369; August 1940. 
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few occasions when the fading could not be accounted for in this simple 
fashion, it was assumed that signal components were arriving from above 
by virtue of reflections from small, relatively abrupt changes in the dielectric 
constant of the atmosphere. The existence of such reflections was demon- 
strated by a frequency-sweep method during propagation studies on the 
70 mile over-water path between Highlands, N. J. and East Moriches, Long 
Island.® 

With microwaves, where, as stated previously, ground reflections are 
usually small or absent, one might surmise that changes in atmospheric 
refraction would have a smaller effect on transmission than at ultra-short- 
wavelengths where strong ground reflections are present, and that fading 
should, therefore, be less. Actually the opposite is observed. Fading is 
found to be more frequent, faster, and deeper as the wavelength is decreased. 
This frequency effect may be explained in a qualitative fashion by a con- 
sideration of the relative sizes of Fresnel zones at ultra-short waves and at 
microwaves. It is known that the dielectric constant of the atmosphere 
usually does not vary with height in a smooth linear manner; on calm nights, 
particularly, very steep gradients in the dielectric constant may exist over 
small vertical ranges measuring only tens of feet. .The effectiveness of these 
steep gradients would be expected to depend on their extent relative'to the 
size of a Fresnel zone. Thus ona path such as that in Fig. I-4, a steep grad- 
ient extending over only a hundred feet would include practically the whole 
first Fresnel zone at 3 centimeters while it would cover only a small part of 
a zone at 3 meters wavelength; the effective gradient, therefore, would be 
considerably less at 3 meters than at 3 centimeters. Analyses based on: 
wave theory show that atmospheric layers, in which the dielectric constant 
has a steep negative gradient, tend to confine or guide the radiation in much 
the same way as a waveguide, and that this “trapping” phenomenon, for a 
given layer thickness, becomes more pronounced as the wavelength is de- 
creased.’ 

The mechanism of microwave propagation is certainly a complicated one, 
and a considerable amount of experimental work in the fields of radio and 
meteorology will be required to unravel it. However, it is very difficult to 
interpret the radio measurements in terms of meteorological data. The 
chief difficulty is that meteorological measurements often do not give an 
accurate picture of the atmosphere, particularly at those times when micro- 
wave fading indicates that rapid changes of some sort are occurring in the 

6 Englund, Crawford and Mumford, “Ultra-Short-Wave Transmission and Atmos- 
pheric Irregularities”, B. S. T. J., vol. i7, pp. 489-519; October 1938. 

7H. G. Booker, in England, was the first to call attention to this pihenonicnon: For 
more recent work see: C. L. Pekeris, ““Wave Theoretical Interpretation of Propagation in 


Low-Level Ocean Ducts, Proc. I. R. E., vol. 35, pp. 453-462; May, 1947. This paper ; 
gives references to other work in this field. 
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atmosphere. The instruments used to measure temperature and humidity 
require a few seconds to reach equilibrium—a length of time comparable at 
times with the period of fading. To measure the variation of dielectric 
constant with height, the measuring instruments are usually carried aloft 
by means of captive balloons. A half hour may be required to measure 
to heights of six or seven hundred feet with the result that the final curve 
represents an unknown combination of variations of dielectric constant with 
height and with time. It is also extremely doubtful that the atmosphere is 
uniform in the horizontal plane—an assumption which is usually made in 
the theoretical treatment of microwave propagation. It seems likely that 
the lower atmosphere is far from being a homogeneous fluid but rather may 
contain small air masses or ‘‘boulders’’ with properties which differ con- 
siderably from those of the surrounding air. Reflections from these boulders 
may be.the cause of radar echoes received from the lower atmosphere.® 
Scintillation fading of microwaves is another evidence of these inhomo- 
geneties in the atmosphere. Scintillation fading, a rapid fluctuation in 
signal level about a more or less steady average value, increases as the wave- 
length becomes less and as the path length is increased. 

In order to evaluate, on a statistical basis, the effect of atmospheric 
changes on a typical microwave circuit, extensive measurements of trans- 
mission were made over a 40-mile.overland path between New York City 
and Neshanic, New Jersey. The tests covered a period of about two years. 
Most of the data were obtained at wavelengths of 10, 6.5, and 3.2 centimeters 
although some data were taken at wavelengths of 42 centimeters and 1.25 
centimeters. The results are described in a recent paper.? In many 
respects, observations were in agreement with those made earlier on the 
39-mile Beer’s Hill-Lebanon path at wavelengths of 4 and 2 meters and on 
the 38-mile non-optical path between Deal, N. J. and Lawrenceville, N. J. 
at a wavelength of 2 meters.’ The same seasonal and diurnal trends in 
fading were found; transmission was generally steady during the midday 
hours and during periods of windy or rainy weather; fading was the same 
on vertical and horizontal polarizations. However, the character of the 
fading was different; the fading at microwaves was much faster and deeper 
than that observed on the ultra-short-wave path. The average daily fading 
range for July on the New Yoik-Neshanic path was 20 db at 6.5 centimeters 
compared with a median daily fading range of 8.5 db for 2.0 meters observed 
in July on the Lebanon-Beer’s Hill path. 

8H. T. Friis, “Radar Reflections from the Lower Atmosphere”, Proc. I. R. E., vol. 
35, pp. 494-495; May 1947 (Correspondence Section). 

9ALL. Durkee, “Results of Microwave Propagation Tests on a 40-mile Overland Path”, 
Proc. I. R. E., vol. 36, No. 2, pp. 197-205, Feb. 1948. 


COR, Burrows, A. Decino and L. E. Hunt, “Stability of Two-Meter Waves”, Proc. 
I.R.E., vol. 26, pp. 516-528; May 1938. 
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Other observations on the New York-Neshanic microwave. path may be 
summarized as follows: While all the wavelengths were affected at times 
of anomalous propagation, the shorter wavelengths faded more severely 
and the character of the fading was different from that observed at the 42 
centimeters wavelength; apparently the 3-10 centimeter range was more 
sensitive to the fine structure of the atmosphere, as pointed out previously. 
During non-fading periods, signal levels were very close to the free-space 
values with the exception of the 1.25 centimeter signal which was usually 
15 db or more below the free space value because of atmospheric absorption 
effects. Some special tests showed that fading was considerably more severe 
when one of the terminals was lowered so that the transmission path was 
grazing slightly below line-of-sight. It was also found that fading was about 
twice as great, in decibels, on the whole path as on either half-section. A 
statistical analysis, on an hourly basis, of all the data on 6.5 centimeters 
showed that only one-half of one percent of the total hours had signal 
minima deeper than 20 db below the free space field. Also during August 
1, the day of the most severe fading, the signal was more than 20 db below 
free space for about one-half of one percent of the time. It was also found 
‘that signals of the order of 10 db above free space were equally probable. 
From a consideration of these statistics, it was decided to engineer the New 
York-Boston repeater circuit with —20 to + 10 db allowance for fading. 
on each link. 


SPECIALIZED EXPERIMENTS 


Much of our more recent work on microwave propagation has been of a 
specialized nature in which apparatus and experiments have been designed 
more for the purpose of studying the mechanism of anomalous propagation 
than for making a statistical analysis of the transmission. Perhaps the most 
informative experiments have been those in which narrow beam scanning 
antennas were used to explore the incident wave fronts. 

The first of these antennas had an aperture of 20 feet and a beam width 
between half-power points of 4 degree at the design wavelength of 3.25 
centimeters. It was built for the purpose of establishing a practical limit 
to the size, and hence the directivity, of microwave repeater antennas from 
the standpoint of variations in the angle of arrival 'of the received wave. 
It had been realized, of course, that variations in the refractivity of the 
atmosphere would cause some deviations in the path of the wave. While 
these deviations should be negligible in comparison with the beam width 
of antennas normally used in the ultra-short-wave region, they might be 
comparable with the beam widths readily obtainable in the microwave 
region. 
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Using this antenria for measurements in the vertical plane and another 
identical antenna for measurements in the horizontal plane, angle-of-arrival 
data were obtained during the summer of 1944 over a twenty-four mile, 
partly over-water, path between Beer’s Hill, N. J. and New York and over 
a thirteen mile over-land path between Beer’s Hill, N. J. and Deal, Ni?J.™ 
In the horizontal plane, deviations in the angle of arrival were rather un- 
common and were not greater than + 0.1 degree from the true bearing of the 
transmitter. In the vertical plane, angles of arrival above the true elevation 
of the transmitter were observed to be as much as 0.5 degree on the New 
York path and 0.3 degree on the Deal path during times of anomalous propa- 
gation. From these measurements it was concluded that microwave re- 
peater antennas could be made highly directive in the horizontal plane but 
should have beam widths somewhat-greater than 3 degree in the vertical 
plane unless means for steering the beams are provided. 

Although the 3 degree beam width of these scanning antennas was sharp 
enough to permit the separation of the direct and the water-reflected com- 
ponents on the New York path, and to demonstrate the anomalous behavior. 
of each, there was evidence that occasionally there were signal components 

‘so close together in angle that a sharper antenna would be required to resolve 
them. Consequently a scanning antenna of the metal-lens type was con- 
structed for operation at 1.25 centimeters. The aperture of this antenna 
was 20 feet in the long dimension; the beam width was 0.12 degrees. Using 
this antenna and also the 3.25 centimeter scanning antennas, angle-of- 
arrival measurements were made in the summer of 1945 on the Deal-Beer’s 
Hill path. The most noteworthy result of these observations was the 
demonstration of multiple-path transmission. Two, three and, at times, 
four distinct signal components were observed simultaneously during one 
night when the transmission was extremely disturbed. These transmission 
paths generally were above the true direction of the transmitter; at one time, 
a weak signal was arriving at an angle of 0.75 degree relative to the line of 
sight. These components varied in angle of arrival and in signal amplitude. 
Wave interference among them caused severe fading on broad beam antennas 
that would accept all the wave paths. 

Another significant result of these angle-of-arrival measurements was 
evidence that the transmission mechanism was very similar for wave- 
lengths of 3. 25 and 1.25 centimeters. Angles of arrival, measured simul- 
taneously at the two wavelengths, agreed very well for times of single-path 
transmission; multiple-path transmission was observed on both wavelengths 
although the 3.25 centimeter antenna was too broad to resolve the com- 

11W. M. Sharpless, ‘“Measurement of the Angle of Arrival a Microwaves”, Proc. 
I. R. E., vol. 34, pp. 837-845; November 1946. 


AB. Crawford and W. M. Sharpless, “Further Gleciva tions of the Angle of Arrival 
of Microwaves”, Proc. I. R. E., vol. 34, pp. 845-848; November, 1946, 
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ponents completely. This result suggested that the 1.25-centimeter scan- 
ning antenna might be a very useful tool for investigating the fading mech- 
anism at 7 centimeters wavelength. 

~ The 22.8-mile path between Crawford Hill and a hill on the Murray Hill 
Laboratory property was chosen for study as a representative link in a 
repeater circuit. (See Profile Map, Fig. I-4.) Transmitters for the 1.25- 
centimeter and 7 centimeter wavelengths were installed in the 100-foot tower 
at Murray Hill. At the Crawford Hill receiving site were the narrow-beam 
scanning antenna and a broad beam antenna for 1.25-centimeter operation; 
also two broad beam antennas, spaced vertically 15 feet, for 7-centimeter 
operation. In addition, a 1.25-centimeter radar could be operated with the 
scanning antenna. A corner reflector target, 55 feet on a side, was located 
at the Murray Hill tower. The signal reflected by this target was about 10 
db stronger than the spurious reflections from other objects at the same 
range. By making use of this target and ground echoes at intermediate 
distances, the radar technique provided a considerable amount of useful 
information concerning the transmission phenomena. | 

Measurements were made on this path during the summer of 1946. As 
had been hoped, the observations showed that transmission on 1.25 centi- 
meters and 7 centimeters was often affected by the same conditions except, 
of course, for atmospheric absorption effects at the 1.25-centimeter wave- 
length. While it was not possible to arrive at explanations for all the fading 
observed, the deep minima in the 7-centimeter signal, ie., fades to levels 
of 15 to 20 db or more below the free space field, usually were the result of 
one of three types of propagation*:— 

Type 1. The 7-centimeter fading was of the rapid, large amplitude type 
characteristic of wave interference. - The 1.25-centimeter scanning records 
showed the presence of multiple-path transmission in which two or more 
readily separable wave paths were observed. While the signals on both of 
the vertically-spaced 7-centimeter antennas faded about the same in ampli- 
tude, their signal minima did not occur simultaneously. A space diversity 
system would be successful in reducing the effects of this type of fading. 

Type 2. The 7-centimeter fading was somewhat slower than in Type 1, 
but still had the appearance of wave interference. The 1.25-centimeter 

_ Scanning records appeared to be of the single path variety. However, close 
inspection showed that, in all probability, more than one transmission path 
was involved but the 0.12 degree beam of the antenna was not sharp enough 
to resolve them. ‘The signals received on the vertically spaced 7-centimeter 
antennas faded together so that space diversity would not be expected to be 
successful unless an extremely large spacing of antennas were used. 

* Recently, on a different overland path having barely one Fresnel zone clearance, 
an important fourth type has been observed when atmospheric refraction gives the ray 


path a curvature opposite to that of the earth, thus effectively reducing the path clear- 
ance. 
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Type 3. The 7-centimeter signal would fade to a low level and remain 
there for a considerable period of time; sometimes for an hour or so. The 
character of the fading was unlike that caused by wave interference. The 
1.25-centimeter signal was simultaneously at a low level and the scanning 
records showed that only one path was involved. Reception was almost 
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identical on the two vertically-spaced 7-centimeter antennas. Radar 
observations suggested that this type of fading was due to attenuation by a 
reflecting layer in the atmosphere at a height intermediate to the heights of 
the transmitters and receivers. It was observed, for example, that while 
the echo from the Murray Hill corner reflector was absent, strong echoes 
were received from the hill directly in front of Murray Hill and some 250 
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feet.lower in height; sometimes multiple paths were observed with this echo. 
Space diversity would fail to improve transmission under these propagation 
conditions, and no other means of improvement is apparent except, perhaps, 
an.alternate path. Fortunately, this type of fading was the least frequent 
of the three types which were characterized by low signal levels. 


RAIN ATTENUATION AND ATMOSPHERIC ABSORPTION 


Attenuation effects due to rainfall and absorption by atmospheric gases 
become increasingly important at the short-wave end of the microwave 
region. Measurements of rain attenuation have been made at the Holmdel 
Laboratory: ™; the results are summarized in Fig. I-5. These curves show 
that for wavelengths above about 5 centimeters, rain attenuation is not very 
serious except for rains of cloudburst proportions. However, at wavelengths 
of one centimeter and less, even moderate rainfall will cause ae attenu- 
ations on paths of the order of 10-20 miles in length. 

Absorption by atmospheric gases, principally water vapor and oxygen, 
becomes important at wavelengths below about 1.5 centimeters. According 
to the theoretical work of Dr. J. H. Van Vleck, Harvard University, water 
vapor has an absorption band at 1.33 centimeters and oxygen has bands at 
0.5 and 0.25 centimeters. Measurements made on the Deal-Holmdel path 
at 1.25 centimeters were in fair agreement with Van Vleck’s results and 
indicated that a typical value of atmospheric absorption for this locality 
In ‘summertime is about 0.4 db per mile." 


SUMMARY 


In the ultra-short-wave region, transmission has been found to be affected 
mainly by ground reflections and variable atmospheric refraction; only 
occasionally are atmospheric reflecting layers and trapping phenomena in- 
volved. These wavelengths ordinarily are not transmitted to great dis- 
tances along the surface of the earth, but are diffracted around obstacles. 
They are used for local broadcasting and mobile radio communication. 

Microwaves are attractive for radio repeater circuits since they permit 
the use of wide transmission bands. Ground reflections are apparently of 
small importance with terrain such as that of the Eastern seaboard and sub- 
stantially free-space propagation is obtained during non-fading periods over 
optical paths which have approximately “first Fresnel region” clearance. 
Atmospheric reflecting layers and trapping phenomena are frequently ob- 
served and signal variations are considerably greater than in the ultra- 

13 Sloan D. Robertson and Archie P. King, ‘The Effect of Rain upon the Propagation 
AOL ioie the 1- and 3-Centimeter Regions”, Proc. I. R. E., vol. 34, pp. 178P-180P; 


14 G, E. Mueller, “Propagation of 6-Millimeter Waves”, Proc. I. R. E., vol. 34, pp. 181P- 
183P; April 1946. 
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short-wave region. Although fading becomes worse as the wavelength is 
decreased, the advantages of increased antenna gain and directivity possible 
at the shorter wavelengths suggest the use of a wavelength just above the 
region where rain attenuation becomes objectionable; i.e. above about 5 
centimeters. ; 

The use of two antennas, operated in space diversity, should reduce the 
effects of fading caused by multiple-path transmission. The use of space 
diversity may be essential in those localities where strong ground reflections 
are present. On the basis of the comparatively weak ground reflections © 
measured on the New York-Boston path it was decided to avoid the com- 
plications that would result from the use of space diversity in this experi- 
mental system. . 


Il. REPEATER CIRCUIT PLANNING 


The diagram in Fig. II-1 shows schematically'a repeater circuit. At the 
input terminal toward the left the signal, S, is fed to the terminal’s trans- 
mitting antenna. The radiated signal is propagated as discussed in Section 
I and produces the signal power s; in the output of the receiving antenna 
of repeater 1. The signal is then amplified Gi times and radiated toward 
repeater 2 and this process is repeated until the signal finally appears in the 
output terminal towards the right. In each repeater the signal is gain- 
controlled automatically for the same level of output powers, ie., S = Si = 
Se =-+:-=S,. It is assumed that the signals are amplitude or frequency- 
modulated C.W. carriers of substantially the same frequency in each link 
and that the repeaters have linear amplifiers. The diagram shows a West- 
to-East circuit only. A circuit for the opposite direction requires duplica- 
tion of all the equipment with the exception of the antenna supporting 
towers. ; 

Some simple formulas for the repeater gain and the signal-to-noise ratio 
at the terminal will be given in this section, without going into any details 
on propagation phenomena, antennas, amplifiers, etc. The formulas will 
orient the reader in regard to the importance of quantities such as: 


d = Repeater separation 
\ = Wavelength of signal same units of length 
A = Effective area of each antenna 
F = Noise figure of each repeater amplifier 
B = Bandwidth of circuit in cycles/sec. 
The free space attenuation (Sz_1/sz) of link “x’? is! 
2, 2 
i= 
A2 


15H. T. Friis, “A Note on a Simple Transmission Formula’, Proc. I. R. E., vol. 34, 
No. 5, pp. 254-256; May 1946. 
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; Allowing the signal power to fade by a factor M below this value, the maxi- 


mum gain of repeater ‘“‘x’? must be 
dzn” eet 
G; = “Az M. é . (IT-1) 


The total maximum gain in the circuit is Gr = Gi X G, X ++ Ga which 
for the same repeater spacings and fading allowances becomes’ ~ 


2 M\" er Did “Be ue 
Gp = G" = (2%) = ( ~ (II-2) 





where D is the total length of the circuit. Because of distortion, original 
costs, and maintenance costs, this total gain should be made as small as 
possible. 
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The following example illustrates the maximum gain required of the ampli- 
fier in a repeater: 

For d = 4 X 10‘ meters (25 miles), \ = 0.075 meters (f = 4000 mega- 
cycles), A = 4.6 meter? (50 sq. ft.) and M@ = 100 (20 db), we have 


G = 4.3 X 10’ (76 db) 


The noise output of a repeater due to noise sources in the repeater itself 
is approximately’® 


Nz=4X 107 FBG, Watts |. 
or from (II-1) | 
dz 
Aa * 


16 H, T. Friis, “Noise Figures of Radio Receivers”, Proc. I. R. E., vol. 32, pp. 419-22; 
July 1944, 





N.=4X107F 
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This noise power is transmitted without gain or loss to the output terminals 
of the repeater circuit. The total noise power at the output terminals is 
therefore 


2 
=N,=4xX107 FB A 52a M;) 


Assuming the same repeater spacings and fading factors M in each link and 
substituting D for nd, 
2 


ZN = 4 x 10 rep > M Watts (II-3) 


The signal-to-noise ratio at the output terminals is S/2N. The circuit 
should be designed for a signal power, S, as low as possible. Therefore, it 
is very important to choose values of the several factors in (II-3) which 
give a low level of output noise. 

Assuming a noise figure F = 20 (13 db) and bandwidth B = 10 mega- 
cycles, eight repeaters of the type described in the above example will have, 


=N = 2.8 X 10 * Watts 


or, assuming a required minimum output signal to noise ratio of 30 db, the 
output power must be S 2 0.28 Watts. 

In this example it has been assumed that the signals in all links have 
faded simultaneously 20 db below the free space value, which may only — 
happen a fraction of a percent of the time. Most of the time the signal-to- 
noise ratio will be higher than the assumed 30 db and under normal trans- 
mission conditions it will be 50 db (fading allowance factor M = 1). 

Assuming the same repeater spacings and fading allowances in each link, 
equation (II-1) and (II-3) give the following formula for the ratio of the 
output power to noise figure of the repeater amplifier . 


S/F =4X%10GB(S/2N)n_ (II-4) 


Equations (II-1) and (II-4) are the important design equations for the re- 
peater amplifier. 

The factors in (II-2) and (II-3) will now be discussed briefly. (II-3) 
shows that the noise figure F should be as small as possible. If by improving 
the equipment the noise figure could be halved, then the signal power S 
could also be halved (unless interference from other microwave circuits 
predominate). Later on the noise figure will be discussed further. 

The bandwidth B is determined by the characteristics of the signal it is 
desired to transmit and by the method of transmission. Our aim has been 
to provide 10-megacycle bands which are sufficient for transmission of stand- 
ard television signals by AM or low index FM. 
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An increase in the effective area, A, of the antennas reduces both the total 
gain required of the amplifier and the output noise power. Crosstalk be- 
tween the several antennas in a repeater station and interference from outside 
sources also decrease as the antennas are increased in size because of the 
increased directivity. Therefore, the antennas should be as large as main- 
tenance and initial costs will permit. Antennas will be discussed in detail 
in Section IIT. 

Equations (II-2) and (II-3) show that the wavelength ) should be small. 
Also more frequency space or signal channels may be had at shorter wave- 
lengths. On the other hand, the fading factor M increases somewhat as the 
wavelength is decreased and, besides, attenuation due to rain sets a lower 
limit for \ in the region of 5centimeters. The status of the apparatus art has 
also been an important factor, but it now permits utilization of the wave- 
length range extending upward from 3 centimeters. Since the war, our work 
has been concentrated on a 10% band around 4000 megacycles or 7} 
centimeter wavelength. The manner in which this 4000-megacycle band 
may be divided up into separate channels is explained in Section IV. 

The effects of varying the repeater separation d will now be discussed. 
d appears in the denominator of the exponent of (II-2) which indicates that 
large separations are favorable, while (II-3) shows that a decrease in separa- 
tion cuts down the noise. Small separations are very costly, the cost being 
almost inversely proportional to the separation. We have concluded from 
propagation studies and site surveys that in the eastern part of the United 
States it is desirable to use separations of about 30 miles, which generally 
provide line-of-sight paths with reasonable tower heights. It should be 
mentioned that the fading allowance factor M is not micependent of d; 
an increased d requires a larger fading factor. 


TII. ANTENNA RESEARCH* 


There are three electrical characteristics which repeater antennas should 
possess. The first is high gain (large effective area), as this will reduce the 
path loss and accordingly the requirements on transmitter power. The 
second is good directional qualities so as to minimize interference from out- 
side sources and also interference between adjacent antennas. The third 
is a good impedance match so that reflections between the antenna and the 
repeater equipment will not distort the transmitted signals. These char- 
acteristics should preferably be attainable without the imposition of severe 
mechanical or constructional requirements. 

It was felt that a 10-foot round or square antenna would be the largest 
that maintenance and initial cost would permit. Propagation studies also 


* This section prepared by W. E. Kock who performed the major part of the work on 
antennas. 
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showed that the variations in angle of arrival of the distant signals would be 
small compared to the beam width of 10-foot antennas. First experiments 
were therefore made with 10-foot diameter parabolic “dish” type antennas. 
The experimental models were made of wood with a metallized reflecting 
surface consisting of silver'conducting paint. Fairly satisfactory tolerances 
were met in these first models, but it was anticipated that trouble would be 
experienced in constructing a permanent metal paraboloid of that size to the 
required tolerances without the use of a heavy and costly supporting means 
for the parabolic sheet. It was also found that an ice coating a quarter 
wavelength. thick on the reflecting surface, when wet, acted as an effective 
absorber of power, since the sheet of water is resistive and is backed up by 
the reflector. Such a condition could produce an intolerable drop in re- 
ceived signal and would have to be prevented by providing the dish with a 
plastic cover. As this cover should preferably house the feed also, it would 
have presented a difficult supporting problem. 

- Two electrical shortcomings of the paraboloid antenna also presented 
themselves. First, it was found extremely difficult to obtain a satisfactory 
impedance match between the antenna and the feed line. This was true 
partly because of energy reflected from the dish re-entering the feed horn, 
(this produced a constant 0.6 db standing wave ratio in the feed line), and 
partly because of the problem of matching the feed horn itself over the de- 
sired 400 megacycle band. Secondly, the mutual interference or “crosstalk” 
between two paraboloids was found to be only 50 to 60 db down when placed 
back to back** (Fig. TIT-1). 

A type of reflector antenna was later investigated, which, Senaceh larger 
physically than.a dish having the same aperture area, overcomes the above 
two objections.’ It is shown sketched in Fig. III-2. The photograph of 
Fig. III-3 shows the antenna lying onitsside. It can be seen that the feed is 
effectively ‘offset’ and reflection back. toward the feed is eliminated; the 
experimental model of Fig. ITI-3 showed only 0.1 db standing wave ratio in 
the feed line over a 10% band of frequencies. Furthermore, a horn or 
“shielded” type feed is used which confines the energy and minimizes stray 
radiation, and measurements indicate that the back-to-back crosstalk sup- 
pression of two such antennas will be high. This long horn is also partly 
responsible for the excellent impedance match. A horn having a large 
aperture ‘‘matches” free space quite well and the slight mismatch at the 
_ throat can be tuned out over a wide band of frequencies. This is not true 

+ A waveguide termination in common use today employs a resistive sheet placed one- 
quarter wavelength in front of a conducting plate; this device absorbs practically all the 
power falling on it. 

** Back to back crosstalk suppression in the order of 125 db would be desirable for 


repeaters receiving and transmitting on the same frequency. 
7U. S. Patent *% 2,236,393, H. T. Friis and A. C. Beck. 
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Fig. I1I-1.—Measuring the back to back crosstalk of two 10 ft. 





paraboloid antennas, 
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of the short, stall aperture horn used in feeding the dish antenna. The 
antenna displayed an effective area which was 66% of its actual area which 
is only 0.9 db below the theoretical maximum. 


PARABOLOIDAL 
SURFACE 





: AXIS OF REVOLUTION [eee FOCAL POINT 
OF PARABOLOID . OF PARABOLA 


Fig. III-2.—Schematic of horn-reflector antenna. 





Fig. ITI-3.—Experimental model of horn reflector antenna. 


The expected gain and directional characteristics of an antenna can be 
realized only if the emerging wave fronts are truly plane. Since deviations 
greater than -t g wavelength can materially impair the antenna perform- 
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ance,!§ reflector antennas have difficult tolerance requirements imposed 
upon them. For example, at 4000 megacycles, the 10-foot reflector must 
conform to parabolic shape to within + § inches, and any twist or warp 





Fig. ITI-4.—Shielded metallic lens antenna. 


of the reflector greater than this would be objectionable. Lenses, however, 
possess the property that a twist or warp in them does not impair their 


18 See, for example, “Radar Antennas”, H. T. Friis and W. D. Lewis, B. ie ea ee 
vol. 26, p. 219, April 1947, Figs. 17 and 28, 
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beam-forming properties, and, with the development of metal lenses for 
microwaves,’® this type of antenna appeared to lend itself very well to re- | 
peater applications. ‘The “shielded” type lens, which is a lens in the mouth 
of a short horn, is shown in Fig. III-4. ‘This antenna, which was developed 
for the New York-Boston circuit,* possesses the property of excellent 





Fig. III-5.—Internal view of lens for the shielded metallic lens antenna. 


crosstalk suppression both side to side (85 db) and back to back (125 db). 
Within the horn, the small amount of energy reflected back from the lens 
is directed away from the feed by tilting the lens, a procedure which does not. 
noticeably affect the radiation characteristics, but which results in a fairly 


19 W. E. Kock, Metal Lens Antenna, Proc. I. R. E., vol. 34, Bs 828, November beasts 
* Developed by W. E. Kock and R. W. Friis. 
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good impedance match (under 0.8 db standing wave ratio) over the desired 
400-megacycle band of frequencies. The lens in the mouth of the horn also 
provies a convenient support for a plastic impregnated Fiberglass sheet 
which acts as a weatherproof cover and protects the lens against ice forming 
between the plates. 

The lens itself, Fig. ITI-5, is based on waveguide principles and causes the 
wave to be refracted by virtue of the fact that waves confined between plates 
parallel to the electric vector acquire a phase velocity higher than their free 
space velocity in accordance with the equation: 


X 2 
Viens = Ufree ica), / = (3) ; (III-1) 


where ) is the wavelength and, a, the distance between the plates. The 
index of refraction is thus less than one, and a converging lens must be made 
concave. 

As seen in Fig. ITI-5, the lens is stepped to reduce its Geielness ‘ha a 
consequence of this stepping, the efficiency at midband of the antenna 
(50%), is a good deal less than the theoretical value of 81%. Furthermore, 
the index of refraction varies with wavelength, as seen from equation III-1, 
and this results in a defocussing of the lens, with a consequent drop in gain, 
at wavelengths different from the design wavelength. This amounts to a 
drop in gain of 1.5 db at the edges of a 400 megacycle band; however, its 
other characteristics of impedance, side lobe suppression (70 db in the two 
rear quadrants), crosstalk, and ease of construction, help to make up for 
the gain deficiency. 

Measurements taken on the antenna when a thick coating of ice had 
formed on the plastic cover indicated that the impedance match was im- 
paired (the maximum standing wave ratio increased from .8 db to 1.6 db), 
but that the gain was not appreciably affected (less than 1 db). Since | 
propagation experiments indicate that severe atmospheric fades are not 
likely to occur during the winter months, some of the fading allowance can. 
be applied against ice loss. | 

_ There was some doubt that the crosstalk figures quoted above could be 
relied upon during heavy rainfalls, as there was indication that the signal 
transmitted from one antenna might be reflected from the rain and thus 
caused to enter an adjacent side-by-side antenna. Measurements during 
a moderately heavy rain proved that this effect was small, but large enough 
so that the 85 db side-to-side nEure was approaching a limit for the 4000 
megacycle band. 

The measurement of antenna characteristics involves microwave tech- 
niques whose development is an important part of a research program. 
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Fig. III-7.—Schematic of metallic delay lens for repeater applications using metal 
strips as the delay elements. 
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The antenna measuring methods which were employed in our repeater 
research follow along the lines of those described in a recent paper.” 
The very large signal ratios of 120 db or more necessitated double detection 
receivers and low noise figure crystal converters. Pattern and gain meas- 
urements of the antennas required measuring sites having large unob- 





Fig. III-8.—A view of the partly assembled lens of Fig. ITI-7. 


structed areas; these measurements were conveniently taken at the Holmdel 
Radio Research Laboratory. Impedance measurements involved the usual 
microwave equipment such as standing wave detectors in waveguide form, 
signal generators and calibrated receivers. 

Research is now underway on an improved metal lens with gain and band- 
width properties which are superior to the lens of Fig. I1I-5. These lenses” 


2C. C. Cutler, A. P. King, W. E. Kock, “Microwave Antenna Measurements”’, 
Proc. I. R. E., vol. 35, No. 12, pp. 1462-1471, December 1947. 

21 Winston E. Kock, “Metallic Delay Lenses”, B.S. T. J., vol. 27, pp. 58-82, January 
1948. 
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employ an artificial dielectric material which duplicates, on a much larger 
scale, processes occurring in a true dielectric. This involves arranging 
conducting elements in a three dimensional array or lattice structure to 
simulate the crystalline lattice of the true dielectric. Such an array responds 
to radio waves just as a molecular lattice responds to light waves, and if the 
spacing and size of the elements is small compared to the wavelength, the 
index of refraction is substantially constant, so that lenses made of this 
material are effective over large wavelength bands. . 

A lens employing conducting spheres as the lattice elements is sketched 
in Fig. IIJ-6. A more convenient structure for large lenses is shown in 
Fig. III-7 and III-8; it uses thin metallic strips, with the width dimension 
parallel to the electric vector. Slotted polystyrene foam sheets support 
the strips and they are stacked up to form the lens. A quarter wavelength 
step in the lens causes the reflections from the lens surfaces to cancel at the 
feed point, which, inthe drawing, is the apex of the horn shield. 

Over a 10% wavelength band, a 6 foot square shielded lens antenna of 
this type exhibited an efficiency of better than 60% and the impedance mis- 
match due to the lens produces only a 0.2 db standing wave ratio in the 
feed line. This antenna thus retains the dimensional tolerance, weight, 
size and crosstalk advantages of the shielded lens over the shielded reflector, 
and has the advantage of higher gain and broader band performance over 
the shielded metal plate lens. | 


IV. FILTER RESEARCH* 


Frequency space for common carrier radio relay systems is available in 
blocks several hundred megacycles wide. Where heavy traffic is to be 
carried such bands must be efficiently exploited. This may in time be 
accomplished by using extremely wide band amplifiers, for example traveling 
wave tubes; however, more immediate success is offered by the possibility 
of operating a number of narrower band circuits of different frequencies. 
This could be done by using a separate transmitting and receiving antenna 
for each circuit. But each antenna, for sound technical reasons, must be 
large and expensive and in addition requires adequate tower support. 
Consequently there is a need for filters which can connect a number of 
individual radio channels to a common antenna. 

The design” of these radio frequency branching filters must be coor- 
dinated with the design of the relay system asa whole. At lower frequencies 
where little or no antenna crosstalk protection can be counted on it is natural 

* This section prepared by W. D. Lewis and L. C. Tillotson who were responsible for 
a large part of the research on filters. 


22 For more detailed discussion see, W. D. Lewis and L. C. Tillotson “A Constant 
Resistance Branching Filter for Microwaves,” B. S. T. J., vol. 27, pp. 83-95, Jan. 1948. 
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to lump transmitting frequencies in one group and receiving frequencies in 
another. When separate microwave shielded lens antennas are employed 
for transmitting and receiving in each direction it becomes practical to use 
a frequency plan in which transmitting and receiving frequencies are in- 
terleaved. Such a plan eases filter requirements considerably and has 
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Fig. IV-1.—Schematic diagram of a possible five channel radio repeater station, 


certain other advantages to be discussed in Section V. A possible re- 
peater employing such a frequency scheme is illustrated schematically in 
Fig. IV-1. 

If a radio frequency branching filter is to fit properly into a repeater it 
must separate or combine channels without excessive loss of signal. In 
addition it must provide an excellent match to the long line which leads 
to the antenna, otherwise troublesome echoes in this line may be caused. 
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Because of IF amplifier band-pass characteristics, suppression requirements 
on the filter, except possibly in the vicinity of receiver image bands, are not 
large. . 

Microwave band-pass filters consisting of one or more cavities arranged 
in sequence along a waveguide have been known for some time. The fre- 
quency, bandwidth and discrimination characteristics of such filters can 
all be chosen within wide limits by appropriate design of the cavities and 
the means for coupling to them. ‘These filters are analogous to lumped- 
circuit channel-passing filters and can in principle, like them, be connected 
in groups to provide a branching network. 
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Fig. IV-2.—Photograph of a branching filter for an experimental radio relay system. 


Several successful two-branch networks have been designed and con- 
structed in this manner. One of these, developed for the New York- 
Boston circuit*, is illustrated in Fig. IV-2. Here two three-cavity filters 
are connected toan E plane Y junction, the waveguide analogue of a series 
connection. The two filters are.tuned to different bands and each is con- 
nected to the junction through a line of length such that it causes no disturb- 
ance in the channel of the other. The electrical characteristics of this 
filter are plotted in Fig. IV-3. 

Problems connected with the design of suitable microwave branching 
filters with more than two branches evidently differ considerably from pre- 
vious filter problems. Channel passing networks which can be connected 
in series or parallel to form a channel branching filter can be designed at_ 
lower frequencies on the basis of lumped circuit theory and built of coils, 
condensers and resistances, but in the microwave region simple elements 

* Developed by the group concerned with high-frequency filter design headed by 
A. R. D’ heedene. : A large part of the research underlying the design of these filters was 


performed by W. W.. Mumford. Prior to the war a considerable amount of research on 
the band-pass type of waveguide filter had already been done by A. G. Fox. 
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and connections do not exist. Where more than two waveguide channel 
passing filters are to be connected to a common junction the design becomes 
complex, since in every channel the sum of the interactions of all the inac- 
tive filters on transmission through the active filter must be zero. It is 
evidently not easy to satisfy this condition, particularly since in doing so 
one must take account of the change with frequency of the effective length 
of all waveguide connecting lines. And even if such a solution is found it 
will be valid for only one set of channels, so that the problem must be solved 
all over again for every change in channel arrangement. 
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Fig. IV-3.—Input standing wave ratios of filter of Fig. IV-2. 


As a result of these difficulties and after a few attempts to overcome them 
it became evident that a more flexible microwave branching filter technique 
should be found. Accordingly a solution on an iterative basis was 
developed. A channel dropping circuit was devised which, when inserted 
in a line, could extract or insert one channel while allowing others to pass 
through without disturbance. This circuit is of the constant resistance 
type; in other words it operates by diverting energy selectively but not 
by reflecting it back to the input. Consequently N such circuits placed in 
sequence do not interact reflectively; they, thus, form an N channel branch- 
ing filter which is also of the constant resistance type. 

An individual constant resistance channel dropping circuit is illustrated 
schematically in Fig. IV-4. It is made up of two hybrid** circuits, two 


23 For a general discussion of hybrid circuits see W. A. Tyrrell, “Hybrid Circuits for 
Microwaves”, Proc. I. R. E., vol. 35, No. 11, pp. 1294-1306, November 1947. 
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identical band reflection filters, and two quarter wavelengths of line. Each 
of the hybrids is analogous to a low-frequency hybrid coil and operates as 
follows. A wave in line C (See Fig. IV-4) incident on the hybrid is divided 
equally and with equal phase into A and B but does not appear in D or 
reappear in C. If waves in A and B are incident on the hybrid a wave pro- 
portional to their vector sum will appear in C, a wave proportional to their 
vector difference will appear in D but nothing will reappear in A or B. 
A wave in the input line incident on the channel dropping circuit will thus 
be divided by the input line into the lines leading to the two band reflection 
filters. These filters are designed to reflect frequencies lying within their 
band and pass all other frequencies. If the frequency is outside of the re- 
flected band the two waves will travel on to connections A and B of the 
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Fig. [V-4.—Schematic diagram of a constant impedance channel dropping filter 
using hybrid junctions and band reflection filters. 


output hybrid. Here they will have equal phase and amplitude, their 
vector difference will be zero and the wave appearing in C of the output hy- 
brid and consequently in the output line will contain all the power. If 
the frequency lies within the band of the reflection filters the two waves will 
be reflected by them and will travel back to the connections A and B of 
the input hybrid. The two waves strike these connections with opposite 
phase since one of them has traveled twice over an extra quarter wavelength 
of line. Their vector sum will consequently be zero and the wave which 
appears in terminal D of the input hybrid and consequently in the dropped 
channel line will contain all the power. The circuit of Fig. IV-4 is therefore 
_ a constant resistance channel dropping network which diverts energy lying 
within the band of the reflection filters but allows all other energy to pass 
through without disturbance. Conversely, by the law of reciprocity, this 


s 
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circuit can insert energy lying within the band of the reflection filters into 
the main line without disturbing any energy passing through it at other - 
frequencies. 
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Fig. [V-.6—Waveguide band reflection filter used in the filter of Fig. IV-4. 


An embodiment of the circuit of Fig. IV-4 suitable for use in the repeater 
arrangement of Fig. IV-1 has been constructed and tested. Figure IV-5 
illustrates the waveguide hybrid employed. Here the waveguide opening 
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Fig. IV-8.—Input standing wave ratios for the filter of Fig. IV-7. 


218 BELL SYSTEM TECHNICAL JOURNAL 


for C is physically parallel to those for A and B and is connected to them 
through a broad-band waveguide taper. Connection D is made through a 
relatively narrow band coaxial and probe arrangement. Figure IV-6 il- 
lustrates one of the waveguide band reflection filters. In this filter reflec- 
tion occurs at three resonant rods, each tuned by an adjustable capacita- 
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Fig. IV-9.—Transmission loss for the filter of Fig. IV-7. 
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tive plug at one end. These rods are placed perpendicular to the electric 
vector of the guide and are coupled to it by means of adjustable screws placed 
over them in the wide wall of the guide. 

Channel dropping units for five channels in the 4,000 megacycle region 
were made up of these components and suitable quarter wavelengths of 
guide. These units were connected in sequence as shown in Fig. IV-7 and 
adjusted systematically. The electrical performance of the resulting five- 
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channel branching filter was measured with a double detection measuring 
set and is plotted in Figs. IV-8 and IV-9. 

These measured electrical characteristics serve as a check on the general 
theoretical ideas concerning constant resistance hybrid branching networks. 
They also indicate that when these ideas are embodied in the form illustrated 
in Fig. IV-7 the result is a branching filter which can be used in currently 
planned radio relay circuits. 

The circuit of Fig: IV-7 provides a satisfactory channel splitting network. 
It does not, however, provide consistent high off-frequency discrimination 
between one of the channel output terminals and the other terminals of the 
filter. When systems requirements* are such that extra discrimination or 
special impedance behavior is required, this can be supplied by inserting 
suitably designed reflecting filters in the branch lines. These added filters 
will not interact reflectively with the branching filter. 


V. THE REPEATER AMPLIFIERT 


In a microwave repeater circuit, Fig. II-1, the signal is amplified at each 
repeater to compensate for the transmission loss in the preceding radio 
path: Since we cannot build perfect amplifiers, the signal will not appear 
at the output terminals as a true replica of the input signal; the circuit will 
distort the shape of the signal and it will also add noise. Therefore, the 
main objectives in amplifier work have been to keep the signal distortion and 
the added noise within certain requirements. 

To be more specific, the repeater amplifier in a relay system must be 
capable of supplying a maximum gain, G, as given by the equation II-1; 
it must have a ratio of output power capacity to noise figure which will meet 
the signal to noise ratio requirements of the system as given by equation 
IT-4; since distortionless transmission is desired, it must have an amplitude 
characteristic as flat as possible and a phase characteristic as linear as pos- 
sible over the essential range of frequencies of the signal it is desired to 
transmit;* and it must be equipped with an automatic gain regulating cir- 
- cuit to hold the output power constant over the expected range of input 
levels. 

The simplest relay amplifier would be one which amplifies the signal and 
sends it on without a change in frequency. However, two major considera- 
tions indicated that early repeater amplifiers could not be so simple. No 

* E.g., the converter may require a reflection in the input line at the image frequency, 
See Section V and Fig. V-4. 

+ Those parts of this section dealing with the general layout, the requirements, and the 
over-all testing of the repeater amplifier were prepared by D. H. Ring, who together 
with A. C. Beck did the work on this phase of the problem. 


*4 Sallie Pero Mead, “‘Phase Distortion and Phase Distortion Correction”, B. S. T. J., 
vol. VII, No. 2, pp. 195-224, April 1928. 
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microwave amplifiers were known which gave promise of yielding an‘ ade- 
quate ratio of output power capacity to noise figure, and there was con- 
siderable doubt of our ability to reduce, sufficiently, the feedback from the 
transmitting antenna to the receiving antenna. 

These difficulties with straight-through amplification can be avoided by 
a repeater amplifier such as is shown schematically in Fig. V-1. The in- 
coming signal is converted to an intermediate frequency, IF, where better 
amplifiers are available and where the major part of the required gain is 
supplied. The amplified IF signal is then converted back to the microwave 
frequency f + Af, where Afis relatively small. The difference Af between the 
incoming and outgoing frequencies permits the use of circuit selectivity to 
counteract feedback troubles, and the radio frequency amplifier following the’ 
transmitting converter can have a relatively large noise figure. 
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Fig. V-1.—Schematic of a repeater amplifier. 


Our initial research on microwave repeaters was directed toward solving 
the problems associated with an amplifier of the type shown in Fig. V-1. 
The gain and level figures shown in this figure apply to the example of an 
eight-link relay system given in section II. They indicate approximate 
minimum objectives for the various components of the repeater amplifier 
to be discussed later. 

Choice of I.F. Frequency—When selecting the intermediate frequency for 
a multichannel repeater circuit utilizing the interleaved radio frequency 
plan of Fig. IV-1 and the intermediate frequency type repeater amplifiers 
of Fig. V-1, the relative position of the various discrete frequencies and 
frequency bands shown in Fig. V-2 must be considered. In order to mini- 
mize the possibility of crosstalk from the image bands and interference from 
the beating oscillators caused by insufficient shielding, it is desirable to 
choose the intermediate frequency in such a way that the image bands fall 
midway between the active bands, and the oscillator frequencies fall midway 
between the image and active bands. ‘These conditions are realized, as 
shown in Fig. V-2, if the intermediate frequency satisfies the relation 
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2IF = nag +P 
or 
If = = (2n + 1) 
where Af is the frequency spacing and is any integer greater than zero. 


In general, better noise figures and circuit stability are obtained with 
low intermediate frequencies, while high intermediate frequencies lead to 
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Fig. V-2.—Frequency spectrum of a multichannel microwave repeater circuit. 


more symmetrical amplitude and phase characteristics. The research on 
the intermediate frequency components described below was conducted in 
the 60 to 70-megacycle range. | 

Frequency Stability—If all receiving and transmitting beating oscillators 
are independently controlled in a multichannel relay circuit of this kind, 
there is a possibility that small variations will add to produce large variations 
at the distant end of the system. This difficulty can be overcome by the 
frequency control system shown schematically in Fig. V-3. The transmit- 
ting and receiving beating frequencies are both derived from an oscillator 
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operating at a frequency suitable for receiving the incoming signal. The 
frequency of this oscillator is controlled by a high Q cavity and a servo 
mechanism to 0.2 megacycles or better.2> One portion of the output of the 
oscillator is used as the beat frequency in the receiving converter. A second 
portion is combined with the output from a crystal oscillator operating at a 
frequency equal to the difference, Af, between the incoming and outgoing 
frequencies. In this way a beat frequency for the transmitting converter is 
obtained which has the same variations as that for the receiving converter 
except for negligibly small variations that may occur in the crystal oscillator. 
As a result of this method of deriving the beat frequencies the outgoing 
frequency always differs from the incoming frequency by an amount equal 
to the crystal oscillator frequency and is not influenced by variations in 
the high-frequency local oscillator. The result is that, except for the small 
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Fig. V-3.—Frequency control system for a microwave repeater. 


crystal oscillator variations, all radiated frequencies in a long circuit, carry 
only the variations in the transmitting oscillator of the originating terminal, 
while the intermediate frequency of each repeater may vary by an amount 
equal to the sum of the variation of its own local oscillator and that of jthe 
terminal transmitter frequency. ; 

Automatic Gain Regulation—The function of the automatic gain control 
circuits is to hold the repeater output constant over the expected fading 
range. As already stated an allowance for fades 20 db down and 10 db up ' 
from free space have been made for 30-mile paths at a wavelength of about 
7 centimeters. In addition to knowledge of the fading range to be com- 
pensated, it is necessary in the design of suitable circuits to know what the 
maximum fading rate is likely to be. Analysis of the records of a number 
of disturbed periods on the New York-Neshanic path indicate a maximum 
rate of 5 db per second. . 


%V.C. Rideout, “Automatic Frequency Control of Microwave Oscillators”, Proc. 
I. R. E., vol. 35, pp. 767-771, August 1947. 
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A conventional delayed automatic gain control circuit has been used in 
which a d.c. voltage, supplied by a peak rectifier in the output of the last 
IF stage, is amplified and fed back to bias several of the IF stages. Analy- 
sis of the transient response of the last repeater in response to a step function 
disturbance in the input to the first repeater, where the number of repeaters 
is of the order of 5 or more, shows that great care must be taken in shaping 
the frequency characteristic of the feedback circuit. 

General Requirements of Components—The important factors bearing on 
the basic layout of the amplifier components have been discussed. Others, 
affecting the design of these various components will now be considered, 
after which a brief description of the research work on each component will 
be given. 

Different repeater circuits will, in general, have different numbers of 
repeaters; also it may be necessary to feed signals into and extract signals 
from them at any point. Under these conditions it is impractical to specify 
only the over-all characteristics for a given number of repeaters. Each 
repeater itself should be individually good and should not depend upon any 
systematic compensation or equalization at any other point in the system. 
The repeater was designed in accordance with this philosophy and, in the 
interest of flexibility and ease of testing, the same line of reasoning was 
extended to cover the design of the components of the repeater as well. 

In accordance with the above considerations major emphasis was placed 
on obtaining a minimum of amplitude variation and phase distortion ‘over 
a 10-megacycle band for each repeater component. However, it was appre- 
ciated that even with the simplest circuits the inherent phase distortion 
would be excessive in long relay systems. Phase equalizers can be designed 
to equalize this distortion, but the difficulties of design and alignment in- 
crease with the magnitude of the distortion to be equalized. Accordingly, 
simple circuits were used wherever gain requirements would permit and at 
the same time parallel research was carried out on the problems of designing 
and testing appropriate phase equalizers. 

While our aim was to provide a repeater 10 megacycles wide suitable for 
any type of modulation, it soon became apparent that it would be un- 
economical to attempt to provide the extreme degree of linearity that would 
be required; for example, for a long relay circuit carrying an amplitude- 
modulated television signal. However, early tests indicated that very satis- 
factory transmission could be obtained with low-index frequency-modulated 
signals, for which reason the later stages of the work were aimed at providing 
an amplifier to be used for such signals. Nevertheless an attempt was 
made to limit the compression in each unit except the R.F. amplifier to 
0.1 db at maximum rated load. 
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A further important requirement placed on all components was that their 
input and output impedances should match the corresponding impedances 
of the components to which they were to be connected with a minimum of 
reflected power over the 10-megacycle band. This requirement was imposed 
on the separate units to provide for flexibility in testing and to permit easy 
patching in of spare units in case of failure. 


RECEIVING CONVERTER* 


The receiving converter, together with the input to the first stage of the 
intermediate-frequency amplifier, occupies a unique position in the repeater 
amplifier in that it is located at that point in the circuit where the signal 
level is lowest. As a consequence, research on receiving converters has been 
directed toward insuring that the least possible noise be added to the signal 
_to be amplified. An extensive background of microwave converter design 
information was available from the work done on converters during the 
war?®.27 which led to the selection of a balanced converter using a waveguide 
hybrid junction and 1N23-B silicon point contact rectifiers. The informa- 
tion already available would probably have been adequate except for the 
two additional requirements imposed by the repeater amplifier: first, that 
the standing wave ratio at the input have a low value and, second, that 
uniform conversion efficiency be maintained over a band of at least 10 mega- 
cycles. 

A receiving converter with its associated input and output circuits is 
shown in Fig. V-4. Frequency conversion is accomplished in a device of 
this kind by virtue of the fact that when two sinusoidal voltages (inthis : 
case the signal and the beating oscillator) are applied to a non-linear imped- 
ance such as a rectifier, new frequencies given by the sums and differences 
of the applied frequencies and their harmonics are generated. The dif- 
ference frequency may thus be selected as the desired output frequency and 
passed through the IF transformers to the IF amplifier. The performance 
of a converter is, however, influenced by the impedances encountered by 
some of the other frequencies generated. For this reason, the separation 
S between the input filter, more properly a component of the channel se- 
lecting network, but here shown as part of the converter, and the converter 
must be given consideration, since it determines the phase of the reflection 
from the filter at the image frequency (the difference between the beating 

-* This section prepared by C, F, Edwards who was responsible for the research on the 
receiving converters. 

25 ‘‘Nevelopments of Silicon Crystal Rectifiers for Microwave Radar Receivers”, J. H. 
Scaff and R. S. Ohl, B. S. T. J., vol. 26, pp. 1-30, January 1947. 

27 Descriptions of several converters developed prior to and during the war as well as- 


a more complete description of the present converter are given in C. F, Edwards’ paper, 
“Microwave Converters”, Proc. I. R. E., vol. 35, No. 11, pp. 1181-1191, November 1947. 
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oscillator second harmonic and the signal) which in turn affects the converter 
‘output impedance and hence the match between the rectifiers ‘and the IF 
transformers. Failure to obtain a proper match results in non-uniform 
transmission over the channel band. 

In addition, to maintain uniform conversion efficiency over a wide band, 
it is necessary to control the impedance encountered by frequencies near the 
beating oscillator second harmonic. ‘This is done by means of the harmonic 
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Fig. V-4.—Receiving converter with input filter. 


filters shown which reflect these frequencies from a point close to the recti- 
fier. Since the point of reflection is close, there is little opportunity for the 
phase of the reflection, and hence the conversion loss, to vary over the band 
of frequencies to be converted. 

The converter in Fig. V-4 was designed to provide as good a termination 
as possible to the incoming waveguide over the band of channel frequencies 
so that a minimum of additional adjustment would be required to match 
the guide accurately at any particular channel frequency. This additional 
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adjustment is provided by the input tuner shown. With it, the reflection 
coefficient may be reduced to zero at the desired channel band center, and 
when this is done the standing wave ratio at the input is less than 2 db over 
a 20-megacycle band. 

The bandwidth of the converter is largely determined by the IF trans- 
formers shown which transform the balanced output impedance to the 
unbalanced 75-ohm coaxial line connecting to the IF amplifier. When 
transformers having a coupling coefficient of about 0.5 are used, transmis- 
sion variations of less than 0.1 db over a 20-megacycle band are obtained. 


SIGNAL 
“INPUT 





Fig. V-5.—Receiving converter with low noise IF preamplifier. 


The noise figure of the repeater amplifier is determined by the conversion 
loss and noise figure of the converter and the noise figure of the IF ampli- 
fier.1® The converter designed for the New York-Boston circuit has a 
conversion loss of about 52 db and its noise figure is 10 db. Thus when it 
is used with an IF amplifier having a noise figure of 7 db, a figure of 14 db 
is obtained for the over-all noise figure of the repeater amplifier. Figure 
V-5 shows this converter with a low noise preamplifier attached.* 


I.F, AMPLIFIERt 


Most of the gain of the IF amplifier is obtained with stages using double 
tuned, symmetrically loaded (or ‘matched’’) interstage transformers, de- 


16 Loc. cit. 

* Developed by H. C. Foreman and B. C. Bellows, Jr. 

{ This section prepared by Karl G. Jansky who, together with V. C. Rideout, did the 
work on IF amplifiers. 
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signed in accordance with the formulas given by a former member of this 
Jaboratory.”® In Fig. V-6, ‘‘a” shows a schematic diagram of such a trans- 
former and ‘‘b” shows the equivalent 7 network generally used. The equiv- 
alent T network shown at “‘c” has also been used. In Fig. V-7 ‘‘a” shows the 
theoretical band-pass. characteristic for this type of transformer with a 
coupling coefficient of 0.5 which was the value used in most stages. The 
circles indicate points measured on a typical transformer. This matched 
network is relatively insensitive to small changes in capacitance so that 
wherever it is used it is possible to change tubes without having to realign 
- the amplifier. The delay distortion for this type of network is, as shown by 
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Fig. V-6.—The double tuned IF transformer. 


“b” in this figure, also relatively small. The gain per stage for a coupling 
coefficient of 0.5 is approximately 6 db for 6AK5 vacuum tubes and about 
12 db for the recently developed WE 404-A tubes. As shown by W. J. 
Albersheim”® it is possible to design circuits which will give more gain than 
the “matched”’transformer, but only at the cost of increased sensitivity to 
capacitance changes. For convenience the IF amplifier is usually divided 
into a preamplifier closely associated with the réceiving converter and a 
main IF amplifier. 

Pre-amplifier—At the time work was begun, noise figures of the order of 


8 V. C. Rideout, “Design of Parallel Tuned Transformers”, B. S. T. J., vol. 27, pp. 
_ 96-108, January 1948. 
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12 to 14 db were obtainable for 65 mc amplifiers of the required bandwidth 
with.6AKS5 tubes and matched input transformers. This was much worse 
than was desired. By using a 6J4 close spaced grounded grid triode in the 
input stage, much lower noise figures were obtained, but the gain with: 
matched input and output circuits was so low (approximately 3 db for a 
coupling coefficient of 0.5) that the noise from the following stage contributed 
considerably to the overall noise figure. 

By removing the loading resistance on the output side of the first inter- 
stage transformer and reducing the coupling coefficient to 0.3, the gain was 
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Fig. V-7.—Band pass and delay characteristics of the IF transformer. 


raised sufficiently to give an overall noise figure of the order of 7 db. Figure 
V-5 shows the preamplifier developed for the New York-Boston circuit. 
This amplifier employs a 6J4 and a WE 404-A and provides a gain of 23 db. 

The 6J4 tube, when used in a grounded grid circuit, has an input imped- 
ance of approximately 85 ohms which is close to the desired 75 ohm imped- 
ance. When a good match is required, it is necessary to use an input 
transformer, but it should be noted that an improvement in the noise figure 
can be obtained by deliberately producing a mismatch in the right direction 
at theinput. Noise figures as low as 4 db have been obtained in this manner 
with recent experimental tubes. 

Figure V-8 is a schematic diagram of an amplifier with a very low noise 
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figure consisting of two of these experimental grounded grid tubes in tandem. 
' The circuits were designed so that there is a mismatch at the input of each 
‘tube. The overall noise figure is 4.0 db with the amplifier connected be- 
tween a 75-ohm generator and a 75-ohm load; the gain is 172 db; and the 
bandwidth at the 0.1 db down points is about 13.5 megacycles. 

Main IF Amplifier—In order to obtain the required output power over 
the desired bandwidth it was necessary, at the beginning of this work, to 
use the WE 367-A tube in the output stage. Since the gain of a stage using 
this tube is very low, it must be driven by a tube similar to the 6AG7 and 
to prevent compression in this latter tube a special high gain, triple-tuned 
interstage transformer was designed, with a relatively low coupling coef- 
ficient. ; 

Methods of paralleling tubes to get more power output or the same power 
output with a much wider bandwidth have been worked out, but recent 





Fig. V-8.—Schematic of a “low noise figure” preamplifier. 


developments in power output tubes and in transmitting converters have 
made them unnecessary, except for applications requiring very wide band- 
widths. 

Automatic gain control can be applied to the IF amplifiers by controlling 
the grid bias of the various stages. However, it is not advisable to apply 
the gain control bias to either of those stages which are preceded by the 
special high-gain transformers. In these stages the slight changes in input 
impedance of the tubes caused by variations in the grid bias would be suf- 
ficient to alter significantly the band-pass characteristics of the transformer. 

Figure V-9 shows an amplifier with about 55 db gain developed for the 
New York-Boston circuit* Automatic gain control bias is applied to the 
grids of the first three stages which employ wide band ie 0.7) matched 
interstage transformers. 

It will be noted that the first interstage transformer of the preamplifier 
and the last interstage transformer of the main IF amplifier require low 
coupling coefficients to obtain the gain desired at these points. For this 


* Developed by A. L. Hopper and B. C. Bellows, Jr. 


230 


BELL SYSTEM TECHNICAL JOURNAL 





Fig. V-9.—The main If amplifier. 
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reason, these two transformers largely determine the band pass and delay 
distortion of the whole amplifier. Typical overall characteristics for a 
complete IF amplifier are shown in Fig. V-10. 
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Fig.-V-10.—Band pass and delay characteristics of complete IF amplifier. 


TRANSMITTING CONVERTER, OR MOoOpDULATOR* 


The transmitting converter problem differed from the first converter 
problem in that high output power was the major consideration, rather than 
low conversion loss and noise figure. . 

There will normally be three frequencies present in the output of a con- 
verter. These are the desired output sideband frequency, /’, the beating 
oscillator frequency, f’ + IF or f’ — IF, and the unwanted sideband fre- 
quency, f’ + 2IF or f’ — 2IF. The strongest of these is the beating os- 
cillator frequency, but fortunately this can be suppressed by 20 to 30 db 
by balance in a balanced converter. Investigation showed that the input 
impedance at the IF terminals of the modulator was affected by the Joad 
impedance of the modulator at both the wanted and the unwanted sideband 
frequencies. The load impedance consists of the input impedance of the 
RF amplifier as seen through the length of transmission line which connects 
the two components. At the wanted sideband frequency, the RF amplifier 
presents a good termination and the load impedance is independent of the 
length of the connecting line. At the unwanted sideband frequency, how- 
ever, the RF amplifier presents a short circuit and hence the load impedance 
is a function of the length of the connecting line. A waveguide filter was 
incorporated in the output circuit of the modulator so as'to reflect the un- 


* This section prepared by W. W. Mumford who did the work on the transmitting 
converters. 
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wanted sideband back into the modulator in proper phase, thereby making 
the length of the line going to the RF amplifier much less critical. The 
proper phase of reflection was obtained by adjusting the length of the 
waveguide between the modulator and the reflecting filter. The IF imped- 
ance of the crystals could thus be varied so as to obtain an impedance 
match between the IF amplifier and the crystals. 





Fig, V-11.—Experimental model of transmitting converter. 


An investigation of the power capacity of various types of standard and 
special silicon crystals indicated that one similar to the 1N28 crystal was 
the best available,for this service, and life tests were made to check the sta- 
bility of the crystal under high level conditions. Special IF input and 
crystal-to-waveguide matching circuits were developed to meet the stringent 
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match requirements. The research on these major factors and many 
subsidiary details resulted in a converter which had an output of 6 milli- 
watts with less than 0.1 db compression, good input and output match, a 
flat amplitude response over more than 10°*megacycles, and a conversion 
loss of about 11 db. 





Fig. V-12.—Unbalanced to balanced coaxial transformer for feeding the 65 mc signal 
to the transmitting converter. 


This converter, or transmitting modulator, had the 1N28 crystals mounted 
in the conjugate branches of a waveguide hybrid junction, as shown in Fig. 
V-11. Adjustable coaxial sleeves surrounded the crystal cartridges in or- 
der to effect an impedance match to the waveguide, and tuning studs were 
provided for trimming adjustments. The beating oscillator was injected 
into the hybrid junction through a broad-band coaxial-to-waveguide trans- 
ducer in the upper branch of the hybrid junction, and the sidebands appeared 
in the conjugate waveguide branch. The 65-megacycle signal was fed onto 
the crystals in push-pull through the unbalanced to balanced coaxial trans- 
former shown in Fig. V-12._ Blocking condensers enabled the crystal cur- 
rents to be monitored separately and RF filters kept the 4000-megacycle . 
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energy from entering the 65-megacycle transformer compartment. The 
converter developed for the New York-Boston circuit, complete with 
sideband filter, spacer, mounting brackets and d-c. leads, is shown in Fig. 
V-13.* 


R. F. AMPLIFIER** 


Prior to the war, a considerable amount of research had been applied in 
the Bell Telephone Laboratories to electron tubes operating on the velocity 
modulation principle, with the expectation that such tubes would find ap- 
plications in radio-relay systems. Although this work was interrupted 
by the war, enough information had been obtained to make it apparent, 
upon resumption of the radio relay work, that such tubes were the only 
ones then known which showed promise of meeting the stated requirements. 

Velocity modulation tubes have been described by several authors,” 3 3! 
and the theory of their operation has been discussed adequately in several 
places.**, 33, 34,35 However, a review in 1939 of the structures then known 
had led to the decision that a new type of construction would be necessary 
to obtain a satisfactory amplifier tube for radio relay purposes. To keep the 
tube voltages within reasonable limits it is desirable to make the input and 
output gaps as small as possible.” Resonant circuits external to the evacu- 
ated envelope are desirable to enable coverage of as large a frequency range . 
as possible with a single tube, and to facilitate addition of broad-banding 
circuits. Grids on the input and output gaps are undesirable because of 
the large interception of current and the difficulty and expense of construc- 
tion. 

With the above considerations in mind, decision was then made to explore 
the possibilities of focussed beams, and of gaps comprised of copper discs 
sealed through a cylindrical glass vacuum envelope. ‘The latter technique 
had been developed at the Bell Telephone Laboratories in connection with 

* Developed by H. C. Foreman and W. W. Halbrook. 

** This section prepared by A. G. Fox and A. E. Bowen. .Messrs. Fox and Bowen, in 
collaboration with A. L. Samuel, A. E. Anderson, and J. W. Clark of these Laboratories, 
did the major part of the research which resulted in this amplifier. 

29 Varian, R. H. and Varian, S. F., “A High Frequency Oscillator and Amplifier”, 
Jour. A pplied Physics, vol. 10, No. 5, PP: 321-327, May 1939. 

30 Hahn, W. C. and Metcalf, G. F., “Velocity Modulated Tubes”, Proc. I. R. E., vol. 
27, No. 2, pp. 106-116, Feb. 1939. 
( 31 oo A. E., “Klystron Tubes”, McGraw-Hill Book Co., New York, 1947. 
Book 

32 Hansen, W. W. and Richtmyer, R. D., “On Resonators Suitable for Klystron Os- 
cillators”, Jour. Applied Physics, vol. 10, No. 3, pp. 189-199, March 1939. 

33 Hahn, W. C., “Small Signal Theory of Velocity Modulated Electron Beams’’, G. E: 
Review, vol. 42, No. 6, pp. 258-270, June 1939. 

34 Hahn, W. C., “Wave Energy and Transconductance of Velocity-Modulated Elec- 
tron Beams”, G. E. Review, vol. 42, No. 11, pp. 497-502, Nov. 1939. 


% Ramo, S., “The Electronic-Wave Theory of Velocity Modulation Tubes’, Proc. 
I. R. E., vol. 27, No. 12, pp. 757- 763, Dec. 1939. 
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the development of water cooled tubes, and showed promise of providing a 
very convenient means for construction of resonant circuits in which a 
part of the circuit was external to the vacuum envelope and part was internal. 

Numerous forms of focussed beam-disc seal velocity modulation tubes 
were constructed for examination of the various factors affecting the per- 
formance of such tubes. Sizes and shapes of gaps, length of drift space, 
voltage and beam current, and other factors were adjusted to optimize the 
performance. Both magnetic and electrostatic focussing of the beam were 
explored. Triple gap tubes, from which gains of more than 30 db were ob- 
tained, were experimented with. 

The end result of this work was the development of a fourstage amplifier 
employing velocity modulation tubes of the disc-seal type (Fig. V-14) 
designed especially for this application.* The electron gun is at the left, 
and the collector is at the right. These tubes employ external circuits in 
the form of resonant cavities assembled around the tube. The electron 





Fig. V-14.—Disc seal velocity modulation tube. 


beams are accelerated by 1500 volts and focussed magnetically by means of 
small alnico permanent magnets placed around the cavity structure as shown 
in Fig. V-15. This shows a single stage of coaxial-coupled amplifier with 
half of each cavity removed to show internal details. The cavities are 
tunable over approximately a 250-megacycle range by means of meta] screw 
plugs located around their periphery. Two different sets of cavities are 
sufficient to cover the entire frequency range of the radio repeater. 

A single amplifier stage of the type shown in Fig. V-15 will exhibit a 
single-tuned type of gain characteristic which is only 5 megacycles wide at 
the 3 db points when operated under matched input and output impedance 
conditions. It was therefore necessary to widen the band by using double- 
tuned circuits throughout. In fact such a design would be difficult to avoid 
because each tube has its own resonant cavities. By running a short length 
‘of waveguide or coaxial transmission line from the output cavity of one stage 
to the input cavity of the next, a double-tuned structure results auto- 
matically. The coupling from cavity to transmission line is made by means 


* These tubes were developed by A. L. Samuel and J. W. Clark. 
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of a window in the side of the cavity when a waveguide line is used, and by 
means of a wire loop, projecting within the cavity when a coaxial line is used. 
Amplifiers employing both types of coupling were designed and tested. The 
coaxial coupling type with links approximately one-half wavelength long 
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Fig. V-15.—A single stage of a coaxial-coupled amplifier using the velocity modulated 
tube of Fig. V-14. 


proved the easier to adjust. Variation of the coupling from zero to a maxi- 
mum is accomplished by rotating the plane of the coupling loop through 
90°. Each cavity is provided with a small plate of resistance film pro- 
jecting into the field in such a way that rotation of the plate about an axis 
lying in the plane of the plate will vary the resistive loading present in the 
cavity from zero toa maximum. Finally, the input and output cavities of 
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the first and last stages respectively are coupled to separate tuned cavities 
in order to provide double-tuned terminals for the amplifier. This results 
in the overall structure shown schematically in Fig. V-16. | 
_ The required bandwidth is obtained by a process of stagger-coupling the 
circuits so that the individual responses are as indicated schematically at 
the top of Fig. V-16. Because there are available continuously variable 
adjustments of tuning, loading, and coupling for each circuit, it is possible 
to obtain a very smooth and symmetrical overall response characteristic as 
shown in Fig. V-17.. The corresponding measured delay distortion charac- 
teristic is shown in Fig. V-18. For this type of tuning the output power is 
about .7 watt at 4 db of compression. More power is, however, obtainable 
if more compression is tolerable; and when used in the repeater for the 
transmission of FM signals, the amplifier is driven to an output of 1 watt. 

Since the circuits are of fairly high Q, the frequency characteritics of the 
amplifier are markedly affected by changes in temperature of the cavities. 
In order to minimize such detuning effects, the amplifier has been placed in 
a temperature controlled compartment. Since about 45 watts are dissipated 
at the collector of each tube by the high-voltage electron beam, the collectors 
must be cooled by a forced air blast. In order to keep the cooling system 
separate from the temperature control system, the collector ends of the tubes 
project through a wall of the temperature controlled compartment into an 
external air duct. Figure V-19 shows a complete r.f. amplifier with the 
cover of the temperature control box removed. This is the amplifier de- 
veloped for the New York-Boston circuit.* The electron-gun ends of the 
tubes face the reader. Short lengths of flexible coaxial cable couple the 
input and output of the amplifier to the associated equipment. 

Testing of Components and Repeater Amplifier—-Many special measuring 
and testing techniques had to be devised for the research work on each com- 
ponent. In addition, standard production tests had to be worked out and 
measuring equipment constructed. 

Impedances were measured at RF with standing wave detectors, and at 
IF with three fixed taps on a coaxial line. It was found that input and 
output SW ratios could be held to 1.7 db or less over the 10 megacycle band. 
Both point by point and swept oscillator methods were used at RF and IF for 
measuring amplitude characteristics. Particularly as a result of the de- 
velopment of swept oscillator measuring techniques it was found practicable 
to adjust each unit to --0.1 db amplitude variation over the 10 megacycle 
band. Noise figures!® of IF equipment were measured by the noise diode 
method.** 


* Developed by F. E. Radcliffe and R. C. Carlton. 

16 Loc. cit. 

38H. Johnson, “A Coaxial Line Diode Noise Source for U.H.F.”, R. G. A; Review, 
vol. VITI, No. 1, pp. 169-185, March 1947. 
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Fig. V-16,—Schematic diagram of four stage coaxial coupled amplifier using velocity modulated tubes. 
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A method*’ of measuring the phase and group delay characteristics of 
the repeater was worked out. This is a rather difficult measurement to 
make because the tolerable phase and delay distortion are very small due 
to the wide band of the system. Relative group delay through the band 
was measured with an accuracy of about -£0.001 microsecond. This cor- 
_ responds to an accuracy in relative phase shift of about +0.35°. The 
measured distortion agreed reasonably well with the distortion calcu- 
lated from the constants of the various circuits. These measurements and 
calculations showed that while the characteristics of an 8-link repeater 
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Fig. V-17.—Frequency response of RF amplifier. 


circuit will give acceptable television pictures, phase equalization will ' 
provide a definite improvement. Longer circuits will therefore require 
phase equalization. Equalizers were devised for both the IF and RF 
circuits which made it possible to equalize the group delay of each compo- 
nent of the repeater amplifier to 0.001 microsecond over the 10-mc band. 

An experimental repeater amplifier was set up so that the output could 
be fed through an attenuator to the input. It was then possible; to break 
the loop and make overall tests of delay, amplitude linearity, and transient 
response with IF measuring equipment. Transient response was measured 
by including a 2000-foot waveguide line between the input and output 
terminals and applying the circulated pulse testing technique.** This 

7D, H. Ring, “The Measurement of Delay Distortion in Microwave Repeaters”, 
elsewhere in this issue of the B.S. T. J. 


383A, C. Beck and D. H. Ring, “Testing Repeaters with Circulated Pulses”, Proc. 
I. R. E., vol. 35, No. 11, pp. 1226-1230, November 1947. 
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testing method permits a study of the shapes of rectangular pulses which 
have passed through a repeater amplifier many times, and thus simulates 
transmission through a relay system with many repeaters. Figure V-20 
shows an example of the results obtained from such tests made on the IF 
amplifier alone. The top row shows a 1.25-microsecond test pulse after 1, 


odie desc le lLa ke ae 
\eSaSSSSSRSEnetE 
Lee lll Gecsodece aneap eM a 























RELATIVE DELAY IN SECONDS 





4052 4054 4056 4058 4060 4062 4064 4066 4068 
FREQUENCY IN MEGACYCLES PER SECOND 


Fig. V-18.—Delay characteristic of RF amplifier. 


10, and 30 trips through the IF amplifier. Part of the distortion appearing 
after one trip through the amplifier was in the viewing oscilloscope. It 
was difficult to detect distortion from a single trip through the amplifier, 
but Fig. V-20 demonstrates how it increases with successive trips. The 
second row of pulse pictures shows the improvement achieved by adding 
phase equalizers to the circuit. It was particularly noticeable in this test 
that without equalizers the details in the shape of the pulse after 10 or 30 
trips were very sensitive to the exact value of the intermediate frequency. 
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When the equalizers were added, the pulse shapes shown in the second row 
were not greatly changed for variations of +3 megacycles or more in the 
intermediate frequency. : 2. 

In concluding this section it should be noted that when the various com- 
ponents were connected together to form an IF-type repeater amplifier, 
it was found that the amplitude and phase characteristics added as expected 
and resulted in a satisfactory amplifier with only a few millimicroseconds 
variation in delay and only a few tenths of a db variation in amplitude over 
a 10-megacycle band. Nevertheless, the equipment is very complicated. 
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AMPLIFIER WITH EQUALIZER 
Fig. V-20.—Results of circulated pulse tests on IF amplifier. 


A straight-through radio frequency amplifier repeater is still to be desired 
and, no doubt, further research will eventually produce such an amplifier. 


VI. THe Complete REPEATER* 


In the preceding section it has been pointed out how the various compo- 
nents that make up the repeater amplifier were added together without the 
introduction of additional distortion. The antennas, filters and amplifiers 
which go to make up one of the complete repeaters shown in Fig. II-1 were 
designed with the same ease of interconnection in mind. However, as will 
be discussed below, the length of the waveguide lines used for these connec- 
tions has an important bearing on the distortion introduced. 

The large amount of equipment in the repeater amplifier makes it 
desirable, from the maintenance standpoint, to locate the amplifier near 
the ground and to provide towers for the antenna where antenna elevation 


* Prepared by D. H. Ring. 
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is necessary. This means relatively long transmission lines between the 
antenna and the rest.of the repeater. - If exact termination of the line by 
the antenna impedance or by the filter input impedance were possible no 
distortion would be produced, but in general there will be a slight mismatch 
and a corresponding reflection of energy at each of these junctions which 
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Fig. VI-1.—Effect of long lines on the amplitude and delay distortion of a repeater. 


will produce variations in the amplitude and delay of the signal throughout 
the desired band. These variations may be greater than those produced 
anywhere else in the repeater. 

The type of distortion originating in this way is illustrated by Fig. VI-1. 
In this figure there is represented an antenna of impedance Z; connected 
by a line of length ¢ and characteristic impedance Zo, toa load, Z:. In actual 
practice this load is the impedance presented by the filter to the line from 
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the antenna. The variations produced in both the amplitude and delay 
of the signal currents are shown by the curves in the figure where relative 
change in the characteristics in question are shown as ordinates and values 


of the quantity leh are shown as abscissas. 

It will be seen from these curves that the amount of variation over a given 
band is a function of the line length. With short lines the amount of varia- 
tion over a band 10 mc wide can be kept very small, but with lines about 
150 feet in length there may be three or four full cycles of variations in 
a band 10 megacycles wide. Table A gives some typical values of the varia- 
tions. Values are given for those degrees of mismatch which would produce 
standing wave ratios of 1, 1.4 and 2 db at the junctions. It will be seen that 
for lines between 100 and 200 feet in length the variations in some cases are 
greater than the limits achieved in the multistage amplifiers and other com- 








TABLE A 
SWR, = SWR, in db (a = 0).......0005-- 1 1.4 2 
ry = re (SD) ives tewes 0.058 0.082 0.115 
Max. Amplitude Variation in db............ 0.058 0.116 0.233 
Max. Delay Variation in 10-9 seconds....... 
P= TOD feet shied Cote ce ales teens 1.33 2.67 5535 
€ = 200 fete cane ced taba wea wees 2.67 5.33 10.66 








ponents of thesystem. Furthermore, it is usually impractical to compensate 
or tune out these variations because they are functions of the electrical 
length of the transmission lines and subject to change with temperature, 
frequency, and other small mechanical and electrical changes. It would 
appear that with present techniques this may be one of the most serious 
sources of distortion in a long relay system. 


VII. CONCLUSION 


Various phases of our work on microwave repeater circuits have been dis- 
cussed. Such a circuit, made up of several repeaters as described in the © 
last section, may be looked upon as a four-terminal network having speci- 
fied amplitude and delay characteristics and should be suitable for the trans- 
mission of any signals for which the characteristics are adequate whether 
they be television signals or those of one of the various forms of multi- 
plex telephony. It is outside the scope of this paper to discuss the uses to 
which such a repeater circuit might be put or the terminal equipment that 
any particular service might require. 

The New York-Boston repeater circuit has been built to provide the ex- 
perimental field trials necessary to answer many remaining questions which 
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deal with the performance of an actual circuit. The results of these trials 
will be reported in a separate paper. 


ACKNOWLEDGEMENTS 


The work which has been discussed in this report of necessity involved 
the combined efforts of many individuals. In addition to those already 
mentioned, practically every other member of the staffs of the Deal and 
Holmdel Radio Laboratories has made valuable contributions to the work, 
as have also many in other departments of the Bell Telephone Laboratories. 

The development of the components of the New York-Boston circuit has 
been under the direction of Mr. Gordon N. Thayer. 

In particular, we wish to acknowledge the support and stimulating advice 
of Dr. R. Bown under whose general direction the work progressed. 


The Measurement of Delay Distortion in Microwave Repeaters* 


By D. H. RING 


Measuring equipment is described which is capable of measuring delay distor- 
tion of the order of 10-® seconds in a wide band microwave television relay 
repeater. Two measuring circuits are discussed. The first is a circuit for meas- 
uring the relative phase shift versus frequency from which the delay distortion 
may be computed. The second circuit gives the delay directly from a single 
measurement. The measuring equipment is designed to work in the intermediate 
frequency range from 50 to 80 megacycles, but by applying suitable conversion 
equipment measurements can be made at microwave frequencies. 


HE successful transmission of broadband television and pulse signals 

over any communication circuit depends upon the preservation of the 
complex wave shapes of the original transmitted signals. Fourier analysis 
tells us that a complex signal wave can be resolved into a spectrum of fre- 
quencies with certain amplitude and phase relationships. It is well known 
that the amplitude relationships of all essential frequencies in this spectrum 
must be substantially preserved. It is equally important that the phase 
relationships of the essential frequencies should be preserved. The in- 
stantaneous value of the received signal is the vector sum of the instanta- 
neous amplitudes of all the component frequencies. Therefore, if the 
relative phase of some frequency component is changed by 180° the sign of 
its contribution to the output is reversed, and it is clear that a closer ap- 
proximation to the original signal could be obtained by suppressing this 
frequency component rather than permitting it to contribute negatively © 
to the output. . . 

It can be shown! that the relative phase relations of the component 
frequencies in a complex signal wave will be preserved if the phase shift in 
passing through a circuit is a linear function of the angular frequency. 
That is 


Ba Te ae (1) 


where 7) is a constant and n an integer. Distortion of the transmitted 
signal will occur if J is not constant over the essential frequency band of 
the signal. We shall not be interested in distortion due to n not being an 
integer” since this case does not occur in carrier circuits where the phase 
shift at carrier frequency, rather than the phase shift at zero frequency, is 
* Presented at National Convention of I. R. E., New York City, March 4, 1947. 
1 Phase Distortion and Phase Distortion Correction, S. P. Mead, B.S.7.J., April 1928. 


199-201. 
2 Phase Distortion in Telephone Apparatus, C. E. Lane, B.S.T.J., July 1930, 494496. 


247 


248 BELL SYSTEM TECHNICAL JOURNAL 


the reference point for phase. Departure of 6 from the linear relationship 


given by (1) is the phase distortion in the circuit. 
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Fig. 1—Typical phase shift curves for various types of circuits 
a. A Transmission line terminated in its characteristic impedance. 
b. A single tuned circuit. 
c. Two tuned circuits with approximately critical coupling. 


The time of transmission® or the delay in passing through the circuit is 
obtained by differentiating (1): 


Delay = = Ty seconds (2) 


dg: 
dw 
Variation in Ty with frequency is the delay distortion in the circuit. 
Figure 1 shows some typical phase curves, and Fig. 2 shows the corre- 
sponding delay curves. In each figure curve (a) represents an ideal dis- 
tortionless circuit with linear phase and constant delay such as a simple 
transmission line. Curves (b) are obtained for single resonant circuits, 


’ To has also been called the group delay and the envelope delay. 
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and curves (c) for coupled double tuned circuits. It is felt that the delay 
curves of Fig. 2 are easier to interpret and give a better physical picture 
of the distortion resulting from phase variations than the phase curves of 
Fig. 1. Therefore, most of the following discussion will be in terms of 
delay rather than phase. 


xWo! 


6.0 T 


DELAY TIME,-T, IN SECONDS 





Fig. 2—Typical delay curves for various types of circuits 
a. A transmission line terminated in its characteristic impedance. 
b. A single tuned circuit. 
c. Two tuned circuits with approximately critical coupling. 


Thus far we have considered the general case and stated that the delay 
must be constant over the essential frequency band for distortionless trans- 
mission. It should be noted, however, that when delay distortion is pres- 
ent different types of signals may be affected differently. For instance, in 
the case of an amplitude modulated carrier of angular frequency wo, Fig. 
2, we note that if the delay curves have arithmetic symmetry about wo, then 
the sidebands at wo + Aw will suffer the same delay and therefore will add 
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in phase upon demodulation. If the delay distortion is not symmetrical 
about wo the sidebands will not add in phase upon demodulation, and the 
demodulated output will suffer both amplitude distortion and some delay 
distortion which differs from the delay distortion at both wo -+- Aw and w9 
— Aw. In the case of frequency modulation dissymmetry introduces 
harmonics in the demodulated output. A detailed discussion of these 
effects is beyond the scope of this paper, but we may note that in general 
a true picture of the delay distortion in carrier circuits is not readily ob- 
tained by observing the demodulated output and that an unsymmetrical 
delay distortion is particularly undesirable in carrier circuits. 


' PRINCIPLES OF DELAY MEASUREMENT 


Delay cannot be measured directly on a steady state basis with a single 
test signal in the simple manner in which amplitude response is measured. 
Instead, the phase shift through the unknown network must be measured 
at two adjacent frequencies and the delay, or slope of the phase shift, 
computed from the relation 


AB 

| ss oe (3) 
Figure 3 illustrates the computation of the average delay in the interval 
Aw from two phase measurements. . 

The steady state phase shift of an unknown network can be measured 
by using the basic circuit shown schematically in Fig. 4.4 This is essentially 
a bridge circuit in which the phase shift in the unknown is balanced by an 
equivalent calibrated phase shifter. The phase comparator is some kind 
of device which will give an indication of a known relationship between the 
phases of the signals arriving over the unknown path and the known refer- 
ence path. . 

The exact form of suitable components and circuit arrangements for 
applying the basic method of Fig. 4 to a particular delay measuring problem 
is largely determined by the order of magnitude of the delay to be meas- 
ured and accuracy desired. In the case of microwave television repeaters 
we are interested in video bands of the order of 5 mc wide. As a rough 
estimate we might say that the highest frequency in the band should not 
be shifted more than one quarter period from its normal phase position. 
In the case of linear delay distortion or a parabolic phase-frequency char- 
acteristic, one quarter period for 5 mc is 0.05 microseconds. In a repeater 
system with 50 repeaters this yields a tolerable systematic delay distortion 
of 10-° seconds per repeater. Therefore we conclude that in developing 


. or a of Phase Distortion, H. Nyquist and S. Brand, B.S.T.J., July 1930, 
526-527. 
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repeaters for this service an accuracy of better than 10-° microseconds in 
‘measuring the relative delay over a band of frequencies will bé desirable. 
The phase shift, AB in Fig. 3, which corresponds to a delay of 10-° seconds 
is a function of the measuring interval Aw = 27Af. If Af is small Ag will 
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Fig. 3—Factors involved in calculating the delay of an electrical circuit 
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Fig. 4—Basic circuit for measuring the phase shift in a network 


be small and difficult to measure. If Af is large the average slope measured 
will not be the true slope in the center of the interval. For a 5-megacycle 
video signal the intermediate and radio frequency bands of interest will 
be in excess of 10 megacycles wide. These considerations led to a choice 
of 1 megacycle as a reasonable value for Af. If Af = 1 megacycle a phase 
shift A8 = 0.36° will result from a delay of 10-9 seconds. 
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Two circuits for measuring delay distortion will be described. The 
first is a phase measuring circuit which, in principle, is an adaptation of 
Fig. 4 to practical operation at intermediate frequencies. The second 
circuit is a modification in which two frequencies differing by Af are fed 
through the circuit under test simultaneously in such a way that the differ- 
ence in the phase shift at the two frequencies is measured. This arrange- 
ment permits the calculation of the delay from a single measurement and 
is, therefore, a delay measuring circuit. | 
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Fig. 5—Schematic circuit for the precise measurement of phase shift in the interme- 
diate frequency range. 


PHASE MEASURING CIRCUIT 


Figure 5 shows a schematic diagram of a phase measuring circuit which 
has been found to be suitable for precision measurements of wide-band 
circuits in the intermediate frequency range of 50 to 80 megacycles. It 
is a double detection system with an intermediate frequency of one mega- 
cycle. The test signal oscillator and the beating oscillator are ganged to a 
single control and adjusted to track so that they maintain a difference of 
approximately one megacycle throughout their tuning range. The test 
signal is fed to a sliding contact on a section of air dielectric coaxial trans- 
mission line. The signal divides at this point to feed the test branch and 
. reference branch. 

The sliding tap on the coaxial line provides a high precision measuring 
phase shifter if the coaxial line is well terminated at each end. The rela- 
tive change A@ in phase of the signals at the two ends of the line when the 
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slider is: moved a distance Ax is two times the phase shift corresponding to 
‘the movement of the slider at the working frequency or 


ie ok a 

ne Besa 4 Ag = degrees | 

where c-is.;the velocity: of light. 

For Ax = 0.1-centimeter and f = 65 megacycles, AG is 0.156 degrees. 
The measuring branch and reference branch feed through the network to 

be measured and a balancing line respectively to identical detectors and 

one megacycle amplifiers. The amplifiers are connected to the plates of a 

cathode ray oscilloscope which is used as a phase comparator. A cathode 





leateaes betaine ‘a — aoe 





X=Acos wt 
YZASIN (wt +84) 
S=2A SINO 


IF A=3"AND @=0.1° 
THEN S ¥ 0.01" 


Fig. 6—Calculation of the sensitivity of a cathode ray oscilloscope as a phase com- 
parator. 


tay oscilloscope has the advantage that the phase comparison with this 
instrument is independent of either the relative or absolute amplitudes of 
the two signals. A straight line on the oscilloscope always indicates that 
the two voltages are in phase (or phase opposition) regardless of the relative 
amplitudes, which determine the slope of the line. The sensitivity, of 
ccurse, is a function of the amplitudes. Figure 6 illustrates how the sensi- 
tivity of an oscilloscope phase comparator can be calculated. If each 
signal alone produces a 6-inch deflection, then 0.1 degree phase difference 
produces an opening of the pattern of 0.01 inches, which is sufficient to be 
detected on a “‘rocking”’ or in-out basis. 

It is obvious that a circuit of this type measures the difference in the phase 
shifts of the two branches. The measured phase shift is, therefore, the 
absolute phase shift in the measuring branch less the phase shift in the ref- 
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erence branch. The balancing line shown in the reference branch of Fig. 
5 can be adjusted in length so that it balances out the average delay in the 
measuring branch. The measured remainder will then be the distortion 
in the measuring branch, since a good transmission line does not have 
phase or delay distortion. The use of a balancing line in the reference 
branch simplifies measurements by reducing the range of movement of the 
slider, and it greatly decreases the errors in the calculation of the delay 
distortion due to inaccuracies in the measuring interval Af. . 

The simplified diagram of Fig. 7 will be used to show how the delay of © 
the unknown network can be calculated. With the signal oscillator set at 
frequency /; the slider is adjusted until the signals at C and E are in phase 
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Fig. 7—Simplified circuit electrically equivalent to the circuit of Fig. 5. 


as indicated by a straight line on the oscilloscope, and the corresponding 
distance x, is measured. The relationships between the phases of the 
signals at the various points in the circuit are: 


gc = ba — bas — Oac 


= $4 — “2a -a) —8 @ 
be = b4 — Gav — hoz 
=, ert a == daft ; | (5) 
C ov 
"dc = oz 


_ where § is the phase shift in the unknown network at frequency f;. Solving 


for B, 


c 


B = oan |% + wy — a| © 
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The signal oscillator frequency is then changed to fe = fi + Af and the slider 
readiusted to the position x2 that again makes ¢@¢ = ¢z. Then as above 


B= ?ais| & epg | | (7) 


where 8’ is the phase shift in the unknown network at frequency fo. The 
delav in the network is, from (3) 





_f -8 | 
 QnAf (8) 
Substituting (6) and.(7) in (8) yields 
_ = : (£- «) + 2%. + — 2h (x2 — x) | (9) 
cL\v Af 


The first term in T is a constant of the set-up and can be dropped when the 
delay distortion only is desired. The second term js small and can often 
be neglected. The third term gives the major part of the difference between’ 
the delay of the reference path and the delay of the network. 

It has been found that the slider position can be easily reset to + 0.05 
centimeters for each frequency. This would mean a maximum error of 
+ 0.1 centimeter for the difference of two readings which corresponds to 
an error in T of about -+-0.4 X 10-® seconds. However, this reset accuracy 
will not be realized in overall accuracy unless a number of precautions are 
- observed in setting up the circuit of Fig.5. In order to avoid stray coupling 
the two branches of the system must be carefully shielded from each other. 
At least 60 db net attenuation must be provided between the detectors via 
the path through the balancing line, phase shifter and unknown network, 
and 40 db attenuation in the path via the BO supply lines. All attenuators, 
plugs, etc., must have voltage standing wave ratios of less than about 
1.015 and the detectors and amplifiers in each branch must have identical 
phase shifts over the range of variation of the IF due to imperfect tracking 
of the oscillators. 


DELAY MEASURING CIRCUIT 


The circuit of Fig. 5 is basically a phase measuring circuit. It can be 
rearranged as shown in Fig. 8 to yield a circuit that will read delay directly 
from a single setting of the slider. In Fig. 8 the signals from the two 
oscillators are both sent through the circuit to be measured, and both are 
sent through the reference branch. Any delay in either path will. alter the 
relative phases of the two signals in that path and this change in relative 
phase shift will be passed on to the beat note formed in the detectors. If 
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the delays in the two paths are different a relative phase shift proportional 
to the difference will appear in the two beat notes. This phase change can 
be measured with a phase shifter in the beat note circuit or with a phase 
shifter which varies the relative phase of the two signals fed to one branch 
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Fig. 8—Schematic circuit for the precise measurement of delay in the intermediate 
frequency range. 
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Fig. 9—Simplified circuit electrically equivalent to the circuit of Fig. 8 


as compared with the relative phase of the two signals fed to the other 
branch. Figure 8 shows how the coaxial slider can be used for the latter 
type of measurement. . 

The simplified diagram of Fig. 9 will be used to show how the delay of the 
unknown network can be calculated in terms of this circuit. It will be 
assumed that the electrical lengths of the paths AF and BF are equal so 
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that this length can be peparace as part of the balancing line £. Then at 
point E we have 





; dr = ba oe ar = ba _ a . 
be = bn — doe = bn — eh 


The primes indicate phase shifts at frequency fe. The phase $”z of the 
beat note, Af = f. — fi, at E is 


ox = bn — be = bs — bs — —Af (10) 


Similarly, at point D 


” , = a x 


oo = $3 — a 8 2r fila == *) = (Q’ _ B) (11) 


If. is adjusted so that the Af’s in the two branches are in phase at the 
oscilloscope then 


"pd = ob" x ‘or a (12) 
where a@ is the difference in the phase shifts in the two beat note circuits. 
Substituting (10) and (11) in (12) and solving for (@’ — 8) yields 


@ -) =22| 2! —sayp-far-g +2] 08) 


Cc 


The delay in the network is 


(Br =p). 
es QnAf 


a dle eae (2x — a) 
7 [4 Uap oe | 
The first two terms in (14) are independent of x and f; and therefore are 
of interest only for absolute measurements. Relative measurements may © 
be made without evaluating these constants. The last two terms are 
functions of x and yield the change in delay as f is varied. It is usually 


most convenient to adjust a.and ¢ so that the average delay in the meas- 
uring branch is given by (14) with 


(2x — a) = 0 (15) 


In general this will minimize the variation of the slider and make the delay 
distortion in the measuring branch roughly proportional to the slider 


(14) 
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movement.’ This is helpful in judging the effect of adjustments of the 
unknown network. The condition (15) corresponds to the optimum 
condition for the circuit of Fig. 5 which requires that the average delays of 
the two branches be equal. In the circuit of Fig. 8, the condition (15) can 
be realized by varying either f ora. Ifa is made zero and (15) is fulfilled 
by adjusting ¢, errors due to variations in Af are minimized. However, if 
Af is held sufficiently constant, the balancing line ¢ can be omitted entirely 
and a 360° variable phase shifter introduced in the beat note circuit of one 
branch to vary a. . 
The measuring interval Af is determined by the difference in the two 
signal oscillators. Since Af has an important influence on the measurement, 
oscillator tracking cannot be relied upon to maintain Af with sufficient 
accuracy. A one-megacycle crystal frequency checker has therefore been 
included in the equipment as shown in Fig. 8. A sample of the signal from 
one of the amplifiers is compared with the crystal oscillator and a trimmer 
on one of the oscillators adjusted for zero beat before measuring each 
point. When T is so large that a balancing line is not practicable, as in 
the case of loop circuits including radio paths or long transmission lines, 
Af must be held constant to about 1 part in 10°. This has been accom- 
plished in a modification of this equipment built by Messrs. W. J. Alber- 
sheim and J. P. Shafer of these laboratories, by deriving the two measuring 
frequencies from a crystal oscillator. 
Tn order to obtain the absolute value of 7, the constants ¢, v, anda must be 
known. The value of {may be found from the physical length of the line; a 
can be measured by measuring the movement of the slider required to rebal- 
ance the circuit when the connections to the two detectors are reversed. The 
absolute delay is particularly sensitive with respect to the difference be- 
tween 2x and a. If a is measured as accurately as possible, then the exact 
setting of the slider corresponding to condition (15) can be found by re- 
versing the connections of the slider at points A and B in Fig. 9. In this 
way a reference value of « may be found which will yield accurate results 
in spite of a small error in the value of a. 
Successful operation of this circuit depends on low standing waves 
- throughout and upon sufficient padding for satisfactory isolation of the 
oscillators. It will be noted that in the analysis it was assumed that f; 
reached the slider via.the path AC, Fig. 9. There must be an attenuation 
for fi in the path AFBC sufficient to render the signal traversing this path 
negligible. Similar unwanted paths exist for fe from B to C and also for 
fi and fe to point F. Attenuation inserted as shown in Fig. 8 can be ar- 
ranged to make these unwanted signals negligible. 
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Rapio FREQUENCY MEASUREMENTS 


The frequency range of a particular measuring equipment using the 
circuits of Fig. 5 or Fig. 8 is determined only by the range of the oscillators 
and the range over which the plugs and jacks and attenuators operate 
with sufficiently low reflection coefficients. While this range is greater 
than .the IF range encountered in microwave repeaters, it does not include 
the actual microwave frequencies. However, it has been found that 
microwave components can be measured satisfactorily by using the circuit 
shown in Fig. 10. In this circuit the measuring equipment is operated at 
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. Fig. 10—Method of using the intermediate frequency measuring circuit of Fig. 8 for 
the measurement of radio frequency networks. 


IF, and the reference branch is unaltered. The measuring branch signal 
is fed first to a converter where it is beat up to the desired microwave 
frequency. The converter is followed by a filter which eliminates the 
beating oscillator frequency and one of the beat frequencies. The filter 
output is then applied to the microwave component under test. The 
output of this is fed to another converter and converted back to IF by 
beating with the same oscillator that was used in the first converter. This 
process removes any variations in phase due to variations in the beating 
oscillator if the connections from the beating oscillator to each converter 
are of equal electrical length. . 

Since considerable extra equipment is included in the measuring branch 
in Fig. 10, it will usually be necessary to make a calibration run with the 
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device to be measured removed in order to eliminate any possible delay 
distortion in the converters and filter. Figure 10 uses the measuring 





ices 


Fig. 11—Photograph of apparatus used in the measuring circuits shown in Figs. 5 
and 8. 


circuit of Fig. 8 rather than that of Fig. 5 because it was found that small 
variations in the transit time in microwave amplifiers cause small variations 
in the phase shift which are the same at all frequencies in the band. These 
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variations cause changes ‘in the relative phase of the two successive meas- 
urements required when using the circuit of Fig. 5 which do not represent 
changes in the delay. In effect the circuit of Fig. 8 makes the two phase 
measurements simultaneously, and thus eliminates effects due to variations 
of phase with time. 
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Fig. 12—Measured curve of the relative delay of an intermediate frequency amplifier 


RESULTS 


Figure 11 shows a photograph of the delay measuring equipment which 
has been described. With the aid of patch cords and switches that have 
been included in the equipment this apparatus can be set up according to 
either Fig. 5 or Fig. 8. The ganged oscillators are on the panel above the 
oscilloscope box. The box contains the dividing attenuators, detectors, 
amplifiers and oscilloscope. A separate power supply is required for the 
oscillators and the output stages of the amplifiers. A number of different 
lengths of flexible coaxial cable are mounted on the panel above the oscil- 
lators and arranged so that various lengths of balancing line can be con- 
veniently obtained by patching. The coaxial phase changer can be seen 
on the shelf in front of the oscilloscope. 
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Figure 12 shows the measured relative delay of a typical intermediate 
frequency amplifier. This amplifier has a substantially flat amplitude 
response over a band of about 12 mc centered on 65 mc. The delay dis- 
tortion over the 10 megacycle band from 60 to 70 mc is about 10 X 10-® 
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Fig. 13—Measured relative delay of an experimental delay equalizer 
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Fig. 14—Measured relative delay of the amplifier of Fig. 12 plus an equalizer 


seconds. Figure 13 shows the measured relative delay of a bridged T phase 
equalizer which can be used to compensate for the distortion in the amplifier. 
Figure 14 shows the measured relative delay for the amplifier of Fig. 12 
and an equalizer measured together. The equalizer has reduced the delay 
distortion over the 10 mc band from about 10-8 to 10-® seconds. These 
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measurements were made with an early model of the measuring equipment. 
Smoother curves are obtained with the apparatus shown in Fig. 11. 

In Fig. 15 the top row shows the distortion of a square top pulse by the 
amplifier of Fig. 12 for 1, 10, and 30 trips through the amplifier without the 
equalizer. The lower row shows a similar set of pictures of the pulse when 
the distortion was reduced by a phase equalizer as shown in Fig. 14. The 
improvement due to the elimination of phase distortion is clearly illustrated. 
These pictures were obtained by the circulated pulse’ technique which 
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Fig. 15—Oscillograms showing the improvement in the square wave response of the 
amplifier of Fig. 12 obtained by delay equalization. The numbers above the traces indi- 
cate the number of times the test pulse has passed through the amplifier and equalizer. 


permits the observation of a pulse after it has passed a number of times 
through the same amplifier. 


CONCLUSIONS 


Two measuring circuits have been described which are suitable for 
measuring the small variations in relative transmission time which are 
present in wide band microwave repeaters. If sufficient care is exercised 
in setting up these circuits an accuracy of better than +10~* seconds can 
be realized in relative delay measurements. The circuit of Fig. 5 measures. 
the relative phase shift as a function of frequency. It has the advantage 
of requiring less signal power and fewer important parameters for absolute 


5 Testing Repeaters with Recirculated Pulses, A. C. Beck and D. H. Ring. Proc. 
I.R.E. Nov. 1947, 1,226-1,230. 
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delay measurements. The circuit of Fig. 8 measures delay directly with.a 
single measurement and has the advantage of ignoring uniform phase vari-' 
ations with time. It is useful for making relative measurements on circuits 
with long constant delay times. A significant improvement in the square. 
wave response of carrier amplifiers has been obtained by as delay. 
equalizers based on measurements made with this equipment. | 


Frequency Shift Telegraphy—Radio and Wire Applications* 
By J. R. DAVEY and A. L. MATTE 


Frequency shift telegraphy is described and compared with amplitude modula- 
tion telegraphy under various conditions found in radio and wire transmission. 
Experimental data are given to demonstrate the influence of various design fac- 
tors on the over-all performance under these conditions. It is shown that the 
most outstanding characteristic of the frequency shift method is its ability to 

. accept large and rapid changes in signal amplitude. Frequency shift telegraphy 
thus proves to be of great advantage for use in the H.F. radio range. Frequency 
shift telegraphy also shows an advantage over amplitude modulation telegraphy 
with respect to noise. For applications where the level variations are small or 
slow the advantage of the frequency shift method over amplitude modulation 
is relatively small. 


INTRODUCTION 


URING World War II, single-channel and multichannel frequency- 
shift radio telegraph systems proved of the utmost importance in pro- 
viding the Allied Powers with a world-wide automatic printing telegraph 
network for handling with precision, secrecy and dispatch the unprecedented 
volume of traffic engendered by a war of global extent. It is expected that 
the next few years will witness a greatly expanded application of this 
method of operation by commercial telegraph companies and others in- 
terested in long distance telegraphy. 

Frequency Shift carrier telegraphy (FS) may be applied to any carrier 
telegraph circuit, but, as will appear below, it provides particularly striking 
advantages in H.F. radio transmission. For some other radio frequency 
ranges and for wire line operation the conditions are such as to limit the 
advantages of the FS method. The main advantages of the FS over the AM 
method are a greater ability to accept rapid level changes, which results in 
better stability and lower distortion, and an improvement in signal-to-noise 
ratio, which permits a reduction in carrier amplitude. It is therefore of 
particular importance where automatic printing is desired over H.F. radio 
circuits. When it is necessary to transmit through very high noise levels, 
low speed AM signaling with aural reception of an audio beat note remains 
the superior method. 

FS is a form of frequency modulation in which signaling is accomplished 
by shifting a constant amplitude carrier between two frequencies represent- 
ing respectively the marking and spacing conditions of the telegraph code. 
Frequency variations in FS telegraphy correspond to amplitude variations in 


*Published in A.J.E.E. Transactions, Vol. 66, p. 479, 1947. 
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AM telegraphy (CW); thus the signal transitions in FS are represented by 
frequency-time transients, while in the AM case they are amplitude-time 
. transients. Since AM telegraph is the more common system, a discussion of 
the FS method involves numerous comparisons between the two systems. 
The merits of a telegraph system must be judged on its ability to combat 
the various adverse conditions encountered in the transmission medium and 
in the terminal apparatus. In general these adverse conditions involve 
variations in amplitude, frequency, and phase of the signals and the presence 
of extraneous signals and noise. 

In the course of the development of a number of FS radio teletypewriter 
systems, a large amount of information concerning the characteristics and 
design parameters of such equipment has been obtained. It is the purpose 
of this paper to abstract therefrom selected data which will furnish a step- 
by-step comparison of the FS and AM methods. Typical terminal arrange- 
ments are described and the effects of varying certain design factors are 
illustrated by experimental data. Although the material presented applies 
largely to H.F’. radio telegraph, much of it is of a general nature and with 
proper interpretation applies to other frequency ranges and transmission 
mediums and to cases in which the telegraph modulated carrier may bea 
sub-carrier or one of several sub-carriers. 


GENERAL DISCUSSION 


Sideband Energy Distribution 


The difference between FS and AM signals as regards distribution of side- 
band amplitudes is illustrated by the following two equations for a carrier 
of frequency w/2r modulated with unbiased square wave dots of fre- 
quence p/2r. 

For AM (On-off) keyed carrier of ahity amplitude! :— 


= 0.5 cos wi foe ~ Icos (w + p)t + cos (w — p)il 
S = ede Gn Bale costa alll 


42 


Se [cos (w + 5p)t + cos (w — Sp)é]...etc. (1) 


For FS keyed carrier of unity amplitude?:— 


2m yt 5 (™ 
e = — | — sin | — ] cos wi 
- a | m 2 
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= — sin (27) (cos (w — 2p)t + cos (w + 2p)é) ~: 


mr 


= —— cos (2) (cos (w — 3p)t — cos (o + 3p)t) 


frequency shift 


where m is the deviation ratio = ———>—__>——_ 
2X signaling speed 


Typical sideband amplitudes calculated from these formulas are shown 
in graphical form in Fig. 1. In the case of FS keying, the relative. ampli- 
tudes of the sidebands vary considerably as the amount of frequency shift 
is changed. For miscellaneous signals these line spectra do not exist but 
they do indicate the general distribution of energy over the band for a given 
signaling speed. 


Methods of Modulating the Carrier 


In AM telegraphy the carrier is usually modulated by simply interrupting 
it for the spacing condition. This is sometimes referred to as “on-off’’ 
keying. For low power and low frequencies the carrier may be keyed di- 
rectly by electrical contacts. A more universally applicable method is to 
use vacuum tubes or other nonlinear elements to effectively interrupt the 
carrier. In some cases it is practical to start and stop an oscillator source 
of carrier. 

The usual radio telegraph transmitter consists of an oscillator followed by 

_a number of cascaded stages of amplification and frequency multiplication 

arranged to reach the desired output frequency and power. For on-off 
keying the carrier is usually interrupted by suitably varying the plate or 
grid voltage of one or more of the stages. 

There are two general methods of obtaining frequency modulation: (a) 
The frequency of an oscillator may be modulated directly by suitably vary- 
ing the frequency-determining circuit, (b) the output of a constant-frequency 
oscillator may be shifted in phase at such rates of change as to produce indi- 
rectly the desired frequency variations. In the latter case the marking and 
spacing intervals of an FS signal would be formed by periods of constant 
rate of phase change versus time. Square wave reversals would therefore 
require a triangular shaped wave of phase versus time. Since the transmis- 
sion of long periods of steady mark or space would therefore involve huge 
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phase swings, the modulator would have to be able to produce a steady phase 
rate of charge. A phase modulator capable of performing in this way while 
not producing undesirable phase discontinuities at the signal transitions 
becomes rather impractical. For this reason FS telegraphy usually utilizes 
the direct frequency modulation method. This may conveniently be accom- 
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Fig. 1.—Amplitude of sideband components for (a) square-wave amplitude modulation 
(b) square-wave frequency modulation. 


plished by the use of a reactance modulator which, by injecting a reactive 
component of current into the tuned circuit of the oscillator, varies the 
resonant frequency thereof. Sucha modulator may be made linear so that a 
frequency shift proportional to the input voltage to the reactance modulator 
is obtained. . 

To apply FS telegraph signals to a radio transmitter the regular exciter 
oscillator is either replaced or modified by an arrangement providing a 
source of R.F. excitation that can be shifted in frequency in accordance 


FREQUENCY SHIFT TELEGRAPHY 269: 


with the telegraph signal. All the stages are operated with full R.F. excita- 
tion continuously to produce a constant amplitude carrier. 

As a matter of expediency frequency shift keying has sometimes been 
provided by switching between two independent sources of carrier current 
separated in frequency by the desired shift. In such a case the frequency 
transitions involve sudden phase discontinuities of random values. This 
results in the instantaneous frequency swinging considerably outside the 
steady-state mark and space frequencies. If the band is wide, such 
as is the case in a radio transmitter, there results a very broad sideband 
radiation capable of causing severe interference to adjacent channels. If 
the band is narrow, as might be the case where sending filters are 
employed, the interference is eliminated but the amplitude transients result- 
ing from the sudden phase SHitts are capable of producing considerable 
distortion. 


Restriction of Transmitted Band 


As seen above, square-wave modulation results in a wide spread of side- 
band components which are of sufficient amplitude to interfere seriously 
with adjacent channels unless greatly attenuated. The transmitted band 
may be restricted either by the use of a band-pass filter centered about the 
carrier frequency or by a low-pass filter to suitably shape the modulating 
wave form. Band-pass filters are usually used if the power level is low and 
the frequency low enough to permit suitable filter construction. For multi- 
channel systems the use of band-pass filters also permits efficient paralleling 
of the transmitting channels. For radio transmitters with a transmitted 
power measured in kilowatts, and with a frequency of several megacycles 
which is frequently changed to suit best the prevailing conditions, shaping 
of the modulating wave is the more practical method of restricting the trans- 
mitted band. 

Insufficient attention has been given in the past to the envelope shape of 
the signals from on-off keyed radio transmitters. With the ever increasing 
crowding of frequency assignments it becomes more and more important to 
restrict the emission of unnecessary sidebands arising from keying. The 
envelope shape in on-off keying may be controlled by properly shaping the 
modulating grid or plate voltage wave. It is important that the stages 
following the keyed stage or stages be nearly linear, otherwise the wave 
shaping will be largely destroyed. In the case of frequency shift keying, 
on the other hand, the wave shaping is preserved after passage through class 
C amplifier or multiplier stages, and these may be operated for maximum 
efficiency. ‘The greater ease of producing and maintaining the desired wave 
shaping, so necessary for close frequency spacing of channels, is one of the 
outstanding advantages of frequency shift keying. 
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APPARATUS 


Typical FS Exciter for Radio Telegraph 


A typical FS exciter arrangement such as is often used with radio telegraph 
transmitters is shown in Fig. 2. A d-c. telegraph wave, after suitable shap- 
ing, is caused to frequency-modulate an intermediate frequency of 200 kc. 
which, in turn, amplitude-modulates a radio frequency from a crystal- 
controlled oscillator. The upper sideband of this latter modulation is an 
FS signal and is selected and amplified sufficiently to drive the first amplifier 
or multiplier stage of the transmitter. The 200-kc. oscillator is frequency- 
modulated by a reactance modulator which, by feeding a leading or lagging 
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Fig. 2.—Block diagram of a typical FS transmitter. 


quadrature component of current into the oscillator tuned circuit, decreases 
or increases the frequency. By operating the reactance modulator within 
its linear range the frequency shift wave form is made the same as the d-c. 
telegraph wave form into the modulator. A d-c. amplifier stage, designated 
“keying circuit”, is provided to furnish a modulating wave effectively iso- 
lated from amplitude and wave front variations of the incoming telegraph 
signals. The d-c. telegraph signals may be polar or neutral and are often 
obtained from a tone demodulator unit which allows keying from a remote 
point by V.F. telegraph. The amount of frequency shift is adjusted by an 
amplitude control in the quadrature feed-back path to the 200 kc. oscillator. 
The shift may thus be varied continuously, or in definite steps to allow for 
subsequent frequency multiplications, by suitable attenuation controls. 
Controlling the shift in this manner keeps the instabilities of the reactance 
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modulator a constant percentage of the frequency shift, which would not 
be the case if the shift were adjusted by varying the amplitude of the modu- 
lating wave. The use of a balanced instead of an unbalanced reactance 
modulator minimizes variation of the mean frequency and also allows the 
shift to be varied without affecting. the mean frequency. 

The frequency-shift signal transitions are wave-shaped, to restrict side- 
band radiation, by means of a low-pass filter in the d-c. telegraph signal path 
to the reactance modulator. The low-pass filtering is made adjustable to 
accommodate a range of signaling speeds. Frequency-versus-time wave 
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Fig. 3.—Frequency shift keyer output wave forms. (a) Low-pass filtering adjusted 
to produce similar wave shapes at dotting speeds of 60, 100, 200, and 240 cycles. (b) 200 
cycle phase modulation superimposed on a 23 dots per second signal. 


forms from an exciter of the type described are shown for several keying 
speeds in Fig. 3a. The effect of wave shaping on the sideband components 
in the R.F. output of such an exciter is shown in Fig. 4. 

Phase modulation may readily be added to the signal in this type of exciter 
by superimposing the desired sine wave modulating frequency on the tele- 
graph signal wave input to the reactance modulator as indicated in Fig. 2. 
Figure 3b shows the keyer output wave form with superimposed phase 
modulation. The use of this type of phase modulation is considered later. 

To obtain optimum results in FS radio telegraph transmission and to 
allow close spacing of channels, a high degree of frequency stability is neces- 
sary. An over-all frequency stability of +100 cycles is desirable in a system 
using a value of frequency shift between 500 and 1000 cycles. A frequency 


272 BELL SYSTEM TECHNICAL JOURNAL 


shift exciter of the type described above, with the crystal oscillator and 
200 ke. FS oscillator located in a temperature-controlled oven, usually has 
a frequency stability such that the mean R.F. carrier frequency may be 
held to within --50 cycles up to frequencies of 20 mc over ordinary periods 
of operation on any one frequency. One of the advantages of this type of 
exciter is that small inaccuracies in crystal frequencies may be compensated 
for by adjusting the mean frequency of the 200 kc. oscillator. 


50 
45 A 
MAXIMUM | 
40 SHAPING 
Je ~a/ 
Wig 
ptr 
O35 
we 
96 
Za 
chad 
ot 7 
PR2s ofS 
ao / MINIMUM 
3220 fi SHAPING 
fa) H 
2315 I 
ag I 
on I 
wn” 


re) 


‘ 
Ui 





10 8 6 4 2 C 2 4 6 8 10 
ORDER OF SIDEBAND 


Fig. 4.—Effect of shaping FS transitions on amplitude of radiated sidebands. 100 
dots per second and 500 cycle frequency shift. The maximum and minimum wave shap- 
ing aon correspond to those used with 60 and 240 dots per second respectively in 
Fig. 3(a). 


Receiving Terminal Arrangements 


Typical receiving terminal arrangements are shown in the block diagrams 
of Fig.5. Up to point ‘‘a” in the arrangements the FS and AM systems are 
identical, being of the usual H.F. superheterodyne type. The output from 
the second frequency-conversion stage may be either in the audio range or 
at a considerably higher frequency such as 50 kc. Following the second 
converter is a band-pass filter (shown at “b” in Fig. 5) which determines 
the final over-all band width before demodulation. The two systems differ 
only in the method of demodulation. The AM (on-off) signals are amplified 
and rectified to give a d-c.telegraph signal. The FS signals are amplitude 
limited and passed through a frequency discriminating network and then 
rectified to give a d-c.telegraph signal. Beyond this point the two systems 
are again identical. The d-c. signals pass through a low-pass filter to remove 
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carrier ripple and higher frequency noise components and are then amplified 
to a suitable level to operate automatic recording or printing apparatus. 
The d-c.signals may also be used to modulate an audio frequency so as to 
pass the signals to a remote point by multichannel voice-frequency carrier 
telegraph methods. 
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Fig. 5.—Block diagram of a typical receiving arrangement for either AM or FS tele- 
graph signals. 


,The radio receiver portion of the terminals up to point “a” should be 
designed to have low noise and good selectivity. Extreme H.F. oscillator 
stability is necessary for either system if narrow band width operation is to 
be maintained without constant attention. An over-all frequency stability 
of +50 cycles is desirable for the receiving terminal over a period of 6 to 
8 hours. Sufficient selectivity and amplifier capacity should be provided 
at all points to prevent overloading by unwanted signals or loss of automatic 
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gain control. In the following discussion those portions of the terminals 
beyond the second frequency converter will be given major attention. 


Experimental Transmitting and Receiving Arrangements 


For the laboratory transmission studies described in the following sec- 
tions the transmitter and receiver were located nearby and connected to- 
gether by means of an amplitude modulator and associated with various 
sources of noise designed to simulate quantitatively and under controlled 
conditions the variations which would be encountered in the actual medium. 

Throughout the tests 7.42 unit start-stop signals were used unless other- 
wise stated, and the speed was 60 words per minute (23 dots per second). 
Their peak distortion and bias were measured on a cathode-ray tube tele- 
graph distortion measuring set. 

An exciter of the type shown in Fig. 2 was used as a source of signals. 
A frequency of 6.4 mc. was employed, with the radio receiver connected to 
the exciter output through an amplitude modulator. This modulator was 
an electronic circuit permitting amplitude modulation of a frequency-shift 
signal to produce unequal mark and space amplitudes. This modulator 
was also used to amplitude-modulate a single frequency for the AM portions 
of the measurements. 

A temperature-limited diode together with a two-stage tuned amplifier 
was used as a source of thermal noise centered around 6.4 mc. A polar 
relay driven by 60-cycle a-c and arranged to produce sharp polar impulses 
from the discharge of small capacitances connected to its contacts was used 
as a source of impulse noise. The noise level was adjusted by an attenuator 
and mixed with the 6.4 mc. carrier of the exciter. A minimum amount of 
wave shaping was used, so that the modulation may be considered as having 
been essentially square-wave. 


Receiving Arrangements 


The experimental data submitted in the following discussion was obtained 
from reception through a laboratory setup essentially like that shown in 
Fig. 5.. The radio receiver proper was a commercial type of H.F. super- 
heterodyne. ‘The output of the second frequency converter was in the audio- 
frequency range, which enabled the use of various band-pass filters at ‘‘b”’ 
of the type used in voice-frequency telegraph systems. The amplitude 
limiter was effective over an imput range of —60 dbm* to above +20 dbm. 
The pass-band characteristics of the radio receiver and of the several band- 
pass filters used in position “b” are shown in Figs. 6 and 7 respectively. 
Unless otherwise stated, the frequency shift signals were centered about 


* The symbol dbm signifies “db referred to one milliwatt’’. 
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Fig. 7.—Attenuation versus frequency characteristics of bandpass filters shown at 
(b) in figure 5, 
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Fig. 8.—Linear discriminator characteristic. 
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2125 cycles and were demodulated by a linear discriminator centered about 
2125 cycles as shown in Fig. 8. The characteristics of the low-pass filtering 
are shown in lig. 9. ‘These low-pass filters were adjusted by oscillographic 
observation of the signal wave form and had cut-off characteristics giving 
very little characteristic distortion’. The d-c. amplifier was a high-gain 
nonlinear type designed so as to have a square-wave output having transi- 
tions established by the passage of the demodulated voltage wave through a 
narrow amplitude range. The amplitude and wave front slope of the de- 
modulated wave thus had no effect on the output wave form and could not 
affect the distortion measuring equipment. | 
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Fig. 9.—Attenuation versus frequency characteristics of low-pass filters. 


EXPERIMENTAL RESULTS 


Band Width Before Demodulation 


The band width before demodulation determines the amount of noise and 
interference which is to be accepted along with the desired signal and thus 
largely determines the signal-to-noise condition at the antenna at which 
the system fails to receive intelligence. The band width at this point 
(point “‘a”’ in Fig. 5) also limits the signaling speed capabilities of the system. 
Tn the following experimental data the values of band width were measured 


between the points of 6 db loss above that at midband. 


Effect on Signaling Speed 


For both methods of signaling considered here the band width must be 
at least twice the maximum signaling speed in dot-cycles per second but it 
is found that signal distortion rises rapidly for a band width less than three 
times the maximum signaling speed, and that a factor of at least four times 
is indicated for a system which is to have reasonably low distortion with any 
margin of safety. The signaling speed capability of a given band width is 
nearly the same for FS signals as for AM (on-off) signals. In Fig. 10 is 
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shown the over-all signaling frequency response to FS and AM signals for 
a nominal band width at point “‘a” of about 740 cycles. It will be noted 
that the FS method is but slightly inferior and that both systems fail at a 
frequency of approximately one-half the band width. 


Effect on Noise: 


The effect of thermal noise on distortion for bandwidths of 115, 230, and 
740 cycles is shown in Figs. 11 and 12. The rms. noise-to-carrier ratios 
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Fig. 10.—Overall frequency response of a 740-cycle band to AM and FS signals. 


PEAK DISTORTION IN PER CENT 


Sp op oS 
-20 -18 -16 -14 -l2 -10 - -6 -4 -2 ) 2 4 6 
- RMS NOISE-TO- CARRIER RATIO IN DECIBELS 





Fig. 11.—Peak distortion versus thermal noise for FS transmission—80-cycle cutoff 
low-pass filter. 


indicated in the figures were measured in the 3300-cycle band of the radio 
receiver in all cases. ‘The actual noise-to-carrier ratio existing in the trans- 
mission band used was lower and may be obtained from the following table. 
The values of correction are from rms. measurements. 


Pass BAND AT Pornt ‘‘a” _ CORRECTION TO BE MADE 
1920 cycles —2.3 db 
740 —6.5 
230 —11.6 


115 a —14.3 
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For AM signals, Fig. 12, varying the bandwidth has little effect when the 
telegraph signal distortion is less than 15%. Although the wider bands 
accepted more noise power, this added noise merely produced high-frequency 
noise components which were removed by the low-pass filter. This added 
noise does, however, cause the peak noise to exceed the signal amplitude at 
a lower noise-to-carrier ratio and cause failure before that for a narrower 
band condition. 

In the case of FS signals, Fig. 11, changing the bandwidth and the fre- 
quency shift simultaneously, and in approximately the same proportion, 
alters the whole distortion characteristic. At low noise levels a wider band 
with a greater frequency shift gives an improved signal-to-noise condition. 
However, as the noise level is increased the wider band causes the peak noise 
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Fig. 12.—Peak distortion versus thermal noise for AM transmission—80-cycle cutoff 
low-pass filter. 


PEAK DISTORTION IN PER CENT 


to exceed the carrier at a lower noise level than witha narrower band. Thus 
a change to a wider band gives less distortion at low noise levels and more 
distortion at high noise levels. This results ina much more sharply breaking 
distortion characteristic for the wider band. This behavior is typical of 
frequency modulation systems in general. 

Although the noise actually passed by the 740-cycle band filter was 
approximately 8 db above that passed by the 115-cycle band filter, the dif- 
ference in the failure points (35-40% distortion) for the two bandwidths: 
will be seen to be only about half this amount in db for both AM and FS. 
This phenomenon is typical of carrier telegraph systems when compared 
at the same signaling speed. This means that it is not particularly beneficial 
to decrease band width to obtain lower distortion under high noise condi- 
tions. The main reason for narrow bands is for more economical use of 
frequency space. 

As to the comparison between FS and AM, the FS method has an advan- 
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tage of 2.5 to 4.5 db at a distortion of 35% to 40%, corresponding to the 
selector failure point of the usual teletypewriter. From a signal-to-noise 
standpoint it is thus seen that the gain in changing to the FS method is 
approximately equal to the resulting increase in average transmitted power 
of about 3 db. A comparison at a lower distortion such as 15% shows 
an advantage of 4 to6 db. At a still lower distortion the 740-cycle band, 
because of the higher deviation ratio, shows an improvement of over 10 db. 
In this region the slopes of the curves make accurate comparisons impossible 
due to the masking effect of other sources of distortion. These large im- 
provements at low noise levels are similar to those associated with wide- 
band FM broadcast systems. However, in carrier telegraph transmission 
the criteria are so different that the difference between a nearly perfect 
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Fig. 13.—Peak distortion versus impulse noise for FS transmission—80-cycle cutoff , 
low-pass filter. 


circuit and one of small distortion is not of great importance except when a 
large number of telegraph sections are to be operated in tandem. From a 
practical standpoint the improvement in signal-to-noise is not more than 
about 6 db for equal band widths. . 

In Figs. 13 and 14 similar characteristics are shown for the case of impulse 
noise. ‘The noise level values of these curves are purely relative since no 
attempt to measure the peak noise was made. The comparisons between 
AM and FS, and between different band widths, agree closely with those 
for thermal noise. 


Effect of Limiter on Signal-to-N oise Ratio 


The limiter in the FS system is a high-gain nonlinear amplifier which 
delivers to the frequency discriminating networks an essentially square wave 
having transitions coinciding with the passage of the instantaneous voltage 
of the input carrier signal through zero. ‘The limiter thus passes only the 
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frequency or phase changes of the signal. Noise voltages which are small 
compared to the signal cause approximately linear phase modulation of the 
signal and this is passed through the limiter. The amount of frequency 
deviation thus imparted to the carrier by a given component of noise is 
proportional not only to its amplitude but also to its frequency separation 
from the carrier. This gives rise to the so-called “triangular noise spec- 
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Tig. 14.—Peak distortion versus impulse noise for AM transmission—80-cycle cutoff 
low-pass filter. 
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Fig, 15.—Effect of the limiter on distortion versus thermal noise for FS transmission— 
740-cycle band, 350-cycle frequency shift, 80-cycle cutoff low-pass filter. 


trum’ when a linear frequency discriminator is used. If the limiter is 
removed from the frequency-shift terminal the noise components in phase 
with the carrier as well as those in quadrature therewith are allowed to reach 
the frequency discriminating network. For a balanced type discriminator 
this increases the demodulated noise for small amounts of noise about 3 db. 
In Fig. 15 is shown the effect of removing the limiter from the circuit. The 
limiter is seen to have little effect on the failure point, but an improvement 
of 2 to 4 db is shown in the 5 to 12% distortion region. 
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Other beneficial effects resulting from the use of a limiter are discussed 
below under “Level Variations’’. 


Demaditate of Frequency Shift Signals 


It is desirable that the frequency discriminating network be of the bal- 
anced type having two branches allowing differential combination of the . 
two rectified outputs. This minimizes the response to amplitude modula- 
tion not eliminated by the limiter. Two general types of networks have 
been in common use for FS telegraph. One consists of two bandpass filters 
centered about the mark and space frequencies respectively and effectively 
dividing the total band into halves. The other consists of a two-branch 
network each branch of which has a varying amplitude characteristic extend- 
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Fig. 16.—Effect of magnitude of frequency shift on distortion versus thermal noise in 
FS transmission—740-cycle band and 80-cycle cutoff low-pass filter. 


ing over the complete transmission band and usually well beyond. The 
amplitude-versus-frequency -characteristics of these two branches have 
opposite slopes and are of such shape that differential combination of their 
rectified outputs results in an approximately linear voltage-versus-frequency 
curve, passing through zero at midband. (Fig. 8) | 

In Fig. 17 is shown the characteristic of a two-bandpass-filter type of dis- 
criminator which was used in early frequency shift terminals. The charac- 
teristic is fairly flat near the mark and space frequencies so that this type 
of discriminator does not produce a triangular noise spectrum. In Fig. 18 
is shown the type of discriminator characteristic obtained by the use of two 
narrow bandpass filters. In this case there is no broad flat region around the 
mark and space frequencies and an intermediate type of characteristic (ap- 
proaching the linear type) is obtained. 

With the linear type of discriminator the demodulated noise has the well 
known triangular spectrum and, as illustrated previously, the signaling 
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speed capability is essentially the same as an AM system of equal band 
width. 

To compare experimentally these two general types of discriminators a 
740-cycle band system with linear discriminator and 350-cycle shift was 
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Fig. 17.—Characteristics of 1920-cycle bandpass filter and associated discriminator 
consisting of two bandpass filters, 


FILTER LOSS FROM MIDBAND 
IN DECIBELS 


\ 
\ 
\ 
\ 
\ 
\ 
‘ 
\ 
\ 
\ 
‘\ 
A 
\ 


DISCRIMINATOR y 
OUTPUT 





oOo 
DISCRIMINATOR OUTPUT IN VOLTS 


1.2. 43 1.4 l5 1.6 1.7 186 49 
FREQUENCY IN KILOCYCLES PER SECOND 


Fig. 18.—Characteristics of 295-cycle bandpass filter and associated discriminator 
consisting of two bandpass filters. 


compared with a system with a bandwidth of about 1900 cycles, a discrimi- 
nator consisting of two 740-cycle bandpass filters, and an 850-cycle shift. 
The results are shown in Fig. 19. The two systems are seen to reach failure 
distortion values at the same signal-to-noise point, with. the linear dis- 
criminator becoming about 3 db superior at distortions around 5%. A 
second comparison was made using roughly equal bandwidths. A 295-cycle 
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band with a discriminator consisting of two bandpass filters and 170-cycle 
shift was compared with a 230-cycle band with a linear discriminator and a 
shift of 140 cycles. The results are shown in Fig. 20. The linear dis- 
criminator in this case appears to fail slightly sooner but shows a superiority 
of about 2 db in the 5% distortion region. Due to the rounded character- 
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Fig. 19.—Linear discriminator versus two-bandpass-filter discriminator having the 
same signaling speed capability. Effect of thermal noise on peak distortion—80- cycle 
cutoff low-pass filter. 
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Fig. 20.—Linear discriminator versus two-bandpass-filter discriminator with equal 
bandwidths. Effect of thermal noise on peak distortion—80-cycle low-pass filter. 


istic of the two bandpass filters used as a discriminator in the second com- 
parison (Fig. 18) the difference in discriminators is less than in the previous 
test. It has been found’ that for @ given signaling speed capability almost 
twice the bandwidth is required if a two-bandpass-filter discriminator is used 
instead of the linear type. This added band width does not appear to cause 
any loss in signal-to-noise capabilities as to the failure point. The less sharp 
breaking point, however, makes the linear discriminator superior for moder- 
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ate and low distortions. For a system occupying @ given bandwidth the 
- two-bandpass-filter discriminator provides some improvement at the failure 
point but is still somewhat inferior to the linear discriminator at moderate 
distortions. More importantly the two-bandpass-filter discriminator im- 
pairs the signaling speed capabilities to an extent which depends upon the 
shape of the cutoff of the filters used. 


Bandwidth After Demodulation (Low-Pass Filtering) . 


In an Am system the low-pass filtering after demodulation can, to a large 
degree, make up for a greater than necessary bandwidth before demodula- 
tion. During marking intervals the added noise admitted by a wide band 
causes noise in the demodulated output at frequencies higher than the signal- 
ing frequency and this can be filtered out, unless the noise is so great as to 
over-modulate the carrier. During spacing intervals there is no carrier and 
hence only the noise is rectified. Added noise admitted by a wide band 
causes not only higher-frequency components in this rectified noise, which 
may be filtered out, but also an increase in the d-c.component. This tends 
to cause marking bias of the received signals as the noise level increases. 

In an FS system, where the carrier is present continuously, the added 
noise from a wider band produces high-frequency noise components in the 
demodulated output which can be filtered out by the low-pass filter if the 
noise level is low. As the noise level increases there are short intervals 
when the noise envelope exceeds the carrier. . The action of the limiter is 
to give preference to the greater signal, in this case the noise, and since the 
noise will appear to the discriminator as a carrier fluctuating around mid- 
band as a center, the demodulated output momentarily dips toward zero. 
As the noise increases, the duration and frequency of these holes in the signal 
increase. The low-pass filter, by excluding frequencies considerably in 
excess of the maximum signaling speed, prevents these holes in the signal 
from producing false or extra transitions in the telegraph signal output. 
The low-pass filtering, however, cannot prevent the true transitions from 
being displaced by this type of noise component since the signal is obliterated 
momentarily. Its most important function is to prevent a breakup in the 
signal output until a fairly high distortion is reached. For noise peaks 
exceeding the carrier the low-pass filter of an AM system also serves much 
the same purposes. 

In Fig. 21 is shown the effect of changing the bandwidth of the low-pass 
filter in an FS and in an AM system in the presence of thermal noise. The 
effect of a narrower low-pass filter is seen to consist mainly in shifting the 
breaking point toward a higher noise level. Similar characteristics for the 
case of impulse noise are shown in Fig. 22. 
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Magnitude of Frequency Shift in Relation to Bandwidth 


A frequency-shift transient in a band of given width has a wave shape 
much like that of an amplitude transient in the same band provided the shift 
is symmetrical and not over 50% of the bandwidth.® If the frequency shift 
approaches the total width of the band the transient is of such shape as to 
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Fig. 21.—Effect of the low-pass filter cutoff frequency on distortion in the presence of 
thermal noise—740-cycle band. 
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Fig. 22.—Effect of the low-pass filter cutoff frequency on distortion in the presence of 
impulse noise—740-cycle band. 


cause distortion and to make the system more susceptible to noise. A very 
small shift results in a low amplitude of demodulated signal, which is more 
readily distorted by noise and biased by frequency drifts. It is of interest, 
however, that the signal-to-noise ratio of the demodulated signal is not 
proportional to frequency shift for high noise conditions. As described 
before, noise peaks tend to reduce momentarily to zero the output from a 
balanced discriminator. The amplitudes of the dips or holes are thus about 
one-half the demodulated signal amplitude for any value of shift. This 
tends to maintain a constant signal-to-noise condition and this characteristic 
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is illustrated by Fig. 16:in which the breaking point with a 100-cycle shift 
occurs only 2 db before that with a 400-cycle shift although the difference in 
actual signal amplitude is 12 db. Figure 23 shows the effect on signal-to- 
noise ratio of progressively varying the frequency shift while the bandwidth 
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Fig. 23.—Effect of magnitude of frequency shift on distortion produced by thermal 
noise—80-cycle cutoff low-pass filter. 
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Fig. 24.—-Effect of frequency drift on distortion produced by thermal noise in FS 
transmission—740-cycle band and 80-cycle cutoff low-pass filter. 
is kept constant. There is a fairly broad region in which the signal-to-noise 
ratio is little affected by the amount of frequency shift. For optimum 
results a frequency shift of 50% to 60% of the band width is indicated, and 
thus has been used in most of the experimental results given herein. — 


Frequency Instabilities 


Frequency drift of the carrier input to the receiving terminal in general 
causes biased signals and if severe enough results in failure of the system. 
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In Fig. 25 is shown the effect of frequency drift on signal bias in an AM 
system. In an AM system little bias is produced until the carrier reaches 
the cutoff region of the filter. The bias then becomes rapidly negative due 
to the increased loss and decreased amplitude of demodulated signal. When 
an automatic gain control arrangement is used the bias becomes positive 
due to a distorted envelope shape. The demodulated wave form deter- 
mines the degree of sensitivity to frequency drift and depends on the band- 
width both before and after demodulation. 

In an /*'S System using a linear discriminator, frequency drift changes 
the d-c. component of the signals and thus changes the operating point on 
the demodulated wave. ‘The amount of bias depends upon the slope of the 
wave front and is thus affected by the amount of low-pass filtering. The 
effect of frequency drift on bias for a number of I'S systems is shown in Figs. 
26 and 27. Ifa two-bandpass filter type of discriminator is used the system 


BIAS IN PER CENT 


w 
re) 





-500 -400 -300 -200 -100 O 100 200 300 400 
FREQUENCY DRIFT IN CYCLES PER SECOND 


Fig. 25.—Signal bias versus frequency drift for AM transmission in a 740-cycle band. 


is insensitive to moderate frequency drifts due to the flat pass bands as illus- 
trated in Fig. 26. The relative shape and amplitude of the signal from a 
linear discriminator does not change appreciably with frequency drift; 
only a d-c. displacement occurs. This makes it desirable to have the low- 
pass filter coupled to the output amplifier: by a network which passes only 
the useful signaling frequencies and blocks the d-c. and very slow drift com- 
ponents. When this is done the effect of frequency drift on bias is not 
greatly different from that for an AM system, as may be seen by comparing 
the dotted curves on Figs. 26 and 27 with Fig. 25. 

The general method of performing this d-c. elimination is illustrated in the 
block diagram of Fig. 29. The output from the low-pass filter is passed 
- through a coupling network which blocks the d-c. and passes the useful 
signaling frequencies. The output of the coupling network is passed 
through a positive feedback nonlinear amplifier which has but two output 
conditions representing the mark and space of the telegraph signal. The 
feedback network passes d-c. and low frequencies so as to just compensate 
for the loss of the coupling network. ‘The time constant of the coupling 
network may be made large enough so that the signal wave form into the 
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d-c. amplifier is practically the same as at the output of the low-pass filter. 
The operating point on the demodulated wave may be readily adjusted by 
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‘ig, 26.—Signal bias versus frequency drift for FS transmission—wide filter bands. 
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Fig. 27.—Signal bias versus frequency drift for FS transmission—narrow filter bands. 


adjusting the bias voltage at the input to the d-c. amplifier. This arrange- 
ment differs from the impulse type of d-c. elimination, in which the demodu- 
lated wave form is effectively differentiated to form pulses but a fraction of a 
unit dot in length. ; 
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When d-c. elimination is used, operation with FS signals decentered in the 
pass band causes little bias but it does cause a loss in signal-to-noise ratio, 
especially at high noise levels. Due to the effect of noise peaks in causing a 
dip toward zero in the demodulator output, the effect of noise becomes 
exaggerated during the signal condition which is farther from midband. 
This change in signal-to-noise condition with frequency drift is shown in 
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Fig. 28.—Signal bias versus carrier level variation for AM transmission. 
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Fig. 29.—Block diagram of dc. elimination and restoration method used to minimize 
signal bias caused by frequency drift. 
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Figs. 24 and 30. .A comparison between a two-bandpass-filter discriminator 
and a linear discriminator with d-c. elimination as regards frequency drifts 
in the presence of noise is shown in Fig. 31. 

Radio telegraph systems operating in the H.F. region require a high degree 
of frequency stability if narrow bandwidths are to be used. Because of 
the several frequency conversions involved a number of different oscillators 
or frequency sources are involved but usually the major burden of frequency 
stability rests on the transmitter exciter and the high-frequency beating 
oscillator of the receiver. 

Various methods of automatic frequency control may be used to hold the 
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carrier input to the demodulator at the correct frequency. In the case of 
FS signals the control may be arranged to operate only on the marking fre- 
quency or to utilize both the mark and space conditions. It is preferable 
to have inherent frequency stability rather than to compensate for the drift 
at the receiving end since it is difficult if not impossible to provide an auto- 
matic frequency control which will not reduce the transmission capabilities — 
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Fig. 30.—Effect of frequency drift on distortion produced by impulse noise in FS 
_ transmission—740-cycle band, 100 cycle frequency shift, 80-cycle cutoff low-pass filter. 
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Fig. 31.—Linear discriminator versus two-bandpass-filter discriminator when frequency 
drift occurs in the presence of thermal noise——4 db rms noise-to-carrier ratio. 


of the system in the presence of noise and other interference. Best results 
are obtained if the frequency stability is high enough to require but a very 
slow correction, which usually dictates some mechanical rather than elec- 
tronic tuning arrangement. Manual retuning may be found satisfactory 
where stability is reasonably good provided care is taken in making the ad- 
justments. Suitable frequency stability with the retention of flexibility in 
frequency adjustment may bé obtained by frequency sources making use of 
a combination of crystal oscillators and high-stability variable oscillators of 
lower frequency. 
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Level Variations 


Extreme and rapid variations of received level exist in H.F. radio trans- 
mission. It is upon the ability to accommodate these level variations that 
the merits of an H.F. radio telegraph system must largely be judged. IFS 
telegraph shows its outstanding advantage in this respect. 

In an AM system in order to obtain zero-bias signals and optimum signal- 
to-noise conditions the operating point must be near the half amplitude point 
on the demodulated wave. This means that complete failure will result for 
a drop in level of 6 db unless some compensating arrangement is provided. 
The slope of the bias-versus-level characteristic depends upon the slope of 
the demodulated wave which in turn depends on the bandwidth of the sys- 
tem and upon the degree of low-pass filtering. The bias-versus-level char- 
acteristics of some AM systems are shown in Fig. 28. Where the level 
variations are relatively slow compared to the signaling speed, autoniatic 
gain control circuits can be used to maintain a nearly constant level into the 
demodulator. However, where large rapid level changes occur, as in the 
H.F. range, it is seen that a narrow band AM system would fail completely 
regardless of the amount of transmitted power. For printer operation 
over an AM system in the H.F. range a fairly wide band and little low-pass 
filtering should be used, so as to keep the wave shape of the signals as 
square as possible and thus obtain a fairly flat bias-versus-level charac- 
teristic. By adjusting the operating point low on the demodulated wave, 
approaching the spacing noise level, the greatest possible range of acceptable 
rapid level change will be obtained. The slower level change components 
may be handled by the usual automatic gain-control circuits. This will 
cause the bias of the signals to average somewhat marking but the peak 
distortions will be kept to a minimum. 

In an F'S system no bias is produced so long as both the marking and spac- 
ing frequencies are affected alike, with their received levels remaining equal. 
Such non-selective fading conditions cause no distortion even when they 
occur at quite rapid rates. If a balanced type of discriminator is used, 
amplitude limiting is not essential to obtaining this immunity from non- 
selective variations in attenuation. It is only when the mark and space 
levels are different that bias results. In Fig. 32 are shown bias versus mark- 
to-space level ratio characteristics both with and without a limiter. More 
bias exists when there is no limiter because the amplitude of the demodulated 
wave is directly affected and consequently the low-pass filtering also becomes 
a factor. With a limiter the amplitude of the demodulated wave is held 
constant and the amount of low-pass filtering has no effect on bias. Some 
bias is still produced, however, due to the differently shaped frequency 
transients in the passband of the receiving system when a level change occurs 
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Fig. 32.—Signal bias versus mark-to-space level ratio in FS transmission—740-cycle 
band, 350-cycle frequency shift. 
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Fig. 33.—Signal bias versus mark-to-space level ratio in FS transmission—-80- cycle 
cutoff low-pass filter. 


at the moment when the frequency changes. In Fig. 33 this bias effect is 
demonstrated for various bandwidths. Tor moderate mark-to-space level 
ratios the bias effect is small and linear, with a slope which usually varies 
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inversely with a change in bandwidth. At extreme level differentials the 
bias may rise very rapidly due to the amplitude and phase characteristics 
of the input passband; transients from the greater amplitude condition may 
severely interfere with the lower amplitude condition. For the types of 
bandpass filters used in the tests it appeared that the amount of level dif- 
ference required to produce 20% bias did not change greatly with band- 
width. Severe wave shaping of the signal at the transmitter was found to 
be an aid in reducing the bias effect due to such transients but the charac- 
teristic distortion became too great to give any practical improvement. 

The fading modulator used in obtaining the data for Fig. 33 caused no 
change in phase. Selective fading over an actual radio circuit would involve 
considerable phase shift and greater distortion might be expected. The 
data of Fig. 32 were obtained by use of the phase control associated with the 
crystal filter of the radio receiver to vary the loss-versus-frequency charac- 
teristics of the receiving pass band and thus cause unequal mark and space 
amplitudes. This method gave an amplitude and phase characteristic for 
the transmission band more like that over an actual radio circuit. 


MULTIPATH PROPAGATION EFFECTS 


The rapid fading conditions prevailing in the H.F. range are brought 
about by multipath propagation. Under such conditions, the signal in- 
duced in a receiving antenna by a distant transmitter may be the resultant 
of two or three separate waves each propagated over a different path. If 
two waves arrive over paths differing in length by an odd number of half 
wavelengths the resulting 180° phase difference causes maximum cancella- 
tion. On the other hand if the paths differ in length by an integral mul- 
tiple of whole wavelengths the waves arrive in-phase and maximum rein- 
forcement results. The difference in path lengths may at times be as great 
as 500 to 1500 kilometers (delay times of 2 to 5 milliseconds) which in the 
H.F. region corresponds to thousands of wavelengths. Under these maxi- 
mum conditions waves at one frequency may arrive in phase while waves at a 
frequency a few hundred cycles away may arrive in phase opposition. Since 
the path lengths are constantly changing, the transmission at a given fre- 
quency is subject.to wide variations in amplitude and phase with time. 
When the difference in path lengths is not great enough to cause frequencies 
in one portion of a communication channel to fade differently from those in 
another portion the term ‘‘non-selective” or “flat” fading is applied. When 
the difference in path lengths becomes great enough to cause considerable 
amplitude or phase distortion over the transmission band the term “selec- 
tive’ fading is used. Since the propagation paths existing at a given 
moment vary for different antenna sites, the fading patterns obtained from 
two or three antennas separated by several wavelengths usually show a 
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considerable phase difference so that a given frequency is not likely to fade 
into the noise level at all antennas simultaneously. By employing separate 
receivers for each antenna and suitably combining or selecting the demodu- 
lated outputs, a system is obtained which is much less susceptible to fading. 
Such a method is called space diversity reception. Inasmuch as fading over 
a given combination of paths is highly selective with respect to frequency 
much the same effect is obtained by frequency diversity reception. When 
this method is employed the intelligence is transmitted on two or more fre- 
quencies simultaneously and then received by separate receivers from a 
single antenna and the resulting demodulated signals combined or selected 
as for space diversity. 

In te'egraph transmission large differences in delay over two separate 
propagation paths cause the telegraph signal transitions to arrive at different 
instants over the two paths. Thus, there are intervals of overlap when a 
marking condition is received over one path and a spacing condition over a 
second path. When two components of nearly equal amplitude arrive at 
nearly 180° phase difference a signal transition may involve large and sudden 
amplitude and phase changes. The resulting transients in the bandpass 
networks of the receiving equipment may cause fortuitous distortions con- 
siderably greater than the difference in delay times over the two paths. 
The wider the pass band of the receiving system the shorter the duration of 
these fortuitous transients and hence the less the distortion. This phe- 
nomenon is one of the determining factors in the selection of bandwidth and 
frequency shift to be used in a given application of FS telegraphy. It 
becomes of increasing importance when the circuits are long and at higher 
signaling speeds such as are used in time-division multiplex methods. 

In an AM system the effect of large differences in path lengths is usually a 
filling in of the spacing intervals with resulting marking bias. In an FS 
system the overlap time and associated transients may add to either mark- 
ing or spacing intervals in a random fashion depending on the amplitude and 
phase conditions at each transition. ‘The overlapping of the mark and space 
frequencies in FS transmission can sometimes be heard in an AM receiver as 
short pips of audio tone at each transition, the audio tone being the beat 
between the two frequencies. 


Use of Superimposed Phase Modulation 

Superimposed phase modulation has sometimes been employed as a 
simple means for achieving a certain amount of frequency diversity both in 
AM and FS telegraph systems. ‘This consists in causing the radiated signal 
to oscillate continuously through a small phase angle at a rate relatively 
high compared to the dotting speed. Phase modulation spreads the energy 
of the signal over a wider frequency band so that the complete loss of the 
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signal through selective fading becomes less probable. The spectra gen- 
erated by sinusoidal phase modulation of 1.0, 1.4, and 2.0 radians are shown 
in Fig. 34. Most of the energy is seen to be concentrated in the carrier and 
first order sidebands. Less than 1.0 radian of modulation results in too 
. little amplitude of the sidebands, while more than 1.5 radians results in too 
wide spread of energy outside the first order sidebands. The center three 
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Fig. 34.—Frequency spectra for sinusoidal phase modulation. 


components are equal at about 1.4 radians. In the case of FS the phase 
modulation frequency appears as a variation in amplitude of the signals from 
the discriminator. For AM no additional amplitude variation is caused by 
the phase modulation if there is no selective attenuation in the medium, 
but if such exists the phase modulation frequency or a multiple thereof 
appears in the rectified signal. To permit these unwanted amplitude varia- 
tions to be removed by the low-pass filter so as not to break up the signals, 
the phase modulating frequency should preferably be ten or more times the 
maximum signaling frequency. 
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A test of superimposed phase modulation on FS signals was made over a 
radio circuit approximately 200 miles in length. A frequency shift of 850 
cycles, and one radian of 200-cycle phase modulation, was used. A 60- 
word-per minute test sentence was transmitted and received without space 
diversity. It was found over a period of several hours that the phase modu- 
lation, on the average, gave a decrease in printed errors of about 50% when 
the error rate was in the proximity of 1 to 2%. For short intervals the re- 
duction in errors was often considerably greater. The use of 200-cycle 
phase modulation when space diversity is used provides little or no improve- 
ment and is therefore undesirable. 

A more effective way of employing phase modulation with FS signals 
would be to use a phase swing of +1.4 radians at a frequency of 2 to 3 times 
the frequency shift and to demodulate separately the three major com- 
ponents of the signal, thus obtaining in effect a triple-frequency diversity 
system. This of course involves quite a wide transmitted band, but it 
might be of use in cases where space diversity is impossible, such as on 
board ships. When a space diversity arrangement is feasible it is much to 
be preferred. 


Diversity Operation 


To obtain reliable operation in the H.F’. range it is common practice to 
employ space diversity reception. The use of frequency diversity, with the 
increase of transmitted power and greater frequency space required, is sel- 
dom justified if space diversity reception can be arranged. For AM radio 
telegraph, double or triple-space diversity receiving arrangements are 
frequently used. Since an FS signal generally covers more frequency 
space, it is even more likely to be mutilated by selective fading than an 
AM signal. It has been found, however, that a double-space diversity sys- 
tem for FS signals usually gives sufficient diversity action provided it is of a 
type that permits switching between channels at signaling speed without 
causing appreciable distortion. This is necessary since it is a frequent 
occurrence that the mark of one channel may fade, leaving a good space, 
while the opposite may occur on the second channel. Since an FS system 
‘can accept rapid level changes, the main purpose of diversity methods is to 
insure that both the mark and space portions of the signal will be received 
above the noise level. In the case of AM telegraph, since it cannot accept 
rapid level changes, diversity operation is important not only in keeping the 
signal above the noise but also in averaging out some of the rapid level 
changes. For this reason AM systems usually show considerable improve- 
ment in going from double to triple diversity. It would be expected that a 
like change would show much less improvement in an FS system. 
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Diversity Channel Selection 


The method employed to combine or select the channels of a diversity 
system is of great importance. For an AM system the relatively simple 
method of using a common load circuit for the diode detectors of the diver- 
sity channels is generally used. By deriving a common AVC voltage from 
the combined output and by properly adjusting the receiver sensitivities a 
fairly constant output is obtained. The parallel connection of the diode 
detectors causes the stronger signal to effectively block the weaker signal 
thus giving a fairly sharp diversity selection characteristic. The problem 
of combining the diversity channels of an FS system is more complicated 
mainly because of the amplitude limiting. If amplitude limiting is used in 
each diversity channel before demodulation, the resulting constant ampli- 
tude signals convey no information as to their relative amplitudes as re- 
ceived from the antennas. Any diversity selection must then be obtained 
by some indirect method. It is necessary to furnish some selecting device 
since the noise from a faded channel, if added directly to a good signal from 
another channel, will cause high distortion. 

In an early frequency shift system employing a two-bandpass filter dis- 
criminator (shown previously in Fig. 17) it was found that for a poor sig- 
nal-to-noise condition the sum of the outputs of the mark and space rec- 

tifiers increased above that for a good signal-to-noise condition. ‘This 
increase was utilized to suppress the output of the poorer channel and em- 
phasize that of the better channel. Although neither the degree of diversity 
selection nor the speed of response was as good as might be desired, fairly 
satisfactory results were obtained. 

Another method which has been used involves the derivation of control 
currents or voltages proportional to the amplitudes of the incoming signals 
which in turn select the better diversity channel by some type of gate action. 
The time constants of the control circuits must be low enough to permit 
switching at signaling speed without introducing considerable distortion. 
The gate circuits must also be of a type which does not introduce interfering 
transients or otherwise allow the control voltages or currents to interfere 
with the signal. This method permits very sharp diversity selection and 
has the capability of approximating ideal results although it becomes some- 
what involved in a practical form. 

A considerably simpler method has been used in some of the more recent 
FS terminals. It is based on the use of a single-amplitude-limiter through 
which pass the signals of both diversity channels. This is made possible 
by arranging the two signals at the input to the limiter to be at different 

frequencies. At the output of the limiter the two signals are separately 
demodulated and then combined. When one of the signals is considerably 
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greater in amplitude than the other at the input to the limiter the relative 
difference in level is increased by an additional amount of about 6 db at the 
limiter output. The limiter output may be considered as the stronger sig- 
nal frequency-modulated by the weaker signal. For small modulation 
indices the amplitude of the first order sideband is approximately one-half 
the modulation index thus explaining the 6 db added difference in level at 
the limiter output. As the input levels approach equality the added level 
difference decreases to zero. A block diagram indicating the arrangement of 
such 4 diversity system is shown in Fig. 35. Tests were made of both paral- 
lel and series connections of the two discriminator outputs. With a parallel 
connection the discriminator having the greater output blocks the rectifier 
output of the other discriminator and thus gives a sharp diversity selection 
characteristic. However, the level ratio of the channels at which a switch 
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Fig. 35.—Block diagram of dual-diversity FS receiving system using a common limiter. 


takes place is affected not only by the ratio at the limiter input but also by 
frequency drift and discriminator slope. With the series connection only 
the input level ratio at the limiter input affects the diversity switching; 
this arrangement was therefore selected as the preferred method although its 
selection characteristic is not sharp. The series and parallel combining 
characteristics are shown in Figs. 36 and 37. 

Various tests were made on a terminal having a 1500-cycle bandwidth 
and using the series combining method to determine the signal-to-noise 
characteristics under different conditions of diversity fading. A frequency 
shift of 850 cycles was used and the midband frequencies of the two diver- 
sity channels at the common limiter input were 30 and 35 Kc. Figure 38 
shows the distortion versus signal-to-noise ratio characteristics of each 
channel separately and in diversity combination for various relative level 
conditions of the two channels. During diversity operation equal noise 
levels were maintained in the two channels and various combinations of 
level differences of the two channels were preserved as the whole signal 
level combination was varied. The level differentials are indicated in 
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the figure foreach curve. The signal-to-noise level scales refer to the highest 
level portion of the diversity signal. Since the amplitude modulator which 
was used to simulate the selective fading did not produce ‘phase shifts or 
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Fig. 36.—Diversity combination characteristic obtained by series addition of discrimi- 
nator outputs—levels measured at output of 400 kc I.F. amplifier. 
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Fig. 37.—Diversity combination characteristic obtained by parallel addition of dis- 
criminator outputs—levels measured at output of 400 kc I.F. amplifier. 


transmission delays as would actually exist over H.F. radio circuits the actual 
distortions shown by the curves are optimistic. 

The ideal diversity selection circuit should theoretically give a signal- 
to-noise characteristic identical to that of a single channel under the signal- 
to-noise condition corresponding to the signal of the best momentary re- 
ception. It will be seen that the test results of Fig. 38 approach this limit 
within 2 or 3 db at a peak distortion of 20%. Part of this difference is 
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due to the dissimilar bandpass characteristics of the two channels of the 
experimental unit used for the tests and part due to the lack of an extremely 
sharp diversity selection. It should be pointed out that the conditions under 
which the theoretical maximum diversity signal-to-noise condition may be 
reached are very hard to obtain in practice. If the noise levels in the two 
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Fig. 38.—Peak distortion versus signal-to-noise ratio characteristics for dual-diversity 
operation with a common limiter. Signal-to-noise measured at output of 400 kc I.F. 
amplifier. 


channels are not equal, and if the diversity selecting method does not at all 
times exclusively select the channel with the greater signal, the distortion 
characteristic will deteriorate accordingly. Because the AVC sensitivities. 
of two radio receivers may differ considerably the noise levels cannot be 
maintained closely the same and usually no provision is made for determin- 
ing the noise level except by ear. The slightly better diversity action which 
can theoretically be obtained is therefore felt to be of doubtful usefulness 
under actual operating conditions. The common limiter method has the 
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advantage of being simple in that the diversity action is obtained in the 
transmission circuits directly and no added switching circuits or adjustments 
are required. Tests of this type of diversity selection in the field have indi- 
cated a marked superiority over earlier FS terminal equipment. 

In connection with the use of a common limiter care must be taken in the 
selection of the two-channel frequencies. Frequencies having nearly 
integral ratios such as 3:5 and 5:7 produce disturbing amplitude modula- 
tion of the demodulated signal. The frequencies should be chosen so as to 
avoid low integral ratios; then all amplitude modulations are negligible or 
easily filtered out. Where frequency drift is to be allowed for, the fre- 
quencies should be chosen so as not to approach a low integral ratio at any 
place in the expected drift range. 

If the radio receivers associated with a space diversity FS system have 
automatic gain control it must be a common control so the receivers will 
change gain equally. The use of common AVC prevents overloading of 
the receivers as the received signal strength varies. If no common AVC 
is available the receivers should be operated in the manual gain-control 
condition. 7 


CONCLUSIONS 


General Comparison of FS and AM Carrier Telegraphy 


The foregoing sections have compared the characteristics of AM and FS 
carrier telegraph transmission under various conditions. Whether or not 
FS would prove to be the preferable method for a specific communication 
use depends largely on the transmission medium and the quality of trans- 
mission desired. As regards frequency space requirements, both methods 
provide essentially the same signaling speed capability for a given band- 
width. . 

As to the ability to transmit through noise, FS has an advantage of 3 to 
4 db at distortions approaching the failure point when equal bandwidths 
are compared. At lower distortions the advantage of FS is 6 db or more 
so that it is attractive in this respect for tandem operation of several tele- 
graph sections where regeneration of signals is not practiced. When fre- 
quency space permits wider bands, with correspondingly increased fre- 
quency shifts, the signal-to-noise advantage of FS over AM increases for 
low noise levels. Wide band FS therefore provides a means of obtaining 
higher quality circuits if the noise level is not too great. 

The AM method is basically less susceptible to frequency variations 
than is the FS method. However, as has been illustrated, frequency drift 
can be compensated for by d-c. elimination so as to make FS comparable 
to AM in this respect. 
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FS transmission is essentially immune from effects of non-selective level 
variation, even when extremely rapid, and in this characteristic displays 
its most outstanding advantage over AM. 


Operation Over Wire Circuits 


A wire circuit usually provides a transmission medium having a low noise 
level with slow and relatively small variations in attenuation. Such cir- 
cuits, when equipped with suitable automatic gain control, allow stable 
operation with AM telegraphy and but little improvement could probably 
be obtained by using FS. The choice between AM and FS under such rela- 
tively ideal conditions becomes one of economic considerations of the ter- 
minal equipment and carrier supply. However, when FS is applied to 
multichannel systems the problem of interchannel interference requires 
attention. For wire circuits having high noise levels or sudden changes in 
attenuation the use of I'S instead of AM provides considerable improvement 
and in severe cases the FS method may be a necessity for satisfactory opera- 
tion. Wide band FS operation with its sharper breaking distortion-versus- 
noise-level characteristic gives a low value of rms-to-peak distortion which 
would be especially advantageous for tandem operation: However, the 
necessary frequency space for wide-band operation is not usually economi- 
cally justified for wire line operation. 


Operation Over Radio Circuits 


For operation over radio circuits providing stable conditions similar to 
those on wire circuits the FS method does not show a great advantage over 
the AM method. In the case of long distance telegraphy in the H.F. range, 
however, I'S shows a marked advantage over AM. This is because of the 
rapid fading and high noise conditions which commonly prevail in the 
H.F. region. The amount of rapid variation in marking level that an AM 
system can accommodate is less than the difference between marking and 
spacing levels that an FS system can tolerate. In the worst case of selective 
fading the level differences between the mark and space frequencies might 
approach values equal to the short time level swings of a single frequency, 
but in general would be less. A given condition of selective fading thus 
causes less distortion in an FS system than in an AM system. FS allows 
the use of narrow bands without much loss in signal quality in the presence — 
of fading, whereas AM does not. FS therefore is essential for satisfactory 
operation of closely spaced narrow band H.F. radio channels. Where fre- 
quency space is not restricted and wider bands are used to permit consider- 
able frequency drift, the improvement afforded by FS over AM is materially 
less. To obtain optimum results from an AM system, however, requires 
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more careful adjustment and more attention than does an FS system. This 
is partly due to the amplitude limiter in the FS system which results in a 
constant amplitude of signal from the discriminator and partly due to the 
fact that an FS signal is no more subject to noise interference during the 
spacing condition than during the marking condition. ‘Therefore the operat- 
ing point on the demodulated wave may be set and left for long intervals 
even though transmission conditions vary widely. This greater ease in 
maintaining good adjustment of the equipment probably accounts for some 
of the apparent improvement in changing from an AM to an FS system. 

It should be noted that a system may fail either because of level variations 
well above the noise level or because of the signal becoming submerged 
in noise. Ifa system fails because it can accept only moderate level varia- 
tions, an increase in transmitted power will provide no improvement since 
the level variations will remain the same as before. On the other hand, a 
system which can accept very wide variations in level will show improve- 
ment upon increasing transmitted power up to the point where no failures 
occur due to an unfavorable signal-to-noise ratio. © 

The over-all improvement obtained in changing from AM to FS radio 
telegraph is sometimes expressed as a ratio of transmitted powers required 
to give equivalent transmission results over the two systems. Such a ratio 
fluctuates widely depending upon the prevailing conditions. With little 
fading the improvement ratio will be mainly due to the better signal-to- 
noise obtained with FS and may be less than 5 db. Under severe fading 
conditions no amount of power may give good results with AM while FS 
may be satisfactory. Thus the power ratio would become infinite. By 
making a long-time comparison an average power ratio figure may be 
- found which gives equal average error rates in the printed copy from each 
system. Such tests’ between a triple space diversity AM system and a 
double space diversity FS system have indicated a power ratio of 11 db 
in favor of the latter when the error rate was 0.1 to 0.5 per cent. 

When the two systems are thus made equal by adjustment of transmitted 
power, more errors due to the signal becoming submerged in noise occur in 
the FS system to compensate for a larger number of errors in the AM system 
due to rapid level changes. Often the reason for changing a radio telegraph 
system from AM to FS is to increase the reliability of the circuit and not 
just to save transmitted power. To insure a definite improvement in such 
cases the carrier level should not be decreased more than about 6 db. 
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Reflections from Circular Bends in Rectangular 
Wave Guides—Matrix Theory 


By S. O. RICE 


A method of computing reflections produced by circular bends in rectangular 
wave guides is presented. The procedure employs the theory of matrices. Al- 
though the matrix equations are quite simple, a considerable amount of calculation 

. is necessary before quantitative results may be obtained. Fortunately, the ap- 
proximate formulas pertaining to gentle bends hold surprisingly well for rather 
sharp bends. These formulas are obtained by a limiting process from the matrix 
equations. The approximate formula for reflection from an H-bend (in which the 
magnetic vector lies in the plane of the bend) generalizes an earlier result due to 
R. E. Marshak. The corresponding formula for the E-bend appears to be new. 


INTRODUCTION 


NUMBER of investigators have studied the propagation of electro- 
magnetic waves in a bent pipe of rectangular cross-section, the bend 
being along an arc of a circle. H. Buchholz’, S. Morimoto’, and W. J. 
Albersheim’ have employed Bessel functions to express the field in the bend. 
The form assumed by the field when the radius of curvature of the bend 
becomes large has been obtained by K. Riess" and R. E. Marshak’ who use 
approximations suited to this case. Marshak also obtains expressions for 
‘ various reflection and transmission coefficients. A discussion of the subject 
using rather simple but approximate analysis is given on pages 324-330 of 
a text book’ by S. A. Schelkunoff. The Bessel function approach is also 
sketched in the same section. 

Here we study the disturbance produced when a wave goes around a 
circular bend (of some given angle) in a rectangular wave guide, the guide 
being straight on either side of the bend. Especial attention is paid to the 
dominant mode reflection coefficients gip and doi corresponding to H-bends 
and E-bends, respectively. As equations (4.2-6) and (4.44) show, these. 
reflection coefficients (which are of the nature of voltage rather than power 
reflection coefficients) vary inversely as the square of the radius of curva- 
ture of the bend when the bend is gentle. The substance of (4.2-6) has 
been given by Marshak’ for the important case in which only the dominant 
mode is propagated and the angle of the bend not too small. 

When the bend is so sharp that the formulas mentioned above do not 
apply the reflection coefficients may be computed from the rather simple 
looking matrix expressions (2.3-3) together with (2.3-4). However, their 
appearance is deceptive and, as is shown by the numerical work in Part V, 
considerable labor is necessary to obtain an answer. 
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The gentle bend formulas were obtained from the matrix equations by the 
limiting process described in Part III. It seems likely that the matrix 
method, which is similar to the method used in an earlier paper’ on trans- 
mission line equations, may be applied to other wave guide problems. 
With this thought in mind, the development of Parts II and III has been 
couched in general terms. . 

The matrices used in the present theory are of infinite order since the 
guide may support an infinite number of modes of propagation. This fact 
makes it difficult to justify all the steps in our analysis, and we do not at- 
tempt to do so.* Despite this lack of rigor, I believe that the procedures 
given here lead to the correct results since they yield, for gentle bends, 
expressions obtained by Buchholz and Marshak. Moreover, although 
numerical results tabulated in Part V were obtained by using matrices of 
only the second and third order, they indicate a rapid convergence as the 
matrix order is increased. 


PART I 


PROPAGATION OF WAVES IN GUIDE 
1.1 Propagation in a Straight Wave Guide 


Rather general expressions for the electric and magnetic intensities H and 
H ina field are (see pp. 127-128 of Refezence’) 


E 


II 


—iwpA + _ grad div A — curl B 

WE 

, (1.1-1) 
H 


I 


curl 4 + a grad div B — iweB 
Tops 


The field is assumed to vary with the time ¢ as e“’, w is the radian fre- 
see pw the Derma Duity and ¢ the dielectric constant (for free space 

= 1.257 x 10° henries/meter, e = 8.854 x 10°” farads/meter). The 
vector potentials A and B satisfy the wave equations 


VA=cA, VB=0B 
V’ = Laplacian operator (h1-2) 
2 2 
og = Wie 
In dealing with bends, it is convenient to choose A and B normal to the 
plane of the bend. In our notation, this plane is always taken to be the x, 
z plane so that A and & are parallel to the y axis. The ¢ axis is parallel 


* Similar questions arise in the rigorous treatment of an infinite set of linear 
equations. A discussion of this subject is given in Chap. ITI of Reference’. 
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to the guide axis and, for the straight guide of the section, the guide walls 
are sections of the planes + = 0,4 = a,y=0,y = 0. 
Thus, a general wave traveling in the positive z direction may be. de- 
scribed by the two functions (which represent the magnitudes of A and B) 
eye sin (xmx/a)cos (riy/b) 
ne | (1.1-3) 
Paeie ee na =0,1,2,--- 
25 dun €™" cos (rmax/a)sin (riry/b) 
(1.1-4) 
m= 0,1,2,-°:; m= 1,2,3,--- 
where the coefficients re and di, are constants and the plus signs indicate 
propagation in the positive ¢ direction. 
The propagation constant I,, is obtained from 


Tan = ot + (xm/a)? + (rn/b)*, o = i2m/Xo, 


Ao = wavelength in free space. 


(1.1-5) 


oe (1.1-5) arises when the typical term in (1.1-3) is substituted for A 
in the equation 
aA ,dA , OA 
Ox? © Ay? a? 
This and a similar equation for B are the forms assumed by (1.12) for the 
rectangular coordinates of our straight guide. 
The electric and magnetic intensities in the guide are given by 


=¢'A (1.1-6) 





2 2 
een i ae oe 
twe OXOY =: OB Oz | twp Oxdy 
2 
1 — 
Rear el. gee ge pls Gal, 
iwe Oy" iwp dy? 
pets OB Cy 8A 1B 
we O20Y = OX Ox = tw, O20 


each follow from (1.1-1). 

It is seen that the wave is completely specified by the gn1’s and df ns 
These may be arranged as (infinite) column matrices in any cnirehicnt 
order. Thus in dealing with (1.1-3) and (1.1-4) we may write 


fio] dh 
B20 | : e . 
gt = | op] am ee (1.1-8) 
§30 | 
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In our work we shall consider only those modes corresponding to a fixed 
value of m (or of 2) and the order is almost automatically fixed. | 

The factors which determine the propagation of the typical terms in the 
summations (1.1-3) and (1.1-4) for A and B are | 


Omn(3) = Bi eo Brinl2) = aa ee , (1.1-9) 


The column matrices obtained by arranging these quantities in the same 
order as in (1.1-8) will be denoted by a(z) and B(z). We may write 


a(z) = eM gt, Be) = oF dt (1-10) 


where exp(—2zI',) and exp(— I's) are square matrices defined by power series 

each term of which is a square matrix: 

= 38F as aT et" 
1° 2h 3b 


_—2.r 
e° =I] 


“ores (1.1-11) 


/ is the unit matrix and I’, is the diagonal matrix* 


‘Ty  O 0 ; ] 
_ |0 To =O - 
Pale : Ta: | (1.1-12) 


in which the order of the diagonal elements is the same as the order of the 
elements in the column matrix gt. Similarly I'g is a diagonal matrix whose 
elements are I'o1, T'o2, P11, T's, -- + , the order being fixed by dt. When TI’ is 
replaced by Ig in (1.1-11) it is easy to obtain r2, Ty, etc. and sum the 
‘resulting series to obtain . 


re #l io 0 | 
2, 0 gras 0 ; 
gle i ae, /(1.1-13) 


L 


A similar expression exists for exp(—zlI's). The expression (1.1-10) for 

a(z) is seen to be true when the square matrix (1.1-13) is multiplied, by 

matrix multiplication, into the column gt. 

_ It turns out that the field in a circular bend (in a rectangular guide) may 
be represented by a generalization of the foregoing expressions. In this: 

generalization, which will be studied in the following sections, the square 

matrices I, and I'g no longer have the simple form of diagonal matrices. 


* That is, a square matrix in which all of the elements other than those in the principal 
diagonal are zero. 
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1.2 Propagation in a Circular Bend 


In dealing with a circular bend we choose cylindrical coordinates (p, ¢, ¥) 
as shown in Fig. 1. With these coordinates we associate new coordinates, 
shown in Figs. 1 and 2, (x, y, s) which have approximately the same signifi- 
cance as in the straight guide. is the distance measured along the axis of 





Fig. 2 


the guide (defined as the locus of the centers of gravity of the transverse 
cross-sections of the guide), and x and y are the transverse coordinates. 

Let p = pi = (p2 + ps)/2 = pe + a/2 be the radius of curvature of the 
guide axis, and let the origin of the polar coordinates be taken at the center 
of curvature. Then z is equal to —piy where the minus sign is necessary to 
make (x, y, 2) a right-handed coordinate system. Since the vertical (in 
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Fig. 1) walls are to be specified by x = Oanda = aweset x = p—p1 a a) 2, 

Thus, the two sets of coordinates are related by ; 
p=xXx+pi— a/2=x+ po 
g = — 2/p (1.2-1) 
ony 

where pi, p2 and a are ‘constants. 


We again choose Ajand Bin (1.1-1) to be parallel to the y axis. In the 
cylindrical coordinates, 








eae See are a) y =194 1 6B 
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where now, from (1.1-2), A satisfies the wave oe 
1 0| OA 10°A 2 
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and likewise for B. 
One method of dealing with (1.2—3) which is sometimes used is to assume 


A = e'?? X (sine or cosine function of y) X f(p) (1.2~4) 


where f(p) turns out to be a Bessel function of order with its argument 
proportional to p. However, we shall proceed in a different direction. 
The change of coordinates (1.2-1) transforms (1.2—2) into 
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0 
een ry rane he eee i (4,2-5) 
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and (1.2-3) into 
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where p; is a constant and p = «+ p; — a@/2 is to be considered a function of 
x. To solve (1.2-6) and the corresponding equation for B we assume 


A = > a@m,(2) sin (mmx/a) cos (rity/b) 


(1.2-7) 
m= 1,2,3,---; n= 0,1,2,-:- 
B=)» Bmn(3) cos (nmx/a) sin (riy/b) 
bial _ (1.2-8) 
m= 0,1, 2,---; n = 1, 2, 3,--- 


these expressions being suggested by (1.1-3) and (1.1-4). . The expressions 
(1.2-5) for the electric intensity show that this choice of A and B make its 
‘tangential component vanish at the walls of the guide. Thus the boundary 
conditions are satisfied. | . 

In order to determine dmn(z) so that the differential equation for A is 
satisfied, we substitute (1.2-7) in (1.2-6). The resulting left hand side of 
(1.2-6) may be regarded as a function, say f(x, y), of « and y with the a’s 
and their derivatives entering as parameters. We must choose the a’s so 
as to make this function zero. Relations which must be satisfied by the a’s 
may be obtained by expanding f(x, y) ina double Fourier series for which the 
typical term is a coefficient times sin (rmx/a) cos (rny/b), and then setting 
the coefficient of each term to zero. This form of expansion is suggested by 
(1.2-7). However, it should be mentioned that such an expansion is best 
suited to a function which vanishes at « = 0 and x = a, a condition not 
fulfilled by f(«, vy) because of the term p 0A/dx in (1.2-6). This causes no 
real trouble because our region of representation runs only from + = 0 to 
« = a and hence our series is no worse than the Fourier sine series for the 
periodic function (of period 2a) which is —1 for —a < x < Oand + 1 for 
O<x< a. 

To carry out the procedure outlined above, we multiply (1.2-6) (after 
putting in (1.2-7)) by sin(rpx/a) cos(afy/b) and integrate « from 0 to a 
and y from 0 to b. Using the expression (1.1-5) for "jn and reducing gives 


— re Ant (z) + 2X [Pim Oy, 0 (2) = Sain amt (z) | = 0 (1 2-9) 


where p may have any one of the values 1, 2,3, -- - and the double prime on 
denotes the second derivative with respect to z. The P’s and S’s are 
constants given by 


Pum = (2/4) [ (pi/p’) sin (rpx/a) sin (wma/a) dx, — (1.2-10) 


Som = —20ma~ | sin (rpx/a) cos (xmx/a) dx/p —_ (1.2-11) 
0 
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‘The evaluation of these integrals is discussed in Appendix I. Thus (1.2-9) 
is the p™ equation of a set of differential equations to be solved simul- 
taneously for az¢(z), aet(z), --°. 

The customary method of solving a set of equations such as (1.2-9) is to 
assume that all the a’s vary as e”’ so that for each ane(z) we may write 
gmt. This leads to a set of simultaneous homogeneous linear equations 
for the constants gnf. In order that these equations may have a solution 
the determinant of the coefficients must vanish. Since the only derivative - 
of ant(z) contained in (1.2-9) is the second, y appears in the determinant 
only as vy”. Let yi : ve ; v3 y °++ be the values of + which cause the deter- 
minant to vanish and let ky;, ke;,--+ be the values of gi¢, get, +++ cor- 
responding to y’ = yj. The #’s are determined to within an arbitrary 
multiplying constant which, for the sake of convenience, is chosen so that 
Ry; = 1, 

Thus one solution of the differential equation (1.2-6) is 


A = & 7 cos (rly/b) >) Rng sin (rma/a). (1.2-13) 
m1 


This particular solution corresponds to the j™ one of the modes (traveling in 
the positive z direction) for which A is proportional to cos (rfy/b). 

In much the same way it may be shown that the series (1.2-8) assumed 
for B is a solution of equation (1.2-6) (with A replaced by B) provided the 
coefficients Brn(2) satisfy the set of equations 


Tt Bot) + 2. (omBnt(e) — UpmBne(@)] = 0 (1.214) 
for p = 0,1, 2,-+: and f= 1,2,3,+++. Here 


Onn = (€2/a) [ (pi/p') cos (wpx/a) cos (wmx/a) dx  (1.2-15) 


U om = TMéepa i : cos (rpx/a) sin (rmx/a) dx/p (1.2-16) 


where «9 = 1 and ep = 2 for p > 0. These integrals are discussed in 
Appendix I. 

The problem of determining the reflection from a bend in a wave guide 
involves considerable manipulation of equations (1.2-9) and (1.2-14).. The 
introduction of matrix notation in the manner suggested by the work of 
Section 1.1 for straight guides simplifies this work. Although am,(z) and - 
Bmn(%) are no longer the simple exponential functions given by (1.1-9), it 
turns out that the column matrices a(z) and Q(z) are still given by (for a 
wave traveling in the positive z direction) by the matrix expression (1.1-10). 
As mentioned earlier, I’, and I'g are no longer simple diagonal matrices. 
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We now turn to the task of expressing (1.2~9) and (1.2-14) in matrix 
form. 


1.3 Propagation Constani Matrix for Curved Rectangular Guide 


From this point onward in our investigation of propagation in the rec- 
tangular guide we shall assume A to be proportional to cos (rfy/b). Thus 
instead of the general expression (1.2—7) for A we shall deal with the more 
restricted form 


A = cos(rly/b) >) ane(z) sin (rmx/a) (1.3-1) 
m=1 
where ¢ has one of the values 0, 1, 2, 3, . Since the most general dis- 


turbance may be obtained by the superposition of disturbances of the form 
(1.3-1) no real generality will be lost. 

‘The introduction of (1.3-1) is suggested by the fact that the set az¢(z), 
aet(z), +++ may be determined from (1.2-9) (at least to within arbitrary 
constants of integration) without considering the other amn(z)’s, 1 ~ ¢. 
The introduction of (1.3-1) is also suggested by physical reasons. The 
plane of the bend is the z, « plane and there is nothing in the system tending 
to change the field distribution in the y direction. 

Equation (1.2-13) is a special case of (1.3-1). Furthermore the most 
general form of (1.3-1) (corresponding to a wave progressing in the positive 
z direction) may be obtained by multiplying (1.2-13) by an arbitrary con- 
stant c; and summing on j. 

In order to write the set of differential equations (1.29) for the Ame(Z)’s in 
matrix form we introduce the infinite matrices 


Te <0 0. i]... ee 
O [xe O : aot (z) 


FU ge ao Sr lh oz) = 1 ole) 
. . . . (1.3-2) 
Pa. Pes Su Su 

P= 


Po Po + ; S =} Sa Soe 


where the elements of I'p are obtained by setting 7 = in equation (1.1-5) 
for Pmn, and the elements of P and Q are given by the integrals (1.2-10) and 
(1.2-11). The rules of matrix multiplication then show that (1.2-9) is the 
p" element of the matrix equation . 
Pol'(z) — (T5 + S)a(z) = 0 (1.3-3) 
Premultiplying by F-! converts this equation into 
a’’(z) — T% a(z) = 0 (1.3-4) 
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where 
Poe (1G aS) (1.35) 


It may be verified by direct differentiation of the series (1.1-11) defining 
exp (—2zI’) that a solution of (1.3-4) is 


a(z) =e 7'% gt (1.3-6) 


where, as in (1.1-10) for the straight guide, gt is a column of constants (of 
integration). However, now I’, is to be obtained by taking the square root 
of the right hand side of (1.3-5), a process which is not easy since it usually 
requires one to obtain the characteristic roots and modal columns of T°, 
(see equation (1.3-10)). 

As far as (1.3-6) being a solution of the differential equation is concerned, 
I’, may be any matrix whose square is given by (1.3-5). We shall restrict 
it as follows: When £ = a/p1 becomes small, as in the case of a gentle bend, 
it is seen from (1.2-10) and (1.2-11) that P approaches the unit matrix and 
S approaches zero. Hence, I’, approaches the diagonal matrix >. I'a is 
chosen so that it approaches Io, that is, all of the elements in the principal 
diagonal are either positive real or positive imaginary. This makes 
exp(—zI'.) approach the diagonal matrix exp(—zIo). With this choice 
of I, the expression (1.3-6) for a(z).corresponds to a wave traveling in the 
positive z direction. 

The various modes of propagation in the bend may be obtained from I’% 
by expressing, in matrix notation, the steps leading to (1.2-13) (which gives 
A for the7 ** mode). We assume a(z) to be the column matrix obtained by 
multiplying the column matrix g of constants by the scalar quantity e””. 
Setting this in (1.3-4) gives 


WI-Tajg=0 |. (1.3-7) 


where J is the unit matrix. In order that (1.3-7) may have a solution, the 
determinant of the coefficient of g must vanish. This leads to the char- 
acteristic equation™ for +: 


ly — 1| = 0 (1.3-8) 


The vertical bars denote the determinant of the inclosed matrix. The roots 
Vi, 'Y2) °° are therefore the latent (or-characteristic) roots of [%. If we let 
k; denote*™* the column g obtained when y = 7; in (1.3-7) then 


* See Section 3.6 of Reference?. 

‘** We choose this notation in order to adhere as closely as possible to that of Refer- 
ence’, Incidentally, the column k; is proportional to the jt® column of «7 where x is 
the modal row matrix introduced in Section 5.1. . 
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(71 aaa) kj=0 | . (13-9) 


and the elements £1;, o;, -- + of the modal column k ; are the ones appearing 
as coefficients in (1.213). 


Equation (1.3-9) and the methods of matrix analysis lead to 
| =hylak', Ta=kblak? (1.3-10) 


where k is the square matrix whose j™ column is &; and [y’]a, [y]a are diagonal 
matrices having vj) 7; a8 the 7 * elements in their principal diagonals. The 
representation (1.3-10) certainly holds for the rectangular guide since in 
this case no repeated roots occur. 

In analogy with the expression (1.3-1) for A we shall henceforth deal 
with B in the form . 


B = sin (rty/b) 3 Bmt(z) cos (rmx/a) | (1.3--11) 


where ¢has one of the values1,2,3,---. In much the same way as before it 
may be shown that for a wave traveling along the bend in the positive 
direction the 8,,¢(z)’s in (1.3-11) are .given by 


B(z) = e 78 dt (1.3-12) 
where dt is a column of arbitrary constants and 
| =O +0) L813) 
In (1.3-12) and (1.3-13) 

Te O O : . Boe(2) 

=F oO me | =| es 
> Ge ) (1.3-14) 

Qoo Qu - b Un Ur 

U 0 Un Ur 


Q=1Q0 Qu-- 7 
35 oe a | 


where the elements of To, Q and U are given by equations (1.1—-5), (1.2- ie 
and (1.2-16), respectively. 


1.4 Continuity Conditions at Junction of Straight and Curved Rectangular 
Guides 


Electromagnetic theory requires that E,, H,, H, and H, be continuous in 
crossing a plane g = constant which marks the junction of a straight and a 
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curved wave guide (of the same cross-section). Comparison of the first 
equation in (1.1-7) with the first equation in (1.2-5) shows that E, is con- 
tinuous if (1).4 is continuous and (2) if 9B/dz in the straight portion is equal 
(the equality being taken at the junction) to (9;/p) 0B/dz in the curved 
portion, Examination of the expressions for the remaining field com- 
ponents shows that all the continuity conditions are satisfied if, at the 
junction, 


[A in straight portion] = [A in bend] 


oA ce 66 a3 PL JA (14-1) 
| 4 = Ee [7 in bend | 


and likewise for B. 

Let A in the bend be given by (1.3-1) and let a(z) denote the column 
matrix of coefficients shown in (1.3-2). A in the straight portion may be 
represented in the same way except that a(z) has a simpler form as explained 
in Section 1.1. When these expressions for A are inserted in (1.4-1), both 
sides multiplied by (2/a)sin(apx/a) after cancelling out the cos (rfy/b), and 
the results integrated with respect to x from 0 to @ we obtain relations which 
may be expressed as the matrix equations 


[a(z) in straight portion] = [a(z) in bend] 


a lhe 





where V is the square matrix whose p™ row and m™ column (p,m = 
1, 2,3, ---) is 
Vom = (2p1/a) [ sin (rpx/a) sin (rmx/a) dx/p, (1.4-3) 
0 


p being equal to p1 + « — a/2. 

By using expression (1.3-11) for B in the continuity conditions, it may be 
shown in much the same way that the column matrix @(z) given by (1.3-14) 
must satisfy the relations 


[G(z) in straight portion] = [8(z) in bend] 

. 1.44 
EZ ook |- E Be) «ne na | oe) 
dz dz 
where W is the infinite square matrix 


Woo Wo1 
W = Wo Wi ° i (1.4-5) 
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whose elemenst are 
W om = (€ppi/a) : cos (arpx/a) cos (rmx/a) dx/p 
0 


eo = I, ¢>= 2 for p>O 


(1.4-6) 


Both Vz, and W ym are discussed in Appendix I. 
PART II 
THEORY FOR A GENERAL WAVE GUIDE 


21 Mairix Propagation Constant for a Curved Wave Guide of Arbitrary Cross- 
Section - 


In Section 1.3 it has been shown that for a curved rectangular wave guide 
there’ exists a square matrix I’, (or I's) which plays the same role in the 
propagation of a wave consisting of many modes as does the propagation 
constant in a simple transmission line. There I’. was obtained from a 
special form of the wave equation which is suited to bends in rectangular 
guides. Here we adopt a different approach with the idea of showing that a 
matrix propagation constant I’ exists under more general conditions. 

The general theory of wave propagation in tubes shows that a wave 
traveling in the positive z direction may often be represented as 


p~ = 2d Gj e iolx, y) (2.1-1) 
where & is some quantity associated with the field and is analogous to the 
functions A and B of Part I. In (2.1-1) « and y are transverse coordinates 
and zg a longitudinal coordinate. y; is the propagation constant for the gn 
mode and ¢;(x, y) the corresponding eigenfunction. Fora circular bend ina 
rectangular wave guide 9,(x, y) is a combination of trigonometric and Bessel 
functions and y; is proportional to the order of the Bessel functions. 

We assume that we may find a set of functions 6,,.(x, y), m = 1,2,3,... 
such that every ¢;(x, y) may be represented as 


| o;(x, y) = > Ring Om, y) (2.1-2) 


The usefulness of this procedure depends upon our ability to pick a system 
of 0(%, y)’s which is appreciably simpler than the system of ¢;(x, y)’s. In 
the work of Part I 6,.(«, y) was taken to be the eigenfunction of the typical 
. mode of propagation in a straight guide, i.e. the product of a sine and a 
cosine. 

We assume further in (2.1-2) that the square matrix k exists where Ry; is 
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the element in the mm‘ row and je column of k; ie. if a root 7; is repeated, 
say, s times there are s linearly independent columns (k;’s) corresponding 
toy;. Substitution of (2.1-2) in (2.1-1) gives 


& = s: Om (x, ) 3 Ring Gp @° 77 
— ae (2.1-3) 
= a Hm ()On (x, 9) 
where - 
Lm(z) = rE Rng 6@.77 (2.1-4) 


Since 0n(x, y) is analogous to the product of the trigonometrical terms in 
(1.3-1) or (1.3-11) these equations show that y,(s) plays the same role as 
Omnt(z) or Bnt(z). Therefore, in accordance with the discussion given at the 
beginning of this section, we wish to snOW that the column matrix p(s) 
(which is similar to a(s) or B(z)) whose m™ element is pn(2) may be ex- 
pressed as 

pe) = et ft (2.1-5) 


In this equation I’ is a square matrix to be determined and Lo is a column 
matrix of constants similar to gt or dt. 
The rules of matrix multiplication and equation (2.1-4) show that 


p(s) = ke Jac ~ (2.1-6) 


in which [exp (—zy)]a is a diagonal matrix having exp (—sy,) as the j” 
element in its principal diagonal and ¢ is the column matrix formed from the 
c;’s. We introduce the column f+ by defining it as »(0) whence 


ft = ke, C= kt ; (2.1-7) 
maeee sk from (2.1-3), the value of Patz = Vis : 
=p? fi Om(x, ¥) | (21-8) 
where fz, is the m*” element in f*. 
From (2.1-6) and (2.1-7) 
w(s) = Rea kf (2.1-9). 


Tn this eaudtion k [exp (—2y)]a kis a square matrix which may be expressed _ 
as 








3 oe ye Ryd aa pa ail (kl-y] lee" 


—2zr 
= ¢€ 


(2.1-10) 
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Here [y]4 represents the diagonal matrix having ys the j*" term in its prnepat 
diagonal and 


= Rly]a k* (2.1-11) 


Therefore + we have shown that (zs) is of the form (2.1-5) which is what we 
set out to do. 

It is rather difficult to compute I’ from (2.1-11) using only the above 
definitions of k and y; for one must first obtain the functions ¢,(x, y). In 
dealing with the rectangular guide it is easier to use equasons (1.3- ry and 
(1.3-13) to determine I. 


2.2 Reflection at a Single Junction 


Let a straight wave guide extending from zs = — « tos = 0 be joined toa 
curved guide of the same cross-section which extends from s = 0 toz = ~, 
Let an incident wave . 


foe] 


@; = Dy lim ©" On, (22, y) (2.2-1) 
come in from the left along the straight guide. ‘The /,,’s are given constants, 
the 6,,’s are the modal propagation constants for straight guides (for Sts 
tangular guides they are the I’, ;’s given by (1.1-5)), and @n(x, y) is the mz" 
eigenfunction for the straight guide (the product of a sine and a cosine for a 
rectangular guide). _ 

What are the reflected and transmitted waves set up: by (2.2-1)? The 
reflected wave is of the form 


& =D fae One, 9) (22-2) 


m=1 


where the fm’s are to be determined. 
If we assume the representation 


fo) 


6 = py Hm (2)Bm (2, y) (2.2-3) 
to hold for all real values of z then, since 6 = 4; + ®, for z < 0, equations 
(2.2-1) and (2.2—2) show that . 

lm(Z) = Ime?" -+ fn e°™, m= 1,2,3,...32<0 (2.2-4) 


Introducing the column matrices p(z), h, f— and the diagonal matrices 
exp (=2I'o) where Ip is a diagonal matrix having-6,, as the m'™ term in its 
principal diagonal enables us to write (2.2-4) as 


p(s) = eo oA+t re a 2<0 (2.2-5) 
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Equation (2.2-5) is more general than the expression (2.1—5) for u(2) in that 
it contains waves going in both directions, but is more special in that Toisa 
diagonal matrix. 

In the curved guide we take p(z) to be given by (2.1-5) , thus 


uz) =e ft, 2>0 (2.2-6) 


where I is a square matrix whose elements are assumed to be known and ft is 
a column matrix whose elements are to be determined along with those of f-. 

The conditions that the transverse components of the electric and mag- 
netic intensities be continuous at the junction of the two guides are assumed 
to lead to the requirements . , 


[u(z) in straight portion] = [u(z) in curved portion] 


ae ceil | re _ | (2.2-7) 
a u(z) in straight portion | =/|V 77 u(z) in curved portion 


where the quantities within the brackets are evaluated at the junction and V 
is a Square matrix whose elements are constants. When the curvature of the 
curved portion becomes small V approaches the unit matrix. For the 
problem at hand (2.2-7) may be written as 


[u(2)}—-0 = [u(2)]e-40 (2.2-8) 
d 
[Zuo] =v[sue| (2.2-9) 
in which the subscripts ¢ = —0, ¢ = +0 refer to the straight and curved 


portions, respectively, of the guide at s = 

The requirements (2.2-7) have been established for the rectangular guide 
in Section 1.4. Their form is also suggested by the conditions that the 
voltage and current be continuous at the junction of two transmission lines. 
Thus if we let u(z) play the role of the voltage, the current in the first line is 
— Zz" du(z)/dz and the current in the second is — Zz" du(z)/dz where Z; and 
Z, denote the distributed series impedances of the two lines. It is seen that 
this leads to scalar equations which look like (2.2-7), but now V denotes the 
scalar Z;/Ze instead of a square matrix. 

Setting the expressions (2.2-5) and (2.2-6) for u(z) in the conditions 
(2.2-8) and (2.2-9) gives two equations which may be solved simultaneously 
to obtain f— and ft in terms of 4,19, TP and V: 


h+f = ft 
(2.2-10) 
—Toh + Tof~ = —VIft . : 
f- = Cot VI) — VI)k (2.2-11) 


ft = Go + VI) 2P ch (2.2-12) 
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Since f~ and ft specify the reflected and transmitted waves, respectively, 
they give the answer which we are seeking. 

If the curved guide should extend from z = — © toz = 0 and the straight 
guide from s = 0 to z = © the response to an incident wave €*’ h coming in 
along the curved guide would be 


a2) =e hte ft, z2<0 


2,2°13 
p(z) =e? fF, 2 >0 
A procedure similar to that used above shows that 
—= —(Tot VI) (To — VIDA 
f= —(ot VI (Wo — VT) saat: 
ft = Got VF)? 2VTh 


where, instead of condition (2.2-9), we have used 


v|4 ne) |, - E ute) | - 2.2-15) 


2.3 Reflection Due to a Bend 


Let the guide be straight for —«0 <3 < —candforc <z< ©, and let 
these two portions be connected by a curved portion in which the longi- 
tudinal coordinate z runs from —c¢ to +c. As in Section 2.2 we take the 
matrix propagation constants for the straight and curved portions to be the 
square matrices Ip and I’, respectively, and assume an incident wave, 
specified by the column matrix /, to come in from z = — . 

The column matrix y(z) whose m” element appears as the-coefficient of 

Om(x, y) in the representation (2.2-3) for ® is now given yy 


uz) =e" hee fy, z<—e 
u(z) = (cosh zI')p + (sinh 2)q, —¢c<2<e (2.3-1) 
we) = eft, c<z . 


In these expressions f—, f+, p, g are column matrices which may be de- 
termined as functions of the known matrices Io, I, V and h by substituting 
(2.3-1) in the conditions (2.2-7) which must hold at the junctions z = —c¢ 
and ¢ = ¢. 

By straightforward algebra similar to that used for the analogous problem 
in transmission line theory we obtain 


eto eho gt = [-F + (VT tanh oP +1) Tod eek i) 
eto ge — glo rt — [7 + OT coth cI +1) Tol eh 
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In these equations the infinite square matrix tanh cT is defined as (sinh cI’) 
(cosh cI)“ and coth cI as its reciprocal. sinh cI’ and cosh cI’ may be de- 
fined as power series in cl’ and may be expressed as combinations of exp (cl’) 
and exp (—cI’). . 

An expression for the column matrix f~ may be obtained by adding the 
equations in (2.3-2). Before doing this it is convenient to introduce the 
two column matrices x and y defined by 


(Vel tanh cI + clo) x = Moe h 


2.3-3 
(VcI coth cI + cYo) y = Do ae ( 


where the scalar length c has been introduced to make the various terms 
dimensionless. Each equation in (2.3-3) represents an infinite set of 
simultaneous linear equations to be solved for the elements of x or y. 

Once x and y are known the reflected wave is given by 


f= EM (ety) — rh (2.3-4) 
and the transmitted wave by 
fr=éM@@—y -(2.3-5) 
PART III 


GENTLE BENDS—GENERAL THEORY 


3.1 Limiting Forms Assumed for T and V 


It will be shown in Part IV that for gentle circular bends in rectangular 
wave guides the matrix propagation constant I’ is such that 


r="%+F (3.1-1) 


where Ty is the square of the matrix propagation constant for the straight 
guide. I is a diagonal matrix having 6,, (which is one of the I',’s, depend- 
ing on the set of modes under consideration, given by (1.1-5)) for the m™ 
element in its principal diagonal. F is a square matrix of infinite order in 
which the elements /’;; in the principal diagonal are of order £ and the re- 
maining elements /’;;, i # 7 are of order &. Here & = a/p1 is the ratio of the 
guide width to the radius of curvature of the bend. As the bend becomes 
more and more gentle, & — 0. 

The asymptotic expressions given in Appendix I show that, for gentle 
bends in rectangular guides, the square matrix V which appears in the 
junction conditions (2.2-7) approaches a unit matrix as &— 0. In par. 
ticular Vj; = 1 + v1; where v;; is of order &, and V;;, the element in the i™ 
row and j™ column, is of order £ when i # j. 
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Throughout the remainder of Part III we shall assume that T?, F and V 
behave as mentioned above. In addition we assume that. there is no de- 
generacy, i.e. all of the 6,,’s are unequal to each other and to zero. 


3.2 Propagation in a Gentle Bend 


Here we assume that the elements of Ig and F in the expression (3.1- =) 
forl’? are known. We wish to find the modal propagation constant 7; and 
the corresponding eigenfunction ¢,(x, y) for the j* * mode. 

After squaring both sides of the collineatory Fede ctdrination 2. 1- 41) 
connecting I’ and the diagonal] matrix [y]¢ we obtain a relation which may be 
written as k[y?]z — T?k = 0. The left hand side is a square matrix having 
(yl — I°)k; as its j column. Here J is the unit matrix and &; is a column 
matrix having i;, ke;,.. . as its elements (k; is the j*™ column of &). Thus 
we have a system of simultaneous linear equations in which the coefficients 
are furnished by the square matrix y3J — I’? and in which the unknowns are 
hij, Raj, +++ . Accordingly, y} is the 7" latent root of I? and k; is its cor- 
responding modal column just as for the rectangular guide in.Section 1.3. 
In order to apply equations (A216) of Appendix II we set }; = 7; and 

u = Tso that, from (3.1-1), . 


‘ = 5 + Pi, hij Sa Fiz; 4A q (3.2-1) 

Therefore 
= 05+ Fiat 2a ee — 8:) (3.2-2) 

Ri; = I, ke; = i/(83 — SF q (3.2-3) 


where we have neglected terms of order £° in a and of order £'in k, % 
s ~ 7. The prime on the summation indicates that the term s = 7 is to be’ 
omitted. 

When f1;, ko;,...are known the eigenfunction ¢;(x, y) may be written 
as a series in 6,,(x, y) by means of equation (2.1-2). 

In Section 3.3 we shall need the form assumed by the square matrix I’ 
tanh cI’ in a gentle bend. This matrix is used in computing the reflection 
from such a bend, as might be inferred from equation (2.3-3). The formula 
to be used is (A2-18) with w = T?, A; = y; and with the elements of the 
square matrix & given by (3.2-3). In the diagonal matrix of (A2-18) we set 

riQy)) _ i tanh ny = 7; tanh CYyi = ¥ jt; 
(3.2-4) 
= tanh cy; 
and for the elements of k~! we use (A2-19) together with the line above it. 
When the three matrices on the right of (A2-18) are multiplied out the ele- 
ment (T tanh cI’) ;; in the 7” row and j** column of I tanh cI’ is found to be 
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(r tanh cl) i= (y it; ae ¥ ts)Raj, iF j (3.2-5) 


HG tanh V)as = Vil: + pay (y its = Vint) R imRmié (3.2-6) 


where in (3.2-5) and (3.2-6) terms of O(€) (“order of’) and O(é*), re- 
spectively, have been neglected. This is in line with the fact that the terms 
in the principal diagonals of our matrices must be accurate to within O(é?*) 
while the remaining terms need be accurate only to within terms of O(6).. 
The summation with respect to m runs from m = 1 to © with 
m = i omitted. : 


3.3 Reflection from a Genile Bend 


. When the bend is gentle so that V and I’ behave according to the descrip- 
tion given in Section 3.1, the matrix expressions for the reflection coefficients 
given in Section 2.3 may be evaluated. The results stated in Appendix IT 
for “almost diagonal” matrices furnish the principal tools for this work. 

It is assumed that the incident wave coming in along the straight guide 
from the left is 6; = exp(—26,) 0,(x, y) and hence contanis only the pt 
mode. Comparing this with the general expression (2.2-1) for ®; shows 
that hy» = 1, lin = 0, m ¥ p, and all the elements of the column matrix / are 
zero except the pt which is unity. 

We start by writing the first of equations (2.3-3) as 


(f tanh cP + V-Ty)x = VIM eh (3.3-1) 


Since V approaches a unit matrix as  — 0, the element (V-),; in the i* row 
and 7*® column of Vis 


(V)i3 = —Viz, 0G 


(V>)ie = 1 — 0 + DO! Vin Vins ee 
where V;; = 1+ 04:, 7,7 = 1,2,3,...and the summation with respect to m 
runs from 1 to ~ with the term for m = i omitted (as indicated by the prime 
on =). In omitting this term we are neglecting v7; because it is of order 
| These results follow from equation (A2~2) of Appendix II. As usual, 
the elements in the principal diagonal are accurate to within O(é?) and the 
remaining elements to within O(é). 

It follows that V-1T) e°*° h = 7 isa column matrix whose it" element is 


Ni = —Vip bp ee 1 F p 


: ; (3.3-3) 
Np = (1 — upp + Dy’ V pm Vip) bp, Baga 
: ™ 
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Likewise, the element in the it row and j** column of the square matrix 
V—Tois —V;;6; when i ¥ j and is ; 
(1 — v6 + D0’ Vi Vins) 85 (3.3-4) 
when 7 = 7. 
By combining the approximate expressions for the elements of I’ tanh cI’ 


(given by (3.2-5) and (3.2-6)) and V-! I'p we find that if w,;; denotes the i** 
row and j** column of I tanh cI’ + V-!T) then. 


Us; =Di-—Vgd;, iF i 
us = iti + Do’ Daikim + (1 — 1% + De Gas) (3.3-5) 
= 0 — 830 + Dy (Danskin + 8:Cni) 


In these equations we have set 
04 = OF + Viti, Cmi = VimVmi 
Di = (Viti — Ymbm) nt = (Viti — Ymtm) Fm i/ (8% — 8m) 


where. y; and &,,; are given by (3.2—2) and (3.2-3). 

We are now in a position to identify the matrix equation (3.3-1) for x 
with the set of equations (A2-20). The quantity »; which appears on the 
right hand side of the i equation in (A2-20) is given by (3.3-3). The 
coefficients which appear on the left hand side are the w’s defined by (3.3-5). 
Therefore, from (A2~21), when i ¥ 9, 


(3.36) 


) 
5,¢° 
Xi — 





[—Viplop . bp) Dip] 
Tip ; 


(3.3-7) 
[Vip bp lp a Fin (6; i Op ty) (5; ae 4) ] 


oP 


(1 FAL +) 
where we have neglected higher order terms and in so doing have replaced 
y; by the simpler 6;. 

When i = p (A2-21) yields 


Xp = Up alt - De (tepm Ump Np — Upm Upp tim) / (Umm Upp) (3.3-8) 


In order to combine the second order terms in 1/1», with those in the rest 
-of the expression for «, we assume that o, is the major portion of up,. 
Then, approximately, 


1/tpp = o, ll ar Op Pop Fp a Poe form Ron + 05 Gan Gal (3.3-9) 
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This assumption, which is equivalent to assuming that tanh cy, differs 
appreciably from —1, does not appear to restrict our results since tanh 
CYp is either purely imaginary or else real and positive. 

Substituting the appropriate values in (3.3-8), neglecting higher order 
terms, and using the definition (3.3-6) for Cup leads to 


Xp = Op * bp é *>[1 a (1 = 5p/op) pp = 2 (op om) oF Cnp(op — 5p) (Om — 6m) 


a Pi — Om Rem) + (cp = 8p)Dom View = 6m Dinp Von (3.3-10) 


A reduction similar to that used in going from the first to the second line 
of (3.3-7) alae our final expression for xp 


Ope? ine? 

_ m V mp bm bm) 

0 Fh ark aes 
mae 


+ 2 Seas ers Yb 2 Paltz (3.311) 


_ Vinp bp Pym + Pom Fp (6p i Bin /ty) (mn i= a4] 


The above expressions for x; and «, have been derived from the first of 
equations (2.3-3). The second of equations (2.3-3) determines the column 
matrix y in the same way that the first equation determines x except that 
coth cI’ now replaces tanh cI’. Therefore, we may obtain expressions for 
the elements of y by replacing the /’s (where t; = tanh cy;) by their recipro- 
cals in the expressions for the corresponding ~’s (i.e. in (3.3-7) and (3. 3-11). 
The values obtained in this way lead to, when 7 + p; 


fi = (a + 9) 
= sy ei sinh cy; + Yp)[—Vip dp — Pag) (bi + 5,)] 
fi = @@—-y) 
me 8p OHH e-—® sinh cy; — Yn)[Vin bp — Fip/(5i — 65)] 


where we have used the expressions (2.3-4) and (2.3-5) for f- and tbe 
When 7 = 9, 


fp = O(a + yp) — € 
where . 


Ay = 2typ + (yp — 89)8p° — 2, | Vn Vanp + 


™ 


: 


(3.3-12) 


6 


> = —e"-7»)[Ay (sinh 2cyp)/2 + As] (3-413) 





Vom Pap + | 
bn — Op 
(cosh 2cyp — e~?7™) 
As = ay 26 5 (8 So o,) [V wil On + V. mp Fm 8» + Lm Fup) 


m 





REFLECTIONS FROM CIRCULAR BENDS . 327 


Q 


The expression for {7 may be obtained in the same way but it is slightly 
more complicated. | 


fh = (ay — yp) = [1 — Ay + Ay (sinh 2cyp)/2 + As] (3.3-14). 
where 


As = (1 — & ?)(yp — 8)?/(465) 


VomE mp — Vino F. 2P om tt 
/ —2¢Ym 7 pm+ mp -mp+ pm pm+ mp 
= Die |- Vom View 1 in — dp T @ — 32) i 
(e-°!™ cosh 2cyp — 1) 
Pa ig cee see saisicsaenaes Sam 
; -5 28m Sp(8m — Sp) 
(mm at 5; 3) on Pa 
in 5, 





i | Vom Fao p= Vero Pond a 


There are several points we should mention about these formulas for fpand 
ft: The summations with respect to m run from 1 to © with the term 
m = pomitted. y,; and 6; are the propagation constants of the 7th mode 
in the bend and in the ‘eLrabht portion, respectively. The difference. 

= 9, py may be expressed in terms of the F’s by equation (3.2-2).. In the 
course ‘al obtaining (3.3- re and (3.3-14) relations of the following sort. 
were used. 


tpl + 1,7? = ¢€ (sinh 2cy p)/2 
(lm + p)(L bob) “(L fn)" = "(cosh 2eyy — &°"™) 
The term A3 arises when we subtract ( — £,) (1 +#,)> from | 


5p/(Sp + Yotp) — 6 ote/ (Vp +6 pty) 


Since yp — 5p is O(e?) for a circular bend in a rectangular guide (y»— >)" 
is O(é*) and hence 4A; is negligible in the cases we shall consider. _ 
The reflected wave set up by an incident wave of unit amplitude and con- 
taining only the p™ mode (i.e. the incident wave described at the beginning 
of this section) is given by the column matrix f~ whose elements may be 
“ obtained from (3.3-12) and (3. 3-13). Likewise, the. transmitted wave is 
given Dy fr. 


PART IV 


GENTLE CIRCULAR BENDS. IN RECTANGULAR WAVE 
GUIDES 


4.1. Propagation of Dominant M ode in a Gentle Bend—H in Diane of Bend 


When the magnetic intensity H lies in the plane of the bend, H, = 0, 
and equations (1.2-5) show that B = 0. Thus we have to deal only with 
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A. In order to study the dominant mode we set ¢ = 0 in the cos (rfy/6) 
(A depends on y through this factor) in the formulas of Section 1.3 which 
involve A and assume the dimensions of the guide to be such that b < a. 

We wish to determine 73, the first latent root of I’, defined by (1.3-5), 
from the approximate formula (3.2-2). In our case the elements dx of diag- 
onal matrix Ig are obtained by putting n(= £) to zero in (1.1-5): 


d, =Tro=o + (am/a)’, m=1,2,3,-°° (4.1-1) 
‘so that (3.2-2) becomes . 


i = Tio + Fu — Dy Pim Fmt an (mM = 1 (4.1-2) 


The first task is to find the elements of the matrix F where, from (3.1-1) 
and (1.3-5), 


F=M%—-T =(P"°-DM+PS (4.1-3) 


In the case under consideration P = IJ + R where R is a square matrix 
whose elements are very small. In fact, the asymptotic expressions lead- 
ing to (A1—18) show that R;; and S;; are O(&), with € = a/p1:, while Rj; 
and S;; are O(£) ifi +7 is odd and O(€) ifi-+jiseven. When the approxi- 
mate value of P- obtained from (A2—2) is set in (4.1-3) and the matrix 
multiplications carried out it is found that 


Fis = — Rio + Sis + O/) 


c0 | co (4.1-4) 
Fi = (—R + >) Rim Rei) Po tSis DS Ra Sar OE) 


- The ‘order of’? symbol O( ) will be omitted in the following equations, it 
-being understood that the terms in the principal diagonal are correct to 
within O(£) and the others to within O(£). 

The values of the F’s which enter (4.1-2) may be computed from the 
asymptotic expressions (A1~18) for the R’s and S’s. They turn out to be 


Fim = —4&m[40 for (m? — 1)? + 3a7(n? — 1)7] 


I 


Fm = —4€m[40 ior 7m? — 1)? + am? — 1)7)] (4.15) 
Fu = &(p(1 — 62”) + 6a-*)/12 


In the expressions for F';, and F,1 m is supposed to have the values 2, 4, 6, 
For odd values of m Fi, and Fini are O(). When i = 1 in the ex- 
pression (4.1-4) for fi, the two series therein reduce to S3 and S, where 
So=. Dd wn’? -— 1)” (4.1-6) 


m=2,4,6++° 
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By expanding the typical terms in partial fractions and using the fact that 
the sums of (see, for example, page 238 of Reference!) 


Up= 1+ 374+ 5%4+- 


(4.1-7) 
= (—)° 42" — 1)Bgr"/q! 
for q = 2, 4, 6, are 7’/8, 7°/96, 7°/960, it may be shown that 
Ss = 1/64, Sy = 2°/768 — 2/128, ies 
Ss = (15¢° — 2°)/3072. 


In (4.1-7) B, denotes the g™ Bernoulli number. The values of S', may also 
be computed in succession from the two relations* 


Dp , pti 
27-1 , 21-2 
Ui es CHS a LD Co aS 
t=1 i=1 


where Cm,, iS a binomial coefficient. Still another method is to make 
use of the generating function 


DP Sen = A +2) ps: (m? —1— #7 = 1 — kaxcot ax 
p=0 M==2,4,6° ++ 


where 4x? = 1-+ #. Note that by this definition S; is 3 in contrast to the . 
non-convergent series obtained by putting p = 1 in (4.1-6). 

Substituting the values for the F’s given by (4.1-5) in the expression 
(4.1-2) for yj and using the sums (4.1-8) of the series which occur gives 
| , : Are" 2 

vi = Th fae a Tio(1 — 6m) + (aTio/2)"(5 — '/3)]  (4.1-9) 
When the dominant mode is propagated without attenuation both yj and 
Tio are negative. 

The general form of (4.1-9) has been obtained by both Buchholz! and 
Marshak’ by different methods. In our notation their result is 

2 4 
Ynn = Tan — =. E + a'Tnn(l — 60m”) + (=) (5 — xn’) | 
4a Tm 
(4.110) 


where ‘Ymn is the propagation constant for the m, nt* mode when the mag- 
netic vector is in the plane of the bend. 


*T am indebted to,John Riordan for these relations. 


330 BELL SYSTEM TECHNICAL JOURNAL 


4.2 Reflection Due to Dominant Mode Incident upon Gentle Bend—H in 
Plane of Bend 


Let the system be the one described in the first paragraph of Section 2.3 
and let the incident wave contain only the dominant mode. Then the 
matrix propagation constant is the I’, of Section 4.1 and the column matrix 
h specifying the incident wave has unity for its top element and zero for its 
remaining elements, i.e., = 1 in the formulas of Section 3.3. 

We shall be interested only in the reflection coefficient, {1 of the dominant 
mode. Here we shall denote it by gio, in line with the notation of equa- 
tion (1.1-3), in order to distinguish it from the corresponding coefficient 
(which will be denoted by do1) when £ lies in the plane of the bend. 

Setting » = 1 in the expression (3.3-13) for the reflection coefficient and 
using equation (4.1--1) for 6, gives 


ft = gm =~" A, (sinh 2cy1)/2 + Aa] (4.2-1) 
where 71 has just been obtained in (4.1-9) and 
Ay = 20 + (yi — Tio) 


_ 2 | Vi Vg Be Vim Poa -+- Vint vi 


nw a-?(m? — 1) 


4 (4.2-2) 
tes cosh 2cy, — e °™" 
: mea2 2D no Tio 7 am? ao 1) 
° [Vir Prt Po as Vm Pin Tio ++ ee Fra] 
From (At-18) and Vin = 1 + on it follows that | 
a = ECL — 6n-)/12 
(4.2-3) 
Vim = Vint = 8a-2Em(m? — ‘2 
where m = 2, 4, 6,---. For odd values of m, Vin and Vmi are O(£?)° | 


Substituting these values together with those for the F’s given by (4.1-5)’ 
using the sums (4.1-8) and the expression (4.1-9) for va Tas finally leads 
to (after considerable cancellation) 


A, = — &Tha 7/4 (4.2-4) 
Likewise, for even values of m, 
VimE milo + VinaF mV 10 + FimEm1 = 16£%m2a-4(m? — 1)-? (4.2-5) 


All of the terms in the expression (4.2-1) for gig are now known (the values 
of Ym may be obtained by setting m = 0 in (4.1-10)). We shall make the 
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further approximation of putting Tmo for ym. Since Ino — Ym is O(€?) no 
serious error is introduced and we have . 





= 2 sinh 2cT'19 £28 (cosh 2cT 19 — e% Tmo) ne 
$10 422 ae > 2 re 
8P yo a T m=24,6-+- Tyo T,.0 @ (m2 — 1) 


(4.2-6) 


in which 
almo = [r'(m? — 1) + aid’, —E=a/p. (4.2-7) 


For frequencies such that only the dominant mode is propagated the 
ratio of the power in the reflected wave to the power in the incident wave is 
| gio |2.. Marshak has given an expression for this ratio which is the same 
as that obtained from (4.2-6) when the negligible (for his case) terms g 2ehmo 
are omitted. 

The corresponding expression for the transmission coefficient derived 


from (3.3-14) for ff is not as simple as (4.2-6). 


4.3 Propagation of Dominant Mode in a Gentle Bend—E in Plane of Bend 


When the electric intensity £ lies in the plane of the bend, £, = 0, and 
equations (1.2-5) show that Ad = 0. Here we deal with B in much the 
same way as we dealt with A in Section 4.1. The dominant mode is ob- 
tained by setting ¢ = 1 in the sin(rfy/b) in the formulas pertaining to B 
‘in Section 1.3. It is assumed that 6 > a. 

Examination of the matrices (1.3-14) indicates that, for the sake of con- 
venience, we should call the top row of our matrices the 0” row and the left- 
most column the 0“ column. In line with this we call yo the propagation 
constant of the dominant mode in the bend. The elements 6;, of the diagonal 
matrix Ig are obtained by putting n(= ¢) = 1 in (1.1-5): ? 


d= Ta =o + (rm/a) +7 /b, m=0,1,2,-++ (43-1) 

‘When we make the appropriate shift in the subscripts, equation (3.2-2) 
yields ¥ 
yo = Tu + Fo + pm Pom Fno@ 7 mt (4.3-2) 


in which the elements of the matrix F are to be determined from (1.3-13): 
F=%—-%=-(@' -DM+OU (4.3-3) 
As in (4.1-4) we have, withQ =J-+ T, 
Fy = —TTa + Vis 
| . | a (4.3-4) 
Py = (-7 + D0 Tin Ts) Ta t+ Us — DS Tin Uni 
m=0 


m=0 
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By using the asymptotic expressions (A1-19) for the 7’s and U’s and sum-. 
ming the series with the value of (4.1-7) for g =.4 given in Section 4.1 we 
obtain 


eee +a” 


Fom = 
Fno = —8&0 rm (4.3-5) 
Foo = £T01/12 

Jn these expressions m is supposed to have the values 1, 3, 5, +--+. When 


m is even Fom and Fino are O(E). 
Substituting (4.3-5) in (4.3-2) and summing the series with the heln 
of the values of (4.1-7) given in Section 4.1 gives 


vo = Toi — ETS + 20°T1)/60 (4.346) 


A result equivalent to (4.3-6) has been given by Buchholz who also gives 
the approximation to the propagation constant when m > 0 (and the elec- 
tric vector in the plane of the bend). In our notation his approximation 
is . 


ae = ie + - 2 [3 ae (a=) (10 + v > mm) 


Aa ; 
oa (as) (21 + x?m | 


In writing (4.3-7) we have corrected a misprint in Buchholz’s expression. 
In order to.agree with Buchholz’s equation (5.30a) the leading term within 
the square brackets would have to be changed from 3 to —3. This change 
was indicated by the results obtained when our equation (3.2-2) was used 
to obtain special cases of (4.3-7). Probably the best way of obtaining 
(4.3-7) is furnished by Marshak’s method (WKB approximation, out to 
second order terms, applied to Bessel’s differential equation). If one wishes 
to verify (4.3-7) by using Marshak’s report® as a guide, he should correct 
the misprint in Marshak’s equation (12a). 





(4.3-7) 


4.4 Reflection Due to Dominant Mode I neidenl upon Gentle Bend—E in Plane 
of Bend 


The problem here is the same as that treated in Section 4.2 except that 
now the electric vector lies in the plane of the bend. In line with equation 
(1.1-4), the reflection coefficient fi: for the dominant mode will be denoted 
by di. As in Section 4.3 the subscripts indicating the position of matrix 
elements will be adjusted so as to start with 0 instead of 1. The square 
matrix W given by (1.4-5), and associated with the junction conditions for 
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B in the same manner as V is associated with A, now replaces V. Thus 
our expression (3.3~13) for the reflection coefficient becomes 

fou => dy =) [Ai(sinh 2cT 01) /2 -+- Ag] . (4.4-1) 


where we have neglected the difference in I'o; and yo and where | 


io) 





At = 2uoo + (yo = ia _ » 
| Wo Wao + woe ea as re Fon 
mee (4.4-2) 
dice Ce) cosh 2cI'o1 — eel mi 7 


m=1 22 milo, 7? a? mi? 
- [Wom Fimo Ta + Wo Fom Var + Fom Fmol 


From Woo = 1 + woo and the asymptotic expressions (A1--19) it follows 
that, for m = 1, 3, 5, 


Woo = 7/12 
Wom = 2ém nr”, Wno = 4ém 3? 


For even values of m, Won and W..0 are O(€). Substitution of these values 
together with those for the Fs given by (4.3-5), using the sums (4.1-7) 
and expression (4.3-6) for vi, — Toi leads to 


A; = #/12 
WomF moma + Wino oml'01 + Fomlmo = — 8207 r om 


for m odd. 
Thus the reflection coefficient for the dominant mode when £ lies in the 
plane of a gentle bend of length 2c is approximately 


(4.4-3) 


—2c Dnt 


rc _& sinh 2cTo A E4T or > cosh 2cIo1 — e 


do1 aa 
24 | alee CE ee MAT m1 





(4.44) 
where T'n1 is given by (4.3-1) and 6 > a. 
PART V 


~ NUMERICAL CALCULATIONS 


5.1 Bend in Plane of Magnetic Vector 


Let a/b = 2.25 and Xo/a = 1.400 where Ap is the free-space wavelength 
of the dominant wave striking the bend. The propagation constant 
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Tyo of the dominant wave is obtained by setting m= 1, 2 = 0 in (1.1-5). 


The I’s corresponding to the men modes may be obtained from the 
same formula: 


: . 
= = + a — — 10.272 alo = 1 3.205 


2 742 
C = 
a lio ¥ 
aTeo = aio + 30 = 19.336 =. ang = 4.397. ~~ (S11) 
a'T3o = aT io + 8° = 68.684 also = 8.288 


We shall consider a 90° bend. The approximation (4.2-6) appropriate 
to gentle bends becomes 


gio = i€ [— .0122 sin (5.03/£) + .0087 cos (5. 03/8)] (5.1-2) 


where the exponential terms have been omitted since they are generally 
negligible. In (5.1-2), & = a/p: and the arguments of the sine and cosine 

terms arise from 2cl'9 = walo/(2£). From (4.1-9) the approximate 
change in the propagation constant produced by the curvature is obtainable 
from 


yi — Tio = .294#/a" (5.1-3) 


where 71 is the propagation constant of the dominant mode in the bend. 
The determination of giy by matrix methods will be illustrated ‘for a 90° 
bend in which p;/a = 0.6. This makes c/a = pyr/(4a) = .4712, C9 = 
71.510 and the appropriate equations in (2.3-5) and (2. 3-4) become, upon 
setting f1 = gio and ft 1 = gio ) 


oto = (x1 — 1) = (.061 + 7.998) (a1 — yr) 
gio = (ay + yn) — 7? = (.061 + 1.998)(m + 91) (5.14) 
Aie.O03205401 | 


Here «1, y1 are the top elements in the column matrices x, y. The problem 
is to compute « and y from the matrix equations (2.3-3) with I’ replaced 
by Ta, Io defined by (1.3-2) with ¢ = 0, andha column matrix whose 
elements are zero except the top one which is unity. Since the order of the 
matrices is infinite, an exact solution calls for an infinite amount of work. 
A compromise must be made between the accuracy desired and: amount of 
work one is willing to do. The following numerical work uses third order 
matrices. 

The first step is to compute the square matrix, obtained from (1. 3-5), 


a Te = P(e? + aS) | “(Gas 
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The elements of the diagonal matrix a?" are given by (5.1-1) and those of 
P and S by the equations and tables of Appendix I. 


1.4429. 8812 .6821 |? | —12.292 3.3447 1.6087 


ar’, = | .8812 2.1250 1.3745 — 6.3039 16.369 6.9738 
.6821 1.3745 2.4879 — 3.5034 —11.3031 65.213 (5.1-6) 
— 9.086 — 1.785 —5.178 
= 157 17.218 — 19.362 
.996 — 13.566 38.329 


The next step is to use (5.1-6) to evaluate the coefficients of x and y in 
(2.3-3). The square matrices Tc tanh I'ac and I'ac coth I'.c cause most of 
the computational difficulties. We shall deal with these matrices by using 
Sylvester’s theorem (an account of this theorem is given in Section 3.9 of 
Reference®). This requires the determination of the latent roots and modal 
rows of aT’. However, it is interesting to note that the matrices in ques- 
tion may also be computed from cT% (which is easily obtained from aT?) 
by processes which employ only matrix multiplication, addition, and in- 
version. 

Thus, setting A” for cT%, 


“lo: _ A’. At 
A‘sinh A = I+atyt ae 
AY AS 
cosh A he ge is 
A coth A = (cosh A)(A7 sinh A)! 


A tanh A = A%(A coth A)-. 


Although the series always converge, they do so too slowly to be of use in 
our computations. The same is true of the series 


Atanh A = >, 84° [(2m — 1 rl +4477". 
m=1 


For the matrices we shall encounter it appears best to use Sylvester’s 
theorem even though this requires the determination of the latent roots and 
modal rows of aI. The square matrix formed from the modal rows* 
will be denoted by «x. 


* As has already been mentioned in the footnote associated with equation (1.3-9), we 
shall use the notation and theory set forth in Sections 3.5 and 3.6 of Reference’. 
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We shall use the relations* 
T.c¢ tanh Tac = x [ye tanh yea x 
oo (5.1-7) 
Tac cosh ae = «x [yc coth yclak 


where the subscript d on the brackets stands for ‘‘diagonal’’ matrix, the 
i™ element in the principal diagonal of [yc tanh yc]¢ being yc tanh y;¢ where 


vic = ont*/a (5.1-8) 


and X; is the ith latent root of aT’, . In our applications 7; is either posi- 
‘tive real or positive imaginary. 


From (5.1-6) the \,’s are the roots of 


d + 9.086 1,785 5.178 
—.157 A — 17.218 19,362 (5.1-9) 
— .996 13.566  — 38,329 


=? — 46.461 2 — 101.96 -+ 3464.5 = 0 * 
and have the values | 
A1 = —8.886, Ae = 8.284, As = 47.06 (5.1-10) 


The elements x21 , K31 of the modal row [1, x21, x31] corresponding to \1 
may be obtained by solving the two equations derived from the last two 
elements of 


(1, Kor, Ksil(AW — a T2) = 0 (5.1-11) 
namely, 
1.785 + (Ai — 17.218) Koi + 13.566 x31 = O 


5.178 + 19.362 ker + (Ar — 38.329) x31 = 0 
When the value of \; from (5.1-10) is used these equations yield 
kor = .1593, k31 = .1750 
Likewise, the first and third elements of 
[kia , 1, ksel(Ael — a Ta) = 0 
and the first and second elements of 
[eas , kes, 1st — a TG) = 0 


* This is the modal row matrix analogue of equation (11) in Section 3.6 of the Refer- 
ence’, The modal rows of I. are equal to the modal rows of a*I.. 
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give . 

; Kiy2 = .0465 K32 = £6524 

Kig = .0165 Keg = —.4555 
Thus, the numbers entering (5.1-7) are 


yic = 4712 (—8.886)"? = 71404, yoo = .4712 (8.284)"7 = 1.356 


yic tanh yic = —8.382 yee tanh yoc = 1.187 

yic coth yic = .2354 yee coth yec = 1.549 

vac = 4712 (47.06)? = 3.233 4 1593 1750 
y3c tanh y3c = 3.228 k=] .0465 - 1 6524 
v3c coth y3¢ = 3.243 0165 —.4555 1 


For the purpose of calculation it is convenient to transform (2.3-3) by 
inserting (5.1-7) and premultiplying by xV-!. We obtain 


(fyc tanh yclax + kV Toc) w@ = KV "cD oe? *h 


a i 2 (5.1-12) 
~* (lye coth yclax + xV Toc) y = nV coe" 8h 
in which “s  6 

1 .1593 .1750 1.1204 3911 1629 |? 

KV = 0465 1 6524 3911 1.2833 4946 
00165 —.4555 1 1629 4946 1.3460 
9492 —.1992 .0883 Tyoc = 71.510 . 

=| —.2608 7686 = .2339 | , Toc = 2.072 

1427. —.7900 1.0160 T30¢c = 3.905 


where the elements of V are obtained from the formulas and tables of Ap- 
pendix I. 

The i” equation of the set obtained by writing out the first of equations 
(5.1-12) is _ 


[kjevic tanh yic + (xV>) iT joclx; = (KV) acne’ (5.1-13) 


Me 


P) 


1 


where (xV~) ;; denotes the element in the i” row and j™ column of «V-}, 
kj: 1s the element in the i row and j*" column (note the reversal of the 
usual convention regarding the order of subscripts) of x, xi; = 1, and k 
has disappeared because it is a column matrix whose top element is unity 
while the remaining elements are zero. It will be noted that the only 
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imaginary terms in (5.1-13) occur in the coefficients of > and arise from the 
imaginary quantity Tisc. By making the substitution 

—— a; Tye" 1 
. 1 + UCI 19 


the set (5.1-13) may be reduced to 


(5.1-14) 


3 F 
(avyi¢ tanh y;c)u1 + > [kzeyec tanh yuo + (KV) ij T joclu; = (KV) a 
j=2 
(5.1-15) 


in which the coefficients are all real. It should be noticed, however, that 
nothing is gained by making the substitution (5.1-14) when the frequency 
is so high that other modes in addition to the dominant are propagated. 

The equation for y corresponding to (5.1-15) may be obtained by re- 
placing tanh by coth and u by » where now 


clio 
2; CI 10 é 


w= 


aera -(5.1-16) 


Incidentally, if we set 7 = 1 in (5.1-14) and (5.1-16) and substitute in 
the expressions (5.1-4) for Lin we may show that, since #; and 2, are real, 


leo P+ lem P= 1 (5.1-17) 


Equation (5.1-17) may be obtained at once from the fact that the energy 
of the waves leaving the bend must equal the energy of the incident wave. 
It may also be shown that gio vanishes when uyCT I) = 1. 

When the above numbers are set in the three equations obtained from. 
(5.1-15) we get 


— 8.382 uw, —1.748 uw. —1.122 uz = .9492 
055 41 +2.780 uz 41.688 ug = —.2608 
| 053 u1 —3.105 ae -+7.191 ug = .1427 
from which: 
uw = — .0940, x, = .1400 + 72.0113 


The equations for 2; obtained by substituting coth for tanh are 


2354 v; — .3732 v2 -+ .3861 3 = .9492 
0717 oy + 3.1417 vo + 1.924 v3 = — .2608 
0534 2, — 3.1143 vo + 7.211 03 = = 1427 


from which 


v1 = 3.930, yr = —.1045 + 7.9803 
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When these values are set in (5.1-4) we finally obtain 


gio = .9822 + 1.1858 
gio = 0048 — 1.0255 


The following table lists values of gig obtained by the methods of this 
section*, Here the bend is in the plane of H, a/b = 2.25, \o/a = 1.400, 
where Xo is the free space wavelength of the incident wave. 1 is the 
radius of curvature of the axis of the guide. The smallest possible value of 
pi/a is 0.5. The term “approx.” refers to equations (5.1-2) while ‘1st 
order’’, ‘‘2nd order’’, etc. refers to the order of the matrices used in the 
computations. The amplitude of the reflected wave is gi and the ampli- 
tude of the wave sent forward is gio when the incident wave is of unit ampli- 
tude. 


gio 
p1/a Approx. 1st order 2nd order 3rd_ order 
6 964 +1.267 980 +1.197 .982 +7.186 
7 974 +74.224 994 +7.105 994 +7.111 
8 984 +4.178 997 +7.082 997. +7.082 
9 988 +7.153 997 +14.074 997 +1.073 
1.0 991 +7.135 .998 +7.066 998 -+7.066 
1.2 994 +7.110 998 +7.056 998 +7.056 
1.5 .996 +72.084 999 +7.043 999 +7.044 
£10 . 
6 — 1.0280 .0020 — 7.0074 .0056 — 1.0280 .0048 —7.0255 
eh —1,0068 —.0005 +7.0023 .0013 —7.0131 .0007 — 7.0066 
8 +7.0062 —.0013 +7.0074 — .0003 +7.0039 —.0004 +7.0051 
me) +7.0128 —.0014 +7.0087 —.0009 +7.0123 —.0009 +7.0123 
1.0 +7.0143 —.0010 +7.0075 —.0010 +7.0148 —.0010 +7.0147 
1.2 +1.0079 —.0002 +7.0018 —.0005 +7.0086 —.0005 -++7.0085 
1.5 —7.0040 -+.0003 —7.0034 -+.0002 —7.0041 -+.0002 —7.0042 


It appears that the values obtained from the first order matrices are quite 
far from the true values. On the other hand there is considerable agreement 
between the approximation and the second and third order values, especially 
at the larger values of p:/a. 


5.2 Bend in Plane of Electric Vector 


The calculations for this case are quite similar to those presented in Sec- 
tion 5.1. If we are to deal with the same waveguide it is necessary to 


* The computations were performed by Miss M. Darville. I am also indebted to her 
for the values given in the tables in Section 5.2 and Appendix I. 
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interchange the dimensions a and 6 so that now &/a = 2.25 and, for the 
same frequency, Ao/b = 1.400. , 

For a 90° bend the approximation (4.4-4) for the reflection coefficient 
do, for gentle bends, i.e. for € = a/py small, becomes 


do: = i€[—.0417 sin (2.23/£) + .0209 cos (2.23/8)] 


where the negligible exponential terms have been neglected just as in the 
analogue (5.1—2) for gio . 

The following table, which is similar to the table at the end of Section 
5.1, gives the results of computations for bends in the plane of the electric 
vector. 


diy , | 
p1/a Approx. 1st order 2nd order 
.6 .823 +7.547 975 +7.223 
if .887 +7.447 .994 +7.051 
8 .921 +72.380 .996 +7.042 
oO. 941 +7.332 997 +7.035 
1.0 954 +7.295 .998 +7.031 
12 970 +7.242 .999 +7.023 
Ls .982 +7.190 1.000 +7.017 
dio . 
.6 —1.0996 — .0855 +7.1284 — .0020 +7.0086 
7 —1.0848 — .0520 +7.1031 + .0050 —7.0975 
8 —i.0706 — .0330 +7.0800 + .0033 —2.0792 
9 —1.0575 — .0214 +2.0605 + .0022 —2.0635 
1.0 —4.0457 — ,0137 +7.0443 .0021 —7z.0507 
1.2 —12.0258 — .0051 +7.0204 .0007 —7.0282 
LS —12.0051 -+ .0001 —7z.0004 .0001 —2z.0062 


. The agreement between the approximation for dio and its second order 
matrix value is fairly good from pi/a = .7 onward. 


APPENDIX I 


CALCULATION OF Pym ETC. FOR CIRCULAR BEND 


It is convenient to write Ppm and Qpm as given by (1.2-10) and (1.2-15) 
in the form 


Pm = 82+ Rom, p,m = 1,2, 3, ++ 
Qtm = 82 + Lents p,m =0,1,2, 5+ 


where 6, is unity if p = mand is zero otherwise and 


(A1-1) 
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Ryn = (2/¢a) i (pip — 1)sin (rpx/a) sin (arma/a) dx 
oO. ~ % : 
, (A1-2) . 
Tom = (€p/@) i (pip — 1) cos (rpx/a) cos (rmx/a) dx 
0 
in which €9 = 1, €, = 2,p = 1,2, --- | 
In (A1~2), (1.2-11), (1.2-16), (1.4-3), (1.4-6) we make the substitutions 
p-—m=rn, u=pa/a— 2/2 = por/a 
_ptme=s, v= pp t/a+ a/2 = psr/a (A1-3) 
yeaa, w=pm/a, p=xtpi—a/2=alytu)/r 


Introduction of the integrals 


I, = (1/7) i [w(y + u)-? — 1] cos sy dy 


| 


(1/a) : (pi p72 — 1) cos (asx/a) dx 





a 7 (A1-—4) 
J, = ri ee dy = rf sin (asx/a) dx/p, 
ae af COS SY 
ads y+u 
enables us to write 
Rom = I,-— Ts, Spm = —ma (J, + J;) 
Tom = €or + :1,)/2, U om = me,a-2(J, — J,)/2 - (A1-5) 
Vom = K, — Ka, W om = €p(K, + Ks)/2 | 
where J, and K; are even functions of s and J; is an odd function of s. 
é9 = lande, = 2,p = 1,2,3,---. Since wand u depend only upon the 


ratio p:/a, the values of J;, K,; and J, depend only upon p;/a and the 
integer s. These quantities are tabulated at the end of this appendix. 
Setting y + u equal to / gives 


t= cal sin s(t — u) dt/t 


= m[Si(sv) — Si(su)] cos su — a[Ci(sv) — Ci(su)] sin su (A1-6) 
where Si and Ci denote the integral sine and cosine functions. Integrating 
by parts enables us to express J; in terms of J,. Thus . er 


[ (y + uy cos sy dy = wt— v1 cossr—amtsl, (A1-7) 
0 ; . 
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and 
I, = 2 wu — (—)'0! — wsJ,] (A1-8) 
except when s = 0 in which case . 
Io = tr wu — v1) — 1 = w/(wr) — 1 
= m/(4av) = [(2p1/a)? — i} 


- When pi/a is large, # and » are large, and the asymptotic expansion of 
(A1-6) gives 


(A1-9) 


ae a Pe sy = (2 )'po1] ats sla 3 = ° a + ares (A1~10) 
When (A1-10) is placed in (A1-8) - . 
I, ~~ or 2!ws [a ? — (—)o] — - (A1-11) 


Formulas for K; may be obtained in much the same way. 
= (p;/a) a cos s(t —2) dt/t 
u ; (Ai-12) 
= (p1/a){[Ci(sv) — Ci(sa)] cos sz + [Si(sv) — Si(sz)]sin si} | 
and when s = 0 | 
Ko = (p:1/a) log (1 + 2/2) (A1-13) 
The asymptotic expression is 
api Ks ~~ s [uw — (—)°v ] — 3's “le * — (—) oo] +e 
It is convenient to write the asymptotic expressions in terms of the new 
variable 
E=a/py (A1-14) 
- When s is even and greater than zero | 
Te~B/s, Ip ~ Oats, Ky ~ 28s? (A1-15) 
and when s is odd - 
J; ~ 2&/s, I, ~ 4éa?s~, Ky ~ lis (A1-16) 
When s = 0 . ; 
Ip~ #/4, Korw 1 + £/12 (A1-17) 


We shall need the following asymptotic expressions which may be obtained 
from the above work 


eS 6 —2 2 = 2f1 1 
ra ~ (1 - §), Su~w —@°&/2, Vn =1+é (B-m) 


where m = 2,4,6,---. 


NN PRR RB Re 


UONEWN HF OOMAD UN 
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Rin 


Sim 
Vim 


~~ 


Too a £/ 4, 


Tmo ~-8&mn-? 
Lom ~ 4&ma, 
Wom ~. 26m, 


Values of Is, Js and Ks 


pi/a Io 
5 00 


i 


(6) 


Om 5 


T 


; Uoo = 0, 
xf, Umo = 0 
Uom ~ 2a~*é 


Te 


i) 


is 


ioe] 


2.2723 2.31879 1.82979 1.43755 
1.04166 1.28232 


6 

7 

8  .64103 
9 .44643 
0 .33333 
1 = =.26041 
.2  .21008 
3.17361 
4  .14620 
5 

0 

5 


.06667 
.04167 


[NONE ON 


a 


ooooodcoooooo°oeSe 
RBBB BrP RB poh HY WH W 


12500 | 


Ji 


88256 
68200 
56052 
47844 
41905 
37361 
33780 
. 30876 
. 21803 
. 16980 


. 76232 
-44628 
30111 
. 22008 
. 16916 
13486 
11059 
.09259 
.07872 
04133 
.02564 


Je 


-99809 2.01979 2. 
21624 1. 
72356 
37231 
10615 
89624 
72581 
58448 
46507 
3628 
.0122 
8062. 


3054 1 


73339 1 
70698 
55663 
45091 
37335 
31450 
26878. 
23251 
12817 
.08128 


pre oyoil is 
. 29206 
. 18992 
13637 
10459 
.08413 
.06961 
.05926 
.05147 
.03080 
.02217 


Js 


31576 


98176 
21541 
.99327 
84343 
73508 
.65279 
. 98812 


53573 


49238 
35299 
27660 


= Ray ~ 16€mr (am? — 1)? 
Sim ~ 4a~Em(m? — 1) 
mi ™ 8Ema—?2(m? — 1)-? ~ Rin/2 
Also, if now m = 1,3,5,--:, 
Woo,= 1+ £/12 


Wao ~~ 4&0 2a? 


Ti; 

1.14772 
.37692 
.18995 
.11546 


07818. 


.05772 
.04344 
-03500 
.03168 
.02315 
£01135 
.00675 


Js 


1.48355 
.84936 
. 96683 
.41079 
31379 
. 24849 
. 20254 
. 16843 
14216 
12243 
.06596 
.04140 


ie 


94432 


. 28832 
14101 


08517 
05814 
.04298 
.03361 
.02732 
.02293 
.01969 


yo5 


1.66348 
1.04899 
77644 
62183 
.92170 
45123 
. 39869 
35788 
-92915 
29825 


343 


(A1-18) 


“(A1-19) 


Ig 


78458 
. 22090 
09986 
.05908 
03865 
02759 
02056 
01633 
01297 
01068 


Je 


15716 
-61931 
-40134 
28546 
21578 
. 16981 
13768 
11413 
.09636 
08254 


— 
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Ko Ki Ke Ks . Ki Ks Ke 
ioe) lee} ao oe} toe) ioe] ine) 
1.43874 .61659 .31832 .22547 .15539 .12199 .09275 
1.25423 .42303 .17355..11508 .06971 .05353 .03704 
1.17306 .33141 .11439 .07463 .04093 .03195 02038 
1.12748 .27575 .08261 .05431 .02737 .02214 .01325 
1.09861 .23761 .06305 .04244 .01976 .01675 .00938 
1.07891 .20953 .04994 .03471 .01509 .01341 .00705 
1.06476 .18785 .04075 .02935 .01192 .01126 .00548 
1 
1 
1 
1 
1 


> 
° — 
was 
Qa 


£05421 .17050 .03391 .02539 .00977 .00956 .00443 
.04610 .15626 .02871 .02243 .00807 .00838 .00365 
.03972 .14434 .02467 .02013 .00682 .00745 .00315 
£02165 .10508 .01332 .01330 .00358 
.01366 .08295 .00838 .01010 .00217 


NON RwWNHOCAUID 


Oe NO 


APPENDIX II 


Functions or Atmost DiaconALt MATRICES 


Let £ be a matrix whose elements are small in comparison with unity. 
It is then often possible to approximate a matrix defined as some function of 
the matrix J + £, where J is the unit matrix, - the expansion 


fd +B) = 1) +3 af +5 Ea) +: (A2-1) 


Thus, for example, when we take to be 2~! we obtain 
(ea ee ee ee 3 (A2~2) 


Here we shall give similar formal results for f(D + £) where now Disa 


diagonal matrix 
E QO - O 
D= 0 t : ; 


0 0 0 wy 


(A2~3) 


whose diagonal elements are unequal and the elements £;; and £;; are small 
in comparison with the absolute value of |d; — d;|. We shall restrict 
ourselves to a first approximation of the non-diagonal terms of f(D + £) 
and to a second approximation of the diagonalterms. The results are closely 
related to the ones obtained from the upeceusDabOn theory used in wave 
mechanics. 

We assume that f(D + £) may be ee by the series 


f(D + B) = al +a(D+ 4) +a(D4+ B?+-+-  (A2-4) 
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where a@,, is a scalar and 
(D+ B= (D+ E(D+ 5H) =D'+DE+ED+E 


and so on. The sum of the terms independent of F is f(D). The terms 
_of order £ are 


EinD+E 

DE+ ED in (D+ 5)’ 

DE + DED + ED? in (D+ E)° 
> D'ED" in (D + E)* 


(A2-5) 


where the summation extends over the non-negative integer values of ¢ 
and m for which + m= 2-—1. Theelement in the 7th row and th column 
of D'ED” is d;E;;d; and hence the corresponding element in the summa- 
tion in (A2-5) is 


Ei; Ld; d; = 


j (d? —dj)/(d;-—d), ix j 
= not ; . (A2-6) 
n a; t=]. 


Thus the terms of order E in the ith row and jth column of i + £) are, 
from Ne 6) and (A2-4), 
d = ty (A2-7) 
Exf' (di), t=] 
where the prime on f denotes its first derivative. 
The terms of order E’ in (D + E)” are 


2 D*' ED‘ ED” = 2 (d® Bedlald; Be Gla (A2-8) 
ktm k.tym 


where the summations extend over all the non-negative integer values of 
k, £, m for which & +£+m=n-— 2. On the right [d:E;;]n denotes a 
square matrix whose element in the ith row and jth column is d;E;; . 
Likewise the second factor in brackets is a matrix having d;£;;d; in the ith 
row andjth column. The element in the ith row and jth column of (A2-8). 
is, from the rule for the product of two matrices, 


> > (di Ei.) (ds E40?) = YBa By >» dt af 
kfm S=1 


Iii, s, andj are unequal the sum in k, £, m is 


1 [—4 - a a 
d;—d;_d;—d, d;—4d, 
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and in case of equality the sum may be found by a limiting process. Since 
we are interested in this order of approximation only for the diagonal terms 
we set 7 = j and obtain for the sum. 


ndj sd — ds 
d; pam ds (d; =: d,)?’ 
n(n — 1) 
21 


_ Thus the contribution to the ith elggonals element of f(D + E) from terms 
of type (A2- 8) is 


ae ay fd) — f(di) — f(ds) 
ay (d;) > Eis Ess E a = (d; = d,)? | (A2-9) 


where the prime on Z-indicates that the term s = iis to be omitted. 

Thus, to summarize, we may say that the first approximation to the non- 
diagonal term in the 7th row and jth column (i ¥ 7) of f(D + £) is 

fd) — f(G) . 
Ey; A2-10 
or ne ( ) 

and the second approximation to the diagonal term in the ith row and ith 
column of f(D + £) is 


iAsS 


—2 ‘ 
d;-, t= 5. 











f(a) + Exf'@ + Bi pay) 
, i'd) _ f(di) — fds) 
+d BuBal 7 - (d; “(d; — d,)? (a) = 14) 


where the primes on f denote derivatives and the prime on 2 indicates that 
the term s = 7 is to be omitted. 

Two results obtained from (A2-10) and (A2-11) are of interest. Tor 
the first result we set f(z). = 2-1 and get the following approximations to the 
elements of (D + £)7?: 


= Baldga)-, 1#q 
N 
de = E = 2d ase i= j. 


For the second result we set f(z) = 2”? and obtain the following approxi- 
mations to the elements of (D + BE)"; 


| Bulay 4 i), i#j 
: | 
de hag E = dX Ei. Egi(di? + anny? | Mey: 


In (A2-12) and (A2-13) the’ summations include the term s = i. 


(A2~11) 





(A2-12) 


(A2-13) 
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We shall now state several results related to the above formulas. Let u 
denote the matrix D + E so.that the typical element w7;.= Ei;, 71 ¥ 7, 
and #;;= d;-+ £;;. Then the latent roots A; , A: , --: Aw of w are the roots 
of the equation obtained by setting the determinant of AJ — w to zero: 

A Uy —u12 
[ALJ —u|=|--t1. N— wy +| = 0. (A2-14) 


The modal column &; corresponding to the jth root \; satisfies the matrix 

equation 

At — u)k; = 0- (A2-15) 
Since the non-diagonal elements of « are small, we see from (A2-14) that 
we may label the roots so as to make ),; nearly equal to w,;; , and this together 
with (A2-15) shows that all the elements of k; are nearly zero except the jth 
which we may choose to be unity. When these approximate values are 
taken as a first approximation in the process of solving (A2-15) by successive 
approximations, the second approximation is found to be 

N N 

Ay = yj + pa SEO Uys ay Ujs Re; 


sal Ujz — Usgs 


Reo; (A2-16) 


Us; ° 
; ar ee, 
Use — Uss 


a 


where the 1 in the column for &; occurs as the jth element. This expression 
for \; occurs in the perturbation method often used in wave mechanics. 
For the modal row x; corresponding to A; we have in much the same way 


a 
*~, 
lI 
oa 
3 
| 


(Xf — u) = 0 
kz = [Kaj , Kags ct Kaz, * 7+ Kwa] (A2-17) 
_ _ Mis Ki | | 
kg = ———— 
tj; — Use 


where the last expression is an approximation and where x;; may be chosen 


at our convenience. 
The results (A2-10) and (A2-11) may also be obtained. from (A2-2), 
(A2-16) and the relation* 


fs) .0 + 0 
flu) =k if fe) i kr (A2-18) 
0 0  f(dy) 


* This is equation (12) in Section 3.6 of Reference’. Although proved only for poly- 
nomials it may be verified to be true for the applications which we shall make. 
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where & is the square matrix [ki , k; , --» Rw] composed of the modal columns 
k; , and the latent roots are assumed to be distinct. The element in the ith 
row and 7 column of k&-? is, to our order of approximation, equal to —k,; 
when 7 # 7 and to 


N 
1+ > Ris Re (A2~19) 


when i = 7. Here, as usual, the prime on 2 indicates that the term s = 7 
is to be omitted. . 


Another result which is sometimes useful is the approximate solution of 
UyX1 + texte + ees H iphyp Fees tntn = 1 


UaiX1 + Urge + +++ + Usptp ++ + Uewtwy = m (A2-20) 


ee 


xj = i _ Sete Gp — (A2-21) 


N 

Np 7 Ups 
= 3 (np Usp — Ns Upp). 
Upp s=1 Uss Upp 





The prime on ~& indicates that the term s = p (which happens to be zero 
here) is to be omitted. A somewhat more general result which states that 
the ith element in the column matrix (D + E)~‘7 is 


N 
nd, + yD E;, dy ds"(Esinis — dis) (A2-22) 


may be derived from (A2-12). Here the elements of the column matrix 
need not be restricted by | 7;| «| 1p |,7 # p- 
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- The Approximate Solution of Linear Differential Equations 
By MARION C. GRAY and S. A. SCHELKUNOFF 


Linear differential equations with variable coefficients occur in many fields of 
applied mathematics: in the theories of acoustics, elastic waves, electromagnetic 
waves in stratified media, nonuniform transmission lines, wave guides, antennas, 
wave mechanics. The ‘Wave Perturbation” method described in greater detail 
elsewhere! is particularly useful in those ranges of the independent variable in 
which the “WKB Approximation” is not sufficiently accurate. The present 
paper endeavors to illustrate the remarkable accuracy of this method, particu- 
larly when compared with Picard’s method. 


I. INTRODUCTION | 


N A recent paper! the approximate solution of linear differential equations 
I by a wave perturbation method was described. When the method was 
applied to equations whose exact solutions were known we were greatly 
impressed by the rapidity of convergence of the successive approximations. 
Hence the purpose of this note is to present some illustrations in the hope 
that others may be interested and may find the proposed method an im- 
provement on those now in use. 

In essence the wave perturbation method dates back to Liouville’, but 
in his mémoires he was interested in a problem of heat conduction involving 
a non-homogeneous differential equation with homogeneous boundary 
conditions, whereas we consider a homogeneous equation 


y= F(x)y — 
with non-homogeneous initial conditions 
y(a) = 1, y’(a) = 0 a (2a) 
or , 
wa)=O0 ya) = 1, (2b) 


the solution being desired in an interval a S$ x S b. Since the solution 
for any assigned initial or boundary conditions can be expressed as a linear 
combination of the solutions satisfying (2a) and (2b) we have not imposed 
any real limitation. 


II. THEORY 


Comparison of the wave perturbation method with Picard’s method 
(which is essentially a linear perturbation method) is particularly instruc- 
» 350 . 


LINEAR DIFFERENTIAL EQUATIONS 351 


tive. It will be recalled that in Picard’s formulation the differential equa- 
tion (1) is replaced by an integral equation 


y(x) = y(a) + @& — a)y'(a) + [ : F(u)y(u)(x — u) du (3) 


where y(a) and y’(a) are assigned initial values*. Writing 


L(x) = y(a) + (w — a)y’(a), (4) 
L,(x) = [ F@LAwe — u) du, n=1,2,3,---, 

the series 
yx) = Lo(x) + Li(w) + L(x) + --- (5) 


is shown to converge to a:solution of the original equation. In practical 
applications, unfortunately, it is usually found that the successive approxi- 
mations converge rather slowly unless the interval (a, 5) is small. 

In the wave perturbation method we first rewrite equation (1) in the 
form . 


oy" = — By + [P+ Fly = —By + fa)y, (6) 


and instead of the’ integral equation (3) we use 


v(x) = (a) cos Ble — a) + 5 ¥"(a) sin Bw — a) 
Lp (7) 
+ 3 [ f(u)y(u) sin B(x — u) du. 


The parameter 8 is arbitrary and might be defined in various ways. We 
have found it convenient to use the definition 


f= i F(a) de, (3) 


/ 


so that if F(«) is negative Bis real and our first approximation 
1 : 
W(x) = y(a) cos B(x — a) + 3 y’(a) sin B(x — a) (9) 


is sinusoidal. If F is positive 8 is imaginary and we start with an exponen- 
tial approximation. If F changes sign in (a, 6) the best procedure is to 


* This is not quite the usual form of the integral equation but it is substantially 
equivalent. 
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subdivide the interval and obtain separate approximations, though this is 
not necessary if F is predominantly of one sign throughout (a, 6). To (9) 
we now add the sequence 


W,(«) = ; i ‘ f(a) Wp-1(at) sin B(x — 1) in, | (10) 


and the eee 
y(a) = Wo(x) + Wilw) + We(x) + --- (11) 


is the desired solution. 

The flexibility of the wave perturbation method as compared with Picard’s 
linear method lies essentially in the introduction of the variable parameter 6. 
Since we make @ depend on the length of the interval (a, 6) in which a solu- 
tion is desired the approximations may be extended over much longer 
intervals than is feasible in Picard’s method. If F(x) is a slowly varying 
function throughout (a, 0), so that f(x) is small, it will be found that the 
first approximation IV)(x) is good, and the second, Wy + 1’; is > generally 
adequate. 

Another choice for @ is 


(12) 





However, the integration in (8) will often be simpler than in (12). 
Picard’s method is a special case of the wave perturbation method, with 
8=0. In fact, if P(x) changes sign in (a, b), then in some cases B as defined 
by (8) will reduce to zero. 
If /*(x) is a rapidly varying function, or if the solution is desired over an 
infinite interval, it is usually advantageous to transform equation (1) by 
first introducing a new independent variable 


= [ / —F(x) de, (13) 


and then removing the first order term in the new equation by an appro- 
priate transformation of the dependent variable. 
Tit. EXAMPLes 
For our illustrations we have used mainly the simple equation 
Dp ey _ (14) 


whose exact solution can be expressed in terms of Bessel functions of order 
+ 1/3. Since the Bessel functions are oscillatory in nature it might be 
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suggested that comparison with Picard’s method is weighted in our favor. 
This does not seem to be the case, as will be illustrated in example 4 where 
the exact solution is a monotonically increasing function. It has also been 
suggested that Picard’s solution might be improved by starting from a better 
initial approximation, say JV), rather than from the linear approximation 
Lo, but we have not found any marked improvement in the succeeding . 
approximations (see examples 1 and 2). The various points of interest 
will be brought out in our examples, with the accompanying figures, which 
we shall now briefly describe. In each figure the heavy curve is the ac- 
curate solution while the approximations are indicated by self-explanatary 
letters. 

Example 1, Fig. 1 


y’ = —-xy,0S%82 
y(0) = 1, y'(0) = 0 
Exact solution: y(x) = 1(2)37 8 oll? J_y (37?) 
(a) Wave perturbation 
Wo = cos x 
W, = —ixcosx+ 4(1 + 2x — 2’) si 


(b) Linear perturbation 


lI 
— 


3 
x 
L=-= 


6 


é 


(c) Linear perturbation using initial sinusoidal approximation 
Ino = cosx = Wo’ 
I, =x+«cosx— 2sins 

Example 2, Figs. 2, 3 and 4 


YY 


yo = —xy,2 Sx%56 


y(2) 


0, (2) = 1 
Exact solution: 


y(w) = —.84423a"? J_1j9(a°!) — 0192914"? Fy)3(3x"”) 
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(a) Wave perturbation, Fig. 2 
Wo = #sin xe — 2) 


12 — 8x + x” 
W, —" 6 
Figure 2 exhibits rapid pulling of the successive approximate waves to- 
ward the exact even though the interval has been chosen deliberately 
unfavorable to the straight wave perturbation method [see example (c) 
and Fig. 4 for the improved treatment]. 
(b) Linear perturbation, Fig. 3 


Lye 2 
as (16 — 16x + 4x3 — x4) 


16 lox 207 x of x" 
Lo = —=—_ oO '"——, 
‘ 6 9 7 9 a0 504 


Using IV instead of Lo 





2 sin 2(a@ — 2) + cos 2(% — 2). 


Sy 
| 


I, = 3- 5 + gsin Ue — 2) +3 cos 2(a — 2). 


(c) Preliminary transformation of variables, Fig. 4 
Introduce 6 = 3x°”, y = o* 


and the modified equation is 


4/2 Y 
ey <p = 
gy! = -(1+35)0 3 S65 476. 


Then, using for simplicity B = 
vo = 2160/8 sin (0 — 80), 00 = 40/2/2 


or 
Wo = (2x) sin 2 (8/2 _ 98/2). 


It will be seen that Wo is a very good approximation throughout the 
range (2,6). Adding W, obtained from 


ni = 58 fcos(@ + 09)(Si 20 — Si 26) — sin(@ + 6)(Ci20 — Ci 2H] 


36/2 
the accurate curve y is reproduced. 

In Fig. 2 the third approximation could not be disticuisied from the 
accurate curve though numerically the values are not identical. For 
purposes of comparison the table of numerical values (Table A) may be 
found interesting. . 





BELL SYSTEM TECHNICAL JOURNAL 








LINEAR DIFFERENTIAL EQUATIONS. 359. 




















TaBLe A 
: 1 2 3 
x y Wo SV SIV Pai 
) 0 0 

2.0 Oe. 0. 0. 0. 0. 

2.2 19721 .19471 . 19720 19721 .19721 
2.4 37694 .35868 .37668 .37694 37694 
2.6 .52056 | .46602 51885 52053 .52056 
2.8 .61035 49979 .60442 _ 61020 .61034 
3.0 .63236 | 45465 .61792 .63180 63232 
3.2 .57922 33773 .55169 57781 .57918 
3.4 45287 . 16749 -40907 44995 45726 
3.6 .26584 — .02919 .20603 26085 .26561 
3.8 .04126. — .22126 — .02974 .03408 ~ .04087 
4.0 — .18921 — .37840 — .26229 — .19800 — .18974 
4.2 — .38951 — .47580 — .45326 — .39852 — .39011 
4.4 — .52506 » — .49808 — .56889 — .53240 — .52557 
4.6 — .56943 — .44173 — .58697 — .57328 — .56964 
4.8 — .51062 — .31563 — .50217 — .51000 — .51044 
5.0 — .35548 — .13971 — .32847 — .35076 — .35494 
5,2 — .13068 _ 05827 — .09772 — .12364 — .13006 
5.4 — ,12052 .24706 14547 12725 . 12080 
5.6 34582 39683 .35200 34988 .34536 
5.8 49485 .48396 47807. 49525 .49347 
6.0 .53114 .49467 .49467 52903 52903 

Example 3, Fig. 5 
2 
oe as 7 
yeryr sy ters” 


y(1) = 1, y’(1) = 0 


Exact solution: y(x) = sin (w — 1) + = C08 (7 — 1) 
(a) Wave perturbation, with the initial conditions satisfied exactly: 
Wo = cos (x — 1) _ 
WW, = 2sin (x — 1) — 2 cos (x + 1) (Ci 2x — Ci 2) 
— 2sin (w + 1) (Si 2x — Si 2) 
(b) Wave perturbation, matching the exact solution at infinity 
Wo = sin (x — 1) 
Wi = 2 sin (x + 1) Ci 2x — 2 cos (x + 1)(Si 2a — 2/2) 


This is an example of a solution in an infinite interval, where the per- 
turbation term is not small throughout. It is interesting to note that the 
second form gives good agreement with the accurate solution in most of 
the range of integration. 
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Example 4, Fig. 6 


IIA 
bo 


yi=+ay, O0S2 


y(0) = 1, y’(0) = 0 
Exact solution: 


ye) = 1R)3 TP  Tas(Ba”). 
(a) Wave perturbation 
Wo = cosh x 


W, = —$ cosh a + (a2 — 2x + 1) sinh x 


(b) Linear perturbation 


Io = 1 
3 
Ly = =. 


This is an example in which the exact solution is non-oscillatory yet even 
in the short interval (0, 2) Wo + W, is a better approximation than Ly + 
Dh. sar 
Example 5, Table I 


W tov +y=0, 0<«S« 


Solution required to match the accurate solution 


w(x) = Jo(x) — 7 No(x) 
at infinity: 














. ther 2 
V= = «* 
Tx 
= 1 1. 1. 
W,= pt? ie! Cig = 3 Sidx — 7) I, 
An/ rx 2 
TABLE I 

x Jo 3Ny Wo Wot: WotWitWe 
10 — .2459 —7z .0557 | —.2468 —i .0526 | — .2460 —:z .0557 
9 — 0903 —z .2499 | — .0938 —7z .2489 | — .0903 —i .2500 
8 £1717 —i .2235 .1683 —i .2264 1717 —i .2235 
7 .3001 +7 .0259 .3009 +i .0207 .3001 +7 .0260 
6 .1506 +72 .2882 .1568 +2 .2855 .1507 +7 .2883 
5 — .1776 +7 .3085 | —.1704 47 ,3135 | —.1777 +7 .3086 
4 — .3975 +i .0169 | — .3979 +7 .0291 | — .3973 +7 .0169 
3 — .2601 —i .3769 | — .2765 —i .3684 | — .2601 —i .3772 
2 .2239 —i .5104 .1967 —7i .5288 .2246 —7i .5109 | © 
1 7652 —i .0883 .7796 —i .1699 .7683 —2i .0860 .7051 —é .0882 
0.8 .8463 +7 .0868 .8920 —i .0130 8499 +7 .0916 .8461 +7 .0868 
-0.6 .9120 +2 .3086 ; . 1.0124 +7 .1899 .9152 +72 .3183 .9116 +7 .3084 








LINEAR DIFFERENTIAL EQUATIONS 363 


For values.of « less than 1 We was evaluated numerically. 
Example 6, Table IT 


yt ty = (5 -1)s, eee oe 
y(1) = 1, y'(1) = 0. 


The solution of this equation using Picard’s method and the integraph 
has been described by Thornton C. Fry. We compare his results with those 
obtained by the wave perturbation method. The equation is first reduced 
to normal form by the substitution y = « ”” w, so that 


3 
| un = (-1 +3)x 


and we have 8 = 44/3. Then 


xv? Wo = cos B (x — 1) + esin B(a — 1) 


oe Wy, = a (3 — 1) cos B(u — 1) 


1. ; 
+ a B(u — 1)| sin B(x — i) du. 


While JV, may be evaluated in terms of Cz and Si functions the values 
tabulated below were obtained by numerical integration. The values of 
the accurate solution 


y = 1.4034 Ji(x) — 0.3251 Vi(x), 


and of the third and eighth Picard approximations, are copied from Fry’s 
paper. 




















TABLE II 
x y 93 ys Wo WotWy 
1.0 1.000 1.000 1.000 1.000 1.000 
1.2 - 998 .998 .998 .990 .998 
1.4 985 984 .986 | .961 985 
1.6 .956 951 955 913 .956 
1.8 .908 894 .910 848 | .908 
2.0 842 809 844 .769 842 
2.2 759 694 .760 | .677 758 
2.4 .659 548 .661 | .575 659 
2.6 547 .370 - 549 .466 547 
2.8 425 .156 427 352 427 
3.0 297 —.096 300 236 300 
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ea Coefficients for Ground. Return Circuits 


By W. HOWARD WISE 


This paper is concerned.with the effect of the finite conductivity and dielectric. 


-constant of the earth on the potential coefficient for a 1-wire ground return 
circuit. It has been customary to say that the potential coefficient V/Q is 


pr = 22 log p”/p" , elm units per cm. 


It is generally realized of course that this is just a good approximation to the 
true pio. To see that it is just an approximation one has only-to-imagine the 
earth turning into air, in which case the distance p” will eventually cease to 
have significance. The object of this paper is to derive the complete expression 


for pie. 
It turns out that 


12 = ¢[2 log p”/p’ + 4(M + iN)] 
where 
© eM EVAL cos y Ey/ at 
) VE ei + (€ — i20r0)t™ 
a = 4row, as in Carson’s work on Zj2 and in mine on Z12 at high frequencies!? 
ge) = 4/1 + ie — 1)/2cd0 = s 


dielectric constant in electrostatic units 
conductivity in electromagnetic units 


M-+iN = 


a 


Coz 
= 10-8 to 10- in ordinary soil 
\ = wavelength in centimeters 
c = velocity of light, in cm per sec. 


M + iN vanishes asf—0,f— ©,e— © org > ~, 


Ordinarily 4(Mf + iN) will not be an important correction to 2 log p”/ p’; but if the 
frequency is high and /t or 2 is small it can be a worthwhile correction. For 
example, if a .02535 inch wire be thrown out on the ground to be a 2 mc antenna 
and we assume thato = 107", = 15 and # = 3 cm. then, with a for wire radius, 


pu = [2 log 2h/a + 4(M + iN)] 
= (10.455 + .152 + 7.319}. 


1/p11 is the capacity to ground. If there were two parallel wires the scalar 
potential at the first wire would be Vi = (1101 ot PO. 


DERIVATION OF THE FORMULA 


(7) 


(8) 


E BEGIN with the wave-function for an exponentially propagated 
current in a straight wire parallel to a flat earth. The wave-function 


for a horizontal current-element dipole has been formulated as an infinite 


1 John R. Carson: “Wave Propagation in Overhead Wires with Ground Return,” 


Bell Sys. Tech. Jour. 5, pp. 539-554, 1926. 
2W. Howard Wise: “Propagation of ee Currents: in Ground Return 


Circuits,” Proc. I. R. E. 22, pp. 522-527, 1934. 
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integral by H. von Hoerschelmann’®. The wave-function for the current 
in the wire is obtained by integrating the wave-functions of the current- 
element dipoles along the wire from minus infinity to plus infinity. It is 


[ve] oe tFR1 gia 
T=axe” i re ( a 
— 








Ri Re 


” 2F (vp) et | ; | 
+) eee yedr) de + 5X0 (1) 


= —yx ic yx oO (1 _ 7) Jolvp)v ei “an 
c X 2e ee al imeem) dv-dx. 








Fig. 1 


The time factor ise’. a; b and c are unit vectors pointing in the x, y 
and zg directions. 


1/2 


R= + 9+ (h- 9 
Ro = G+ + (+ 2)” 
p = (+ y)? 

ik = Qn/d 


2 5 é 4 . 
ky = enw? — i4ropw in electromagnetic units 


By supposing ¢ to be measured in electrostatic units we can write 
= Re — i2cdo)y. 


3H. Von Hoerschelmann, Jahrb. der draht. Teleg. 5, pp. 14-188, 1912. - 
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It is assumed that uw is everywhere unity in electromagnetic units. 


L=(P— RB)? m = (2 — B35)” 
w=h-+ sz, 7 = R/ke | 
y = a-+ iis the desired propagation constant 


The electric field parallel to the wire is 


ie, Al oe er 
E; = ia +k 2(2 II, +5 yt & tt.) 


her aV o 
= —-Ihe™ Zy — a 
It has previously been shown that? 
Zy, = t[2 log p’/p’ + 4(Q — 7P)] (3) 
where 
0 = 4P = (fF as + is? — ve” cos y'v- dy, 


al 
w’ = wVaandy’ = y Va. 


To get the potential coefficient for a ground return circuit it is necessary 
to compute the scalar potential. 


V = iwk” (2 I, + ey + 2 I, 
ox oy 02 (4) 


= 0 pir 


_ As in previous work the propagation constant y is assigned the value ik 
as a first approximation. This is an ideal value for y but the following 
considerations make it an imperative choice: (1) to assume that the current 
is propagated down the line with a velocity less than that of light makes the 
integrals very hard to evaluate, (2) to assume that the attenuation is not 
zero on an infinite line amounts to assuming an infinite source of energy and 
makes the integrals diverge. 

It should not be inferred that the resulting formulas are necessarily poor 
if the physical system does not closely approximate the ideal one in which 
yistk. ikis employed as a convenient first approximation in evaluating the 
correction terms in Zy and py. Eventually, if there were but one wire, 
one would compute y = ~/(z + Zu)(G + iw/pu), wherein Z1, and py 
have been evaluated with 7k for y, and this would be a second approxima- 
tion toy. Past experience with the second approximation so obtained has 
justified the expectation that it would bea satisfactory final result. Since 
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the integrals diverge if the attenuation is not zero the use of an infinite line 
formula presupposes reasonably efficient transmission. 


‘Since IL « e7” we have 





oe 0 _ : ee I _ cl 
1wk 5y Ue = z us = 7B o8.= 7 21: 
; Of. 00% fay Pe ae 
Since ae ap ikI = iwQ or I = 0 lus is 
jee . We = Qc? [2 log p/p’ + 4(0 — iP). (5) 
We have next to consider 
ia, 0 _ 2w(l = ae o oo 8 7) ar Jo(vp)e”' v+dy 
k? dz ik? (1 + m)(l + 72m) 
- sma )0 f ey. dp [ie 
7 Qe ae) ome (i+ m)(l + 7m) ax ze) om 


The infinite integral is 


ro) eo ye dy a 
0 + ml + Pm) fe * Jo(vp) 


Since Jo(v ~/ x2 ++ 92) and cos kx are even functions of « and sin kx is an 
odd function of x 





; tel ~Y To(vp) ie]. 


top 
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JolvV/ x22 + ye dx = 2 Tov -V/ 22 + 2) cos hu-dx 
L— 00 8 y 
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O if v<k 


cos y4/p2 — kk? , 
2 ————.. if k< 
VP RB. ae 


and so our integral is 


—wl 


ev cos yl 


rik | (1 + m)(l + 72m) vray on 





or, since 2? = v* — PR, 
a e~” cos yl-dl | 
4 = . 2 
0 C+ VE + PR — B))I+ OVE + PE — B)) 
or, if we put J.= vila aa | 
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and i(k} — k*) = 4row (1 +: it) = as’, 
2k f*  €”’* cos yl vedy 
Vah C+ VP +i\y te Ve + is?) 
where, as in Q — iP, w! = w/a and y’ = y Va. 





. 0 —wily —" wily 
Noting next that or e = —vV/ave 
we have 
tw 0 Ca ssa 





Boz = -004 [ ENT + is)(v+7 2) + is?) 


2 — 
= =-o¢ 4 LE eae — ¢””* cos yy (1 — 7°)vedy. 
“40 vp T v 1S 


Since (1 — oe =v+r Vii + ist Lis? — 72(4/o? + Gs? + 7) 
this 1 1s. 


2 , 27 2 
2 9 ° T* tS —w'ly ? 
— Oc" > [ Vr bis - vy - 9. 77> = le cos y’ v-dv 
ofa [| + pt ore + is? 4 








we _ e-”’” cos y’p- dy 
=o - 4(Q inn+a ly Sere etl (6) 
On adding (6) to (5) we have 
Opry = Qe[2 log p/p’ + 4(M + iN)| | (7) 
where 7 , 
VV 2  e-”'® cos yp 
a = eee 
7 -[ MAAS ACF am 


VB + ie + t/7 


M + 1N vanishes asf > 0, f—- ©, €— © oro > w, 


When #; — &” is minute the leading terms in the approximation (9) 


for M + iN are= — ae log (p” +/(ke — k*)/2) — i . 
8 2 ws 4 
AN APPROXIMATION FOR M + iN 


It is possible to get series expansions for M -++ 7 but those which have 
been obtained do not facilitate computation. A fairly good approximation 
to M + iN is arrived at as follows. 
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-i27 2 _ t(nta/4) 
Let ie" = uv, “u=e 
¢€— 12ck0 = a= 1/7? 
e Gta)iVet cos yo~/al = He" ae Oe a 


= (h+s-iveVo 
a (h+2e+ iy) Vo 
g =(h4t+2f Va 
Since (f° + ey? = (P+ ate — 2iu)"” 
=ttu—tu/titut+.. 
we put 
1/(fe + 2 a at) = 1/t@+ u—tu/(t+ u)+ al) . 
=(¢+a)/(a+ DP + (t+ 1jtut wv] 
= (t+ n+ 7)/(a + 1) + 7)¢ + 12) 
where r, = u(1 — 0/1 — 4/(a + 1))/2. 
and or = ui + V1 — 4/(a + 1))/2. 
Then 








M+iN®& dt 








1 C( rm ——— 
(a+ 1)\(2e-n)Jo \VEtrn t+r es: 


1 
~ 2a + 1)(r2 — 1) 


— ne (eo?) + reli (e?"')] (9) 


[—ne* (Ee? + rye? li (e9'"2) 





When y = 0 this reduces to 


J 
(a + 1)(r2 — 11) 





M+inNw [—ree™ Li (e") + rye’? li (e€°)]. (10) 





* iad —t bead ieee: 
lie“) = -| “dt =C + logs + 2 i 


where C = .577215665 and, if 2 =r". —r< 9 <a 


2 1=0 (—z)! 


li(e*) ~ yar of Ae Ras 


The accompanying charts give M and N for y = Oande = 15. With 
z = h they give M and N for 1. The computed points are indicated by — 
solid dots on the chart for M; they were obtained by numerical integration. 

The approximation (10) was checked against the values obtained by 
numerical integration at a number of points. The discrepancy in each case 
amounted to less than one per cent for both Mand NV. This approximation 
is a much easier way to evaluate the integral than is numerical integration 
but it is a tedious computation with many chances for error. Conse- 
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Fig. 3—N, computed with e = 15 and y = 0. The number associated with a curve isthe value of (k + s)¢Va. 
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‘Fig. 2—M, computed with e = 15and y= 0. The number associated with a curve is the value of ( + z) tva. 
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quently it is important to observe that a coarser kind of approximation 
may often be good enough. Thus, ore y to be zero, 


i} uid te —guia 
dt = —- &!* Hi (e-™*), 
a 





utive [ + at 


If g is very small one might use 


utive [#4 +f aot 


1 a i pee = 
= Hog (1.4 4) Tago (e ie 


| Ordinarily precision is not required in M + iN because 4(M + oN) is a 
small term in pie. 





Abstracts of Technical Articles by Bell System Authors 


Electrochemical Factors in Underground Corrosion of Lead Cable Sheath. 
V. J. ALBANO. Stray current is the principle cause of corrosion failures on 
underground telephone cables in most cities where trolleys are operated. 
To mitigate this condition, the cable sheaths are “drained” to the negative 
return system of the traction system. By this means not only is the stray 
current anodic area largely eliminated, but the cables automatically become 
negative to earth, and therefore are cathodically protected. The cathodic 
protection afforded in this manner prevents other types of corrosion from 
occurring. With the gradual abandonment of trolley systems, and with the 
extension of underground cables into non-trolley areas, the percentage of 
underground telephone plant receiving this protection is decreasing. As a 
result, the problems of lead corrosion due to such causes as galvanic and local 
cell action of various types, and chemical action by substances in the soil 
are becoming more prevalent. It is the purpose of this article to review some 
of the basic principles of corrosion not involving stray currents, and show 
how they apply to the problems of lead cable sheath corrosion. 

PCM Equipment. H. S. Biack and J. O. Epson. PCM, pulse code 
modulation, is a new solution to the problem of overcrowded frequency 
spectrum. It appears to have exceptional possibilities from the standpoint 
of freedom from interference, and seems to have inherent advantages over 
other types of multiplexing. . 

Coaxial-Cable Networks FRANK A. Cowan. This paper discusses the 
general features of the coaxial system, its application for both telephone and 
television, and the future prospects for very-broadband transmission facili- | 
ties in the communication network. 

Parabolic-Antenna Design for Microwaves! C. C. CUTLER. This paper 
is intended to give fundamental relations and design criteria for parabolic 
radiators at microwave frequencies (i.e., wavelengths between 1 and 10 cen- 
timeters). The first part of the paper discusses the properties of the para- 
bola which make it useful as a directional antenna, and the relation of phase 
polarization and amplitude of primary illumination to the over-all radiation 
characteristics. In the second part, the characteristics of practical feed 
systems for parabolic antennas are discussed. 

1 Corrosion, October 1947. 

2 Electrical Engineering, November 1947. 


3 Proc. I. R. E—Waves and Electrons Section, November 1947. 
4 Proc. I. R. E., November 1947. 
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Microwave Antenna Measurements» C.C, Cutter, A. P. Kinc and W. E. 
Kock. A description is given of the techniques involved in measuring the 
properties of microwave antennas. The measuring methods which are 
peculiar to these frequencies are discussed, and include the measurement of 
gain, beam width, minor lobes, wide-angle radiation, mutual coupling be- 
tween antennas, phase, and polarization. The requirements of the antenna 
testing site are taken up, and components of a complete measuring system 
are briefly described. oo 

Microwave Converters.6 C.F. Epwarps. Microwave converters. using 
point-contact silicon rectifiers as the nonlinear element are discussed, with 
particular emphasis on the design of the networks connecting the rectifier 
to the input and output terminals. Several converters which have been 
developed during.recent years for use at wavelengths between 1 and 30 centi- 
meters are described, and some of the effects of the impedance-versus-fre- 
quency characteristics of the networks on the converter performance ‘are | 
discussed. . 

Recent Developments in Relays :" ClasePpclased Reed Relay, W.B. ELLWwoop; 
Mercury Contact Relays, J. T. L. Brown and C. E. PoLtarp, Relays which 
combine high-speed and great uniformity of performance over long periods 
of time are required for some uses in the telephone plant. The relays de- 
scribed possess these qualities to an unusual degree. Detailed description 
is limited to two types, each typical of a generic ces in which the prin- 
ciples involved apply to all. 

These relays are based on the philosophy that a motor element (any dévics 
for conversion of electromagnetic to mechanical energy), which is efficient 
and magnetically and elastically stable and operates contacts sealed in a 
proper atmosphere free from dirt and film, will give reliable performance if 
the contact load is engineered to the capacity of the contact.. The relays 
require no maintenance beyond unit replacement, for there is no possibility 
of a change in adjustment after assembly is completed. 

In one form the contact is provided for by metal in solid form, while in the 
other a mercury film supported on solid metal surfaces provides the con- 
tacting medium. The mercury at the contacting surfaces is replenished 
continuously through a capillary path from a mercury reservoir below the 
contact. 

An Adjustable Wave-Guide Phase Changer’ A. GARDNER Fox, A very 
interesting and useful component of the wave-guide art is the differential _ 
phase-shift section, wherein dominant waves of one polarization are caused 
to. travel through a section of wave guide at a different velocity than waves 


5 Proc. I, R. E., December 1947. 
6 Proc. I. R. E., November 1947. 
? Elec. Engg., November 1947, 

8 Proc. I. R. E., December 1947. 
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polarized at right angles to the first. Particularly useful are the A90-degree 
and A180-degree differential phase-shift sections which produce differential 
delays between the two polarizations of 90 degrees and 180 degrees, respec- 
tively. The properties of these sections are discussed, and it is shown how 
they may be combined to form a phase changer which will transmit substan- 
tially 100 per cent of the incident power with a phase which is readily ad- 
justable. Several different methods of building these sections are finally 
described. | 

Considerations in the Design of a Radar Intermediate-Frequency Amplifier.? 
ANDREW L. Hopper and STEWART E. MILLER. The intermediate-frequency 
amplifier of a microwave radar receiver is commonly required to provide 
approximately 100 decibels amplification in a bandwidth of 1 to 10 mega- 
cycles, centered at frequencies in the 30- and 60-megacycle regions. Meet- 
ing such requirements involves the use of five to ten amplifier stages of the 
highest efficiency that can be suited to production methods. In addition, 
the noise figure of the radar intermediate-frequency amplifier is a significant 
contributor to the over-all radar receiver noise figure, and must therefore be 
maintained at an absolute minimum. By examining a particular inter- 
mediate-frequency-amplifier design (one providing an over-all bandwidth of 
10 megacycles centered at 60 or 100 megacycles), this paper discusses qualita- 
tively the theoretical problems involved in such a design and gives data of 
practical importance to the engineer attempting to build a similar amplifier. 
Measured characteristics of approximately fifty amplifiers are summarized to 
illustrate the end results achieved. 

Historical Note on the Rate of a Moving Atomic Clock}® HERBERT E. 
Ives. The history of the idea of variation of frequency with velocity is 
followed through Goigt, Larmor, Lorentz, and Einstein. The Michelson- 
Morley experiment is explainable by any contraction of dimensions in the 
ratio (1 — v°/c?)*:1 along and transverse to the direction of motion. To 
each contraction corresponds a different value of frequency change. The 


theoretical speculations pointing to the relation am = (1 — v°/c?)’ are 
discussed, together with the significance of the experimental test by means of 
canal rays. 


New Low-Coefficient Synthetic Piezoelectric Crystals for Use in Filters and 
Oscillators! W. P. Mason. Two crystals of the monoclinic sphenoidal 
class have been found which have modes of vibration with zero temperature 
coefficients of frequency, high electromechanical coupling constants, and 
high Q’s or low dissipation. These properties make it appear probable that 
such crystals may have a considerable use in filters and oscillators as a sub- 


9 Proc. I. R. E., November 1947, 
0 Jour. Opt. Soc. Amer., October 1947. 
n Proc. I, R. E., October 1947. 
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stitute for quartz, which is difficult to obtain in large sizes. These crystals 
are ethylene diamine tartrate (EDT) having the chemical formula 
CeHisN2Og, and di-potassium tartrate (DKT) having the formula K2,C,H:O, 
— 4H,0. | 

The paper describes the properties of EDT, since this crystal has been 
found more advantageous than DKT. The 13 elastic constants, the 8 pie- 
zoelectric constants, and the 4 dielectric constants have been measured over 
a temperature range, and from these measurements the regions of low tem- 
perature coefficients and high electromechanical coupling have been located. 
Six low-temperature-coefficient cuts have been discovered and the properties 
of these cuts are given. These cuts are being applied in the crystal channel 
filters of the long-distance telephone system, and may be applied to the 
control of oscillators. 

Multi-Channel Carrier Telegraph.” A. UL. Matte. Discussion of a car- 
rier telegraph system, adapted specifically to railway requirements, to meet 
the needs for high-quality line transmission. 

Reflex.Oscillators for Radar System.2 J. O. McNAtty and W. G. SHEP- 
HERD. The advantages to be gained in the operation of radar systems at 
very high frequencies have led to the use of frequencies of several thousand 
megacycles. Operation at these frequencies has imposed serious problems 
in obtaining suitable tube behavior. Because of the difficulty in obtaining 
amplification at the transmission frequency, the r.f. section of the usual radar 
receiver consists of a crystal converter driven by a beating oscillator and 
operating directly into an if. amplifier. Since the midband frequency of 
the latter has commonly been either 30 or 60 Mc., it has been necessary to 
provide beating oscillators operating at frequencies differing from ‘those of 
the transmitter by only a few per cent. | 

For radar systems intended to operate at approximately 3000 Mc., which 
were under development in the early days of the war, it was found that tri- 
odes then available gave unsatisfactory performance. Attention shifted to 
the possibility of using velocity-modulated tubes, and the particular form 
known as the reflex oscillator came into general use. 

In this paper the requirements on beating-oscillator tubes for radar sys- 
tems are discussed, and the design features which have made the reflex 
oscillator eminently satisfactory in this application are pointed out. Prob- 
lems encountered in such oscillators,are outlined, and the solution in a num- 
ber of cases is indicated. In some instances military requirements and ex- 
pediency were in conflict with the optimum performance, and hence certain 
compromises were necessary. 


2 Railway Signaling, December 1947. 
13 Proc. I, R. E., December 1947. 
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Space-Charge and Transit-Time Effects on Signal and Noise in Microwave 
Tetrodes.* L.-C. PETERSON. Signal and noise in microwave tetrodes are 
discussed with particular emphasis on their behavior as space-charge condi- 
tions are varied in the grid-screen, or drift, region. The analysis assumes 
that the electron-stream velocity is single-valued. For particular conditions 
the noise figure may be substantially improved by increasing the space- 
charge density in the grid-screen region until an entering electron encounters 
a field of a certain magnitude. The noise reduction is largely due to the 
cancellation in the output of the noise produced by the random cathode 
emission. ‘The method of noise reduction described is applicable only when 
the transit angles of both input and drift regions are fairly long. 

In a forthcoming paper, H. V. Neher describes experimental results which 
broadly agree with the theory. 

“Cloverleaf” Antenna for F. M. Broadcasting.» Puri H. SMrre. The 
radiation requirements and general design considerations for transmitting 
antennas suitable for f.m. broadcoasting are briefly discussed, and an ex- 
planation of the design and operation of the arrangement of radiating ele- 
ments and associated feed system employed in the “cloverleaf” antenna is 
given. Both calculated and measured data are included, showing field-in- 
tensity distribution, gain, impedance-frequency characteristics, etc. De- 
sign features which are discussed include a simple coaxial impedance-match- 
ing transformer developed initially for microwave application, and the 
method and facilities provided for the removal of sleet. 

Hybrid Circuits for Microwaves.6 W. A. TyRRELL. The fundamental 
behavior of hybrid circuits is reviewed and discussed, largely in terms of 
reciprocity relationships. The phase properties of simple wave-guide tee 
junctions are briefly considered. Two kinds of hybrid circuits are then 
described, the one involving a ring or loop of transmission line, the other 
relying upon the symmetry properties of certain four-arm junctions. The 
description is centered about wave-guide structures for microwaves, but the 
principles may also be applied to other kinds of transmission lines for other 
frequency ranges. Experimental verification is provided, and some of the 
important applications are outlined. 

- M4 Proc. I. R. E., November 1947, 


_ & Proc. I. R. E-Waves and Electrons Section, December 1947. 
' 16 Proc, I. R. E., November 1947, 
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